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PREFACE 


The  work  presented  in  this  report  was  performed  by  Kaman  Aerospace  Corpor¬ 
ation  under  Contract  DAAK51 -78-C-001 7  for  the  Applied  Technology  Laboratory, 
U.  S.  Army  Research  and  Technology  Laboratories  (AVRADCOM),  Fort  Eustis, 
Virginia.  The  program  was  implemented  under  the  technical  direction  of 
Nicholas  J.  Calapodas  of  the  Aeronautical  Technology  Division. 

The  report  presents  a  general  description  of  the  analytical  and  experimental 
techniques  for  implementing  system  identification  methodology  with  specific 
application  to  the  AH-16  airframe. 


SUMMARY 


The  purpose  of  this  program  was  to  perform  experimental  and  analytical 
work  to  verify  the  structural  system  identification  methodology  developed 
at  Kaman  Aerospace  Corporation.  To  accomplish  this,  system  Identification 
techniques  were  applied  to  a  U.  S.  Army  AH-1G  airframe  (S/N  57-15683)  to 
create  mathematical  models  from  dynamic  test  data  and  to  improve  a  reduced 
version  of  the  existing  NASTRAN  model  of  the  AH-1G  helicopter  using  dynamic 
test  data.  An  additional  objective  of  the  program  was  to  test  the  effect¬ 
iveness  of  the  mathematical  models  developed  in  predicting  the  effects  of 
structural  stiffness  and  mass  changes  to  the  helicopter. 

The  analytical  models  that  were  developed  proved  to  be  effective  in  pre¬ 
dicting  the  effects  of  the  stiffness  and  mass  changes  which  were  incorpor¬ 
ated  in  the  helicopter.  Additionally,  the  improved  version  of  the  reduced 
NASTRAN  model  successfully  predicted  the  frequency  response  of  the  AH-1G 
airframe  for  the  frequency  range  of  interest. 

A  significant  part  of  the  program  involved  development  and  refinement  of 
test  procedures  for  acquiring  impedance  type  test  data  which  would  yield 
modal  data  suitable  for  implementation  of  the  system  identification  method¬ 
ology.  Complex  modes  resulting  from  elastomeric  isolators  and  nonstruc- 
tural  access  panels  were  found  to  exist  in  the  AH-1G  airframe.  These  com¬ 
plex  modes  were  successful ly  measured,  which  represents  an  advancement 
in  the  state  of  the  art  of  dynamic  testing  and  modal  parameter  extraction. 
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INTRODUCTION 


High  vibration  levels  experienced  by  current  production  helicopters  result 
in  fatigue  damage  to  airframe  structure  and  components  with  ultimate  re¬ 
duction  in  reliability.  Additionally,  high  vibratory  structural  response 
and  accompanying  high  noise  levels  in  the  cabin  area  degrade  the  efficien¬ 
cy,  health,  and  comfort  of  both  the  crew  and  passengers.  Currently,  the 

solution  to  this  problem  relies  on  flight  test  information  to  effect  vibra¬ 
tion  reduction  procedures,  but  these  procedures  result  in  high  cost  and 
schedule  delays.  An  alternate  solution  is  to  develop  analytical  methods  to 
accurately  predict  the  vibration  response  characteristics  of  the  helicopter 
and  allow  structural  modifications  during  prototype  development  to  reduce 
vibration  levels  with  inherent  lower  development  time  and  cost. 

Thus,  the  Army  has  a  continuing'  goal  to  improve  analytical  methods  capable 
of  predicting  the  vibration  levels  experienced  by  a  helicopter  in  flight. 
One  aspect  of  this  involves  the  analytical  modeling  of  rotorcraft  airframes 
and,  in  particular,  use  of  the  finite-element  method  to  make  analytical 
predictions  of  airframe  vibration  frequencies  and  responses.  The  finite- 
element  method,  which  is  the  only  method  in  general  use  capable  of  the 
detailed  mathematical  modeling  of  helicopter  structures,  has  experienced 
difficulty  in  predicting  airframe  vibration.  Consequently,  improvement 
in  the  modeling  of  helicopter  structures  for  the  purpose  of  vibration 
prediction  is  essential. 

Once  a  prototype  has  been  built,  system  identification  methodology  can  be 
used  to  generate  mathematical  models  or  to  improve  existing  analytical 
models  using  ground  vibration  test  data.  Thus,  system  identification  can 
be  applied  to  combine  test  data  and  analytical  models  into  an  engineering 
tool  useful  in  developing  new  rotorcraft  and  in  incorporating  modifica¬ 
tions  to  existing  airframes.  System  identification  unifies  analyses  de¬ 
veloped  from  idealizations  regarding  the  mass,  stiffness,  and  damping 
characteristics  of  the  structure  with  dynamic  test  methods  resulting  in 
improved  mathematical  models  which  can  successfully  predict  the  dynamic 
response  characteristics  of  the  helicopter. 
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Kaman  Aerospace  Corporation  has  developed  system  identification  methodology 
through  two  separate  techniques.  The  Incomplete  Model  and  the  Truncated 
Model  methods  have  progressed  successfully  from  theoretical  formulation, 
computer  simulations,  application  to  small  structures,  and  finally  to  full- 
scale  craft  applications. 

The  purpose  of  the  present  program  was  to  apply  system  identification  method 
ology  to  the  AH-1G  airframe  to  create  a  mathematical  model,  to  improve  the 
existing  NASTRAN  model  using  dynamic  test  data,  and  to  test  the  effective¬ 
ness  of  the  mathematical  models  in  predicting  the  effects  of  mass  and  stiff¬ 
ness  changes  to  the  airframe.  The  program  involved  extensive  shake  tests  of 
several  mass  and  stiffness  configurations  and  extraction  of  natural  frequen¬ 
cies,  orthonormal  mode  shapes,  and  modal  parameters  for  the  respective 
conditions. 

Ground  vibration  tests  of  the  AH-1G  helicopter  revealed  that  the  airframe 
is  characterized  by  complex  mode  shapes  which  are  caused  by  the  nonpropor¬ 
tional  distribution  of  damping  resulting  from  the  elastomeric  isolators  and 
nonstructural  access  panels  in  the  structure.  These  complex  modes  were 
successfully  measured  during  the  test  phase  of  the  program;  nowever,  both 
the  Incomplete  Model  and  the  Truncated  Model ,  which  rely  on  mode  shapes  in 
their  formulation,  cannot,  currently,  accommodate  the  complex  modes.  Conse¬ 
quently,  the  complex  modes  had  to  be  transformed  to  the  real  domain  for  use 
in  the  theories.  It  is  essential  that  the  Incomplete  Model  and  Truncated 
Model  methodologies  be  extended  to  include  the  complex  mode  feature  to 
improve  the  accuracy  of  their  vibration  prediction  capability. 
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NASTRAN  MODEL 


EXISTING  NASTRAN  MODEL  OF  AH-IG 


An  important  aspect  of  helicopter  vibration  prediction  is  related  to  the 
analytical  modeling  of  the  airframe,  and,  in  particular,  the  use  of  the 
finite  element  method  to  make  analytical  predictions  of  the  dynamic  re¬ 
sponse  characteristics  of  the  vehicle.  The  finite  element  method  is  the 
only  method  in  general  use  that  is  capable  of  the  detailed  mathematical 
modeling  required  for  the  analysis  of  the  complex,  redundant  airframe 
structures  of  modern  helicopters.  The  basic  NASTRAN  model  of  the  AH-IG 
helicopter  used  in  the  current  program  is  documented,  in  detail,  in  Refer¬ 
ence  1.  The  NASTRAN  model  was  developed  to  represent  the  vibration  re¬ 
sponse  of  the  AH-IG  airframe  for  excitation  frequencies  below  30  Hz.  Fig¬ 
ure  1  illustrates  the  actual  airframe  structure  with  skins  removed  and 
Figure  2  presents  the  NASTRAN  model,  which  consists  of  a  built-up  ideal¬ 
ization  of  the  fuselage  and  an  elastic  line  representation  for  the  tail 
boom,  vertical  fin,  and  horizontal  stabilizer.  As  described  in  Reference 
1,  several  regions  of  the  structure  are  not  effective  in  providing  stiff¬ 
ness  and  are  not  included  in  the  NASTRAN  model.  These  components  are  the 
canopy,  cowling  around  the  engine  and  main  rotor  pylon,  access  panels  at 
fuselage  station  (FS)  61.25  in.  to  213.94  in.,  doors  on  the  ammunition 
bay  at  FS  95  in.  to  135  in.,  top  access  door  on  the  nose  at  FS  33  in.  to 
43  in.,  and  the  drive  shaft  connection  of  the  engine  to  the  main  trans¬ 
mission.  The  model  was  intuitively  developed  using  structural  elements 
available  in  the  NASTRAN  library  and  consists  of  bars,  rods,  shear  panels, 
triangular  and  quadralateral  plate  elements,  and  scalar  spring  elements. 


T. Cronkhite,  James  D.,  et  al.  A  NASTRAN  VIBRATION  MODEL  OF  THE  AH-IG 
HELICOPTER  AIRFRAME,  Vol .  1,  USAAMRDL  R-TR-74-045,  Rock  Island 
Arsenal,  U.  S.  Army  Air  Mobility  Research  and  Development  Laboratory, 
Rock  Island,  Illinois,  June  1974,  AD-A009482. 


15 


NASTRAN  MQDF.L  USED  IN  SYSTEM  IDENTIFICATION 


Originally,  the  NASTRAN  model  contained  2940  degrees  of  freedom  and  was 
reduced  to  2699  degrees  of  freedom  through  the  use  of  Multi-Point  Constraint 
(MPC)  equations.  Use  of  the  MPC  equations  allows  description  of  a  degree 
of  freedom  as  a  linear  combination  of  the  values  of  two  or  more  degrees  of 
freedom.  The  model  was  further  reduced  to  1714  degrees  of  freedom  by  the 
application  of  Single  Point  Constraints  (SPC),  which  forces  a  degree  of 
freedom  to  zero.  Use  of  the  GUYAN  reduction,  which  is  a  static  condensation 
procedure,  further  deleted  the  degrees  of  freedom  to  98,  and  eliminating 
six  support  degrees  of  freedom  reduced  the  model  to  the  final  92  degrees 
of  freedom. 

The  existing  NASTRAN  model  of  the  AH-1G  was  reduced  to  92  degrees  of  free¬ 
dom  for  the  purpose  of  the  Structural  System  Identification  Technology  Veri¬ 
fication  Program.  The  92  degrees  of  freedom  selected  were  consistent  with 
the  locations  of  the  accelerometers  placed  on  the  airframe  as  described 
in  Reference  2  and  consistent  with  the  requirements  of  the  subject  con¬ 
tract. 

GUYAN  Reduction 


The  mathematics  of  the  static  condensation  procedure  that  yields  the  re¬ 
duced  mass  and  stiffness  matrices  in  a  NASTRAN  dynamic  analysis  is  briefly 
described.  Since  the  technique  is  a  static  condensation  it  is  theoretically 
applicable  only  at  zero  excitation  frequency.  However,  it  is  generally 
used  independent  of  excitation  frequency.  The  formulation  is  similar  to 


2.  Jones,  R. ,  et  al.,  EXPERIMENTAL  VERIFICATION  OF  FORCE  DETERMINATION 
AND  GROUND  FLYING  ON  A  FULL-SCALE  HELICOPTER.  Kaman  Aerospace  Cor¬ 
poration  Report  Number  R-1625,  USAAVRADCOM-TR-81 -D-l 1 ,  Applied  Tech¬ 
nology  Laboratory,  U.  S.  Army  Research  and  Technology  Laboratory, 
Fort  Eustis,  Virginia,  23604,  May  1981,  AD  A100182,. 
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that  of  Reference  3  and  is  presented  hare  only  for  completeness.  The 
equilibrium  equations  for  a  structure  may  be  written  as 


[K]  (x)  -  (F)  (1) 

The  structural  equations  can  be  rearranged  and  written  In  partitioned 
form  as 


t  \  r 
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K22  . 

>  *2 

■rf 

Assuming  that  no  forces  are  applied  at  the  coordinate  set  x2  the  two 
resulting  matrix  equations  developed  from  equation  (2)  become 

K11  X1  +  K12  x2  *  F1  ^ 


+  K22  x2  =  0 


(4) 


The  coordinate  set  x2  can  be  obtained  from  equation  (4) 


x2 


x 


1 


(5) 


Substituting  equation  (5)  into  equation  (3)  will  eliminate  the  coor¬ 
dinates  x2 

-  ICj2  K.-,  ^12  ^  X1  =  F1  (®) 

where  [K]^  =  J_^n  ~  ^12  ^22~^  *^12^]  defined  as  the  reduced  stiff¬ 
ness  matrix.  The  original  coordinate  system  has  effectively  been 


3.  Guyan,  R. ,  Reduction  of  Stiffness  and  Mass  Matrices.  AIAA  JOURNAL, 
vol.  3,  No.  2,  February  1965,  pp.  380. 
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transformed  to  the  coordinate  system  at  which  forces  are  considered  to  act 
thus  /  \  h*  , 


{»}  ■  h  ■)  ■ F-  -  -  -  ]  {*,} 

*•  J  x,  v  T  'J 

l  c  ‘*22  *12 


[T]  {x,} 


The  kinetic  energy  can  be  written  as 

V  =  ?  {*}  [Ml  {*} 

Substituting  equation  (7)  Into  equation  (8)  yields 

V  *  l  {x,}  T  [T]T  [M]  CT]  {x,} 
The  reduced  mass  matrix  Is  then  defined  as 


[m]red  *  [t]'  CM]  [T]  do; 

Substituting  for  [T]  from  equation  (7)  and  writing  equation  (10)  in  parti¬ 
tioned  form,  the  result  is 


*^RED  T  !  “  K22  K12 


-i  „  t  n i 


M12  I 


-1  ..  T 


and  finally 


The  reduced  mass  and  stiffness  matrices  developed  In  the  Guyan  Reduc- 
tior  procedure  are  nondlagonal. 

NASTRAN  Model  Gross  Weight  Configuration 

The  NASTRAN  model  of  the  AH-1G  as  described  In  Reference  1  was  developed 
for  a  production  configuration  Identified  as  FY71  S/N  21123.  The  AH-1G 
S/N  57-15683  was  used  In  the  Structural  System  Identification  Tech¬ 
nology  Verification  Program.  As  a  result,  the  existing  Government-furnish¬ 
ed  NASTRAN  model  of  the  AH-1G  had  to  be  modified  for  compatibility  with 
the  aircraft  tested.  This  required  modifying  large  weight  items  and  useful 
load  items.  The  gross  weight  configuration  of  the  original  NASTRAN  model 
was  8933  lb,  consisting  of  a  total  weight  empty  of  5760  lb  and  a  useful 
load  of  3*173  lb.  In  addition,  the  main  rotor  blade  weight  and  main  rotor 
hub  weight  were  included  in  the  model  as  concentrated  weights.  The  NASTRAN 
model  of  the  aircraft  that  was  shake  tested  during  the  subject  contract 
was  8750  lb;  the  main  rotor  blades  were  omitted  and  only  231  lb  of  hub 
weight  were  included.  The  following  tabulations  describe  the  concentrated 
mass  items  that  were  changed  from  the  original  NASTRAN  AH-1G  model. 


LB 

Original  NASTRAN  AH-1G  Model 

8933 

Weight  deleted 

XM-28  Gun  Turret  Launcher 

41 

Mini -gun 

48 

40-mm  Drum/Pa llet 

75 

40-mm  Grenades  (250  rounds) 

190 

7.62-mm  Mamee  Drum 

63 

7.62-mm  Ammunition  (4000  rounds) 

260 
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Forward  Fuel  780 

Aft  Fuel  830 

Pilot  200 

Gunner  200 

Main  Rotor  Hub  490 

Main  Rotor  Blades  458 

XM-157A  Pod  (Right  Wing)  57 

M-151/XM  429  Rockets  (Right  Wing)  146 

XM-157A  Pod  (Left  Wing)  57 

M-151/XM  439  Rockets  (Left  Wing)  146 

Smoke  Grenades  and  Dispenser  (Right  Wing)  40 

Smoke  Grenades  and  Dispenser  (Left  Wing)  40 

mr 

NASTRAN  model  of  AH-1G  S/N  57-15683 
Weight  added _ 

Test  .’instrumentation  195 

M2 00  A/1  Inboard  Rocket  Pod/19  Rockets  (Right  Wing)  529 

Inboard  Pylon  Assembly  (Right  Wing)  31 

M2 00  A/1  Inboard  Rocket  Pod/19  Rockets  (Left  Wing)  529 

Inboard  Pylon  Assembly  (Left  Wing)  31 

M2 00  A/1  Outboard  Rocket  Pod/12  Rockets  (Right  Wing)  385 

M2 00  A/1  Outboard  Rocket  Pod/12  Rockets  (Left  Wing)  385 

Main  Rotor  Hub  231 

Turret  Ballast  150 

Gun  Sight  Ballast  25 

Pilot  230 

Gunner  180 

Fuel  933 

Instrumentation  (accelerometers, wire, etc. )  TOO 


ms 


Thus,  the  NASTRAN  model  weight  for  the  configuration  tested  was 
WT  =  8933  -  4121  +  3934  =  8746  lb 

This  compares  with  the  actual  weight  of  the  test  vehicle  of  8784  lb 
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Coordinate  Selection 


Selection  of  coordinates  on  the  AH-1G  fuselage  was  based  on  require¬ 
ments  for  determining  the  modal  acceleration  coefficients  of  the  struc¬ 
ture.  In  NASTRAN,  the  Guyan  reduction  Is  the  technique  used  for  reduc¬ 
ing  the  number  of  static  degrees  of  freedom  to  be  used  In  a  dynamic 
analysis.  It  Implies  a  meaningful  redistribution  of  Inertia  properties 
to  a  subset  of  the  total  number  of  degrees  of  freedom  required  In  a 
static  analysis  of  the  structure.  One  of  the  primary  requirements  In 
selection  of  the  dynamic  analysis  set  of  coordinates  Is  that  they  should 
be  uniformly  distributed  throughout  the  structure.  This  Insures  that 
the  modal  acceleration  coefficients  and,  subsequently,  the  orthonormal 
modes  will  be  representative  of  the  deformation  patterns  of  the  structure. 

Ideally,  the  transformation  from  the  full  size  model  to  the  smaller  size 
required  for  an  operational  model  can  be  achieved  with  the  reduced  model 
capable  of  accurate  analytical  prediction  of  the  dynamic  response 
characteristics  of  the  vehicle.  Practically,  however,  the  reduction  pro¬ 
cess  introduces  error  in  the  smaller  size  model.  Minimizing  the  error 
Inherent  In  the  reduction  procedure  Is  of  paramount  Importance.  Guidelines 
for  developing  the  reduced  model  through  Implementation  of  the  Guyan  Re¬ 
duction  procedure  are  described  in  Reference  1.  These  principles  were 
carefully  followed  in  selecting  degrees  of  freedom  pertinent  to  the  AH-1G 
reduced  model.  A  total  of  92  degrees  of  freedom  were  used  in  the  reduced 
NASTRAN  model  used  in  the  study.  The  government-furnished  NASTRAN  model 
consisting  of  2940  degrees  of  freedom  was  reduced  to  92  degrees  of 
freedom  using  SPC,  MPC,  and  the  Guyan  technique. 

Accelerometer  Location 


A  total  of  92  accelerometers  were  available  on  the  AH-1G  fuselage  for  use 
In  the  current  program.  Of  this  total,  58  accelerometers  were  positioned 
at  locations  pertinent  to  the  system  identification  program  and  34  accel¬ 
erometers  had  been  installed  for  previous  use  as  described  in  Reference  2. 
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Generally,  the  locations  for  placement  of  the  accelerometers  on  the  AH-1G 
fuselage  were  selected  within  the  constraints  specified  In  Reference  1. 
Ideally,  the  fuselage  station,  butt  line  station,  and  water  line  station 
of  each  accelerometer  should  coincide  with  a  particular  grid  point  in  the 
NASTRAN  model.  However,  since  the  accelerometers  were  rigidly  attached 
to  the  structure  via  mounting  T-locks  to  minimize  extraneous  localized 
response  effects  it  was  not  possible,  In  every  Instance,  to  locate  the 
accelerometers  coincident  with  a  NASTRAN  grid  point.  The  accelerometers 
were  placed  with  direction  of  major  sensitivity  oriented  appropriately 
with  respect  to  the  fuselage  station,  butt  line,  and  water  line  orthogonal 
axes  system.  In  addition,  accelerometers  oriented  with  major  axis  In  the 
water  line  direction  but  offset  from  the  butt  line  zero  axis  could  provide 
a  measure  of  the  torsional  deformation  of  the  fuselage.  The  positive 
directions  in  the  axes  system  were  aft,  up  and  to  the  right,  as  viewed  from 
the  pilot  seat.  Table  1  presents  a  listing  of  the  accelerometer  designa¬ 
tion  and  location  on  the  AH-1G  in  terms  of  fuselage,  butt  line,  and  water 
line  stations. 

The  accelerometer  identifications  are  prefixed  by  either  the  letters  X, 

Y  or  Z,  which  refer  to  the  direction  of  acceleration  measurement  in  the 
longitudinal,  lateral  or  vertical  directions,  respectively.  Additionally, 
the  letters  R,  L,  T  and  B  are  used  as  suffixes  to  describe  the  location 
of  the  accelerometer,  with  R  and  L  referring  to  right  and  left  distances 
from  the  butt  line  zero  reference,  and  T  and  B  referring  to  top  and  bottom 
surface  of  the  fuselage.  As  an  illustration  of  the  accelerometer  identi¬ 
fication,  Z46R  refers  to  an  accelerometer  that  is  oriented  for  measurement 
of  acceleration  along  the  vertical  axis,  the  number  46  refers  to  the  nominal 
fuselage  location,  and  the  suffix  R  indicates  the  right  side  or  positive 
lateral  axis  direction.  The  actual  position  of  accelerometer  Z46R  is  given 
in  the  table  as  FS  46  in.,  BL  9.5  in.,  and  WL  65.5  in.  Figure  3  presents 
geometric  data  and  the  reference  axes  system  for  the  AH-1G  test  program. 
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TABLE  1.  ACCELEROMETER  DESCRIPTION. 


No. 

Accelerometer 

Identification 

Fuselage 

Station 

Butt  Line 
Station 

Water  Line 
Station 

1 

Y46L 

46.00 

-10.10 

45.00 

2 

Z46R 

46.00 

9.50 

65.50 

3 

Y50 

40.25 

0.00 

52.00 

4 

Z50 

38.75 

0.00 

52.00 

5 

X60R 

59.75 

9.50 

51.50 

6 

Y60L 

60.25 

-  9.50 

59.50 

7 

Z60L 

60.00 

-17.50 

60.00 

8 

Z61R 

61.25 

17.50 

60.00 

3 

Y90 

103.00 

0. 00 

52.00 

10 

Z90R 

100.75 

20.00 

50.00 

11 

Z90L 

100.75 

-20.00 

50.00 

12 

Z92L 

92.00 

-18.00 

46.00 

13 

Z93R 

93.00 

19.00 

46.00 

14 

Y95R 

95.00 

9.50 

46.50 

15 

Z100T 

102.75 

0.00 

100.00 

16 

Y115R 

115.00 

10.00 

45.00 

17 

Z115R 

115.00 

8.75 

62.00 

18 

Y139R 

139.00 

18.00 

46.00 

19 

X140 

138.00 

0.00 

50.00 

20 

Y140 

138.00 

0.00 

50.00 

21 

Z140R 

133.50 

20.00 

50.00 

22 

Z140L 

134.50 

-20.00 

50.00 

23 

Z147R 

147.00 

8.70 

55.50 

24 

Y152L 

152.45 

-12.50 

27.50 

25 

Z155R 

155.50 

12.00 

28.50 

26 

Y156 

156.41 

0.00 

35.00 

27 

Z156R 

156.41 

19.00 

36.00 

28 

Z156L 

156.41 

-19.00 

36.00 

29 

X180T 

209.00 

-  6.50 

125.00 

30 

Y187 

187.50 

0.00 

46.00 

31 

X188 

187.75 

0.00 

56.95 

D.O.F. 

No. 
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ass 


No. 

Accel erometer 
Identification 

Fuselage 

Station 

Butt  Line 
Station 

Water  Line 
Station 

D.O.F. 

No. 

32 

Z189 

189.25 

0.00 

36.50 

64 

33 

Z195B 

195.10 

0.00 

65.00 

37 

34 

Z195T 

195.10 

0.00 

75.00 

25 

35 

ZLONG 

196.00 

12.00 

59.00 

36 

36 

ZLATR 

196.00 

-12.00 

59.00 

38 

37 

ZCOLL 

207.00 

-12.00 

57.00 

34 

38 

X200R 

193.50 

55.00 

65.00 

55 

39 

X200L 

193.50 

-55.00 

65.00 

59 

40 

Z200R 

193.50 

55.00 

65.00 

56 

41 

Z200L 

193.50 

-55.00 

65.00 

60 

42 

Z201R 

201.10 

42.50 

65.00 

54 

43 

Z201L 

201.10 

-42.50 

65.00 

58 

44 

Z202R 

202.00 

60.00 

63.10 

E7 

J  1 

45 

Z202L 

202.00 

-60.00 

63.10 

61 

46 

XMHUB 

197.00 

0.00 

157.05 

29 

47 

YHUB 

198.38 

2.00 

157.05 

30 

48 

ZHUB 

198.38 

0.00 

158.18 

31 

49 

MZHUB 

201 . 50 

0.00 

149.55 

32 

50 

Z210T 

207.00 

-  7.00 

125.00 

28 

51 

Z215R 

214.70 

18.59 

55.58 

52 

52 

Z215L 

214.70 

-18.59 

55.58 

53 

53 

X220L 

212.00 

-20.00 

35.00 

51 

54 

Y220B 

212.00 

0.00 

35.00 

50 

55 

Y220T 

208.00 

-  8.00 

125.00 

27 

56 

Y226R 

226.41 

10.18 

29.15 

46 

57 

Z226 

226.40 

0.00 

35.12 

48 

58 

Z226L 

226.41 

-10.67 

29.09 

49 

59 

Y229R 

229.12 

7.86 

26.00 

47 

60 

Z250 

250.00 

0.00 

64.07 

35 

61 

Y251 

251.37 

0.00 

52.59 

39 

62 

Z251 

251.00 

0.00 

79.50 

33 
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TABLE  1.  CONTINUED. 


No. 

Accelerometer 

Identification 

Fuselage 

Station 

Butt  Line 
Station 

Water  Line 
Station 

D.O.F 

No. 

63 

Y267 

267.25 

1.00 

64.50 

41 

64 

Z268R 

268.25 

12.00 

64.50 

40 

65 

Z268L 

268.25 

-12.00 

64.20 

42 

66 

Y300 

300. 00 

6.00 

50.00 

43 

67 

X301 

300.70 

-  0.50 

30.75 

44 

68 

Z301 

301.20 

0.00 

30.00 

45 

69 

Y339L 

338.61 

-13.00 

57.02 

6 

70 

Z340 

341.50 

0.00 

42.00 

7 

71 

Y380 

379.00 

0.00 

50.00 

4 

72 

Z381 

381.09 

0.00 

42.25 

5 

73 

Z396R 

396.50 

21.50 

60.00 

2 

74 

Z396L 

396.50 

-21.50 

60.00 

3 

75 

Z400 

400.00 

0.00 

48.00 

9 

76 

Z404R 

404.00 

31.25 

56.50 

1 

77 

X405 

405.33 

-  0.50 

44.75 

8 

78 

Y422L 

422.24 

-  9.00 

57.21 

10 

79 

Z422L 

422.24 

-  8.50 

57.71 

12 

80 

Z423 

423.08 

0.00 

46.25 

11 

81 

Y440 

448.00 

0.00 

50.00 

13 

82 

Z460 

465.40 

0.00 

65.00 

18 

83 

Y464L 

464.10 

-  7  00 

59.81 

14 

84 

Z464 

464.10 

0.00 

51.75 

15 

85 

X479 

479.23 

0.00 

54.25 

16 

86 

Y480 

480. 23 

-  4.00 

72.40 

17 

87 

Z485 

485.50 

0.00 

93.25 

20 

88 

Y490 

489.00 

0.00 

75.00 

19 

89 

Y517 

517.50 

-  5.50 

105.25 

21 

90 

Y521 

520.67 

-  4.00 

118.27 

23 

91 

X540 

545.00 

0.00 

135.00 

22 

92 

Z540 

545.00 

0.00 

135.00 

24 

27 
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MODE  SHAPES 
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CLASSICAL  MODES 

Classical  modes,  frequently  referred  to  as  normal  modes,  exist  in  systems 
which  have  proportional  damping  of  any  magnitude.  This  implies  that  the 
mode  shapes  for  the  damped  system  are  the  same  as  those  for  the  system 
with  zero  damping.  The  motion  in  any  mode  is  characterized  by  the  exist¬ 
ence  of  fixed  nodes  on  the  structure,  and  all  points  in  the  system  attain 
their  maximum  or  minimum  response  simultaneously.  Therefore  the  elements 
of  each  node  will  either  have  the  same  phase  or  be  180  degre  e  out  of 
phase,  describing  a  standing  wave  pattern. 

In  addition,  if  the  modes  are  reasonably  separated  in  frequency,  the  real 
part  of  the  acceleration  mobility  will  display  two  turning  points  in  the 
vicinity  of  a  natural  frequency  and  the  imaginary  part  of  the  acceleration 
mobility  will  show  a  single  turning  point  at  the  natural  frequency. 

Figure  4a  presents  acceleration  mobility  measurements  obtained  on  the  AH-1G 
airframe.  The  figure  shows  a  transfer  function  relating  vertical  accelera¬ 
tion  measured  at  a  point  in  the  nose  of  the  airframe  (Z50)  for  excitation 
applied  vertically  at  the  tail  skid  (Z485).  Two  close,  but  distinguish¬ 
able,  modes  are  present.  The  dominant  mode  is  essentially  classical,  dis¬ 
playing  two  turning  points  in  the  real  acceleration  mobility  and  a  single 
turning  point  in  the  imaginary  acceleration  mobility  measurements. 

Figure  4b  shows  the  acceleration  mobility  data  in  the  complex  plane  common¬ 
ly  referred  to  as  the  Nyquist  diagram. 

COMPLEX  MODES 


Complex  modes  exist  in  a  structure  if  it  is  characterized  by  the  presence 
of  nonproportional  damping.  In  the  classical  situation,  the  modes  are 
distinguished  by  a  shape  wherein  all  points  in  the  system  attain  their 
maximum  or  minimum  motion  at  the  same  time.  In  contrast,  a  system  with 
nonproportional  damping,  and  therefore  complex  modes,  is  not  characterized 
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a.  Measured  acceleration  mobility. 
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b.  Nyquist  diagram  (frequency  5. 5-9.0  Hz) 
Figure  4.  Acceleration  Mobility  of  the  AH-1G. 
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by  the  presence  of  fixed  nodes  In  each  mode  of  vibration.  Rather*  the 
points  on  the  structure  undergo  motion  at  the  same  frequency,  but  with 
different  phase  relationships.  Hence,  complex  modes  can  be  described 
as  traveling  waves  as  opposed  to  standing  waves,  which  distinguish  the 
normal  modes  of  vibration. 

For  the  general  case  of  nonclasslcal  or  complex  modes  the  real  and  imagin¬ 
ary  acceleration  mobility  may  not  exhibit  the  characteristic  shape  displayed 
In  the  classical  mode  situation.  It  is  possible  for  the  Imaginary  mobility 
to  display  two  turning  points  in  the  region  of  a  natural  frequency  while 
the  real  mobility  may  show  a  single  turning  point  at  a  natural  frequency. 
Figure  5a  shows  an  example  of  this  phenomenon  in  data  measured  on  the  AH-1G 
airframe.  The  data  presented  are  for  a  transfer  function  describing  the 
acceleration  response  measured  at  the  nose  of  the  airframe  for  force  excita¬ 
tion  applied  vertically  at  the  tail  skid.  It  is  also  possible  that  either 
or  both  of  the  real  and  the  imaginary  parts  of  the  acceleration  mobility 
will  display  two  turning  points  proximate  to  a  natural  frequency. 

Figure  5b  shows  che  aforementioned  measured  acceleration  mobility  in  the 
complex  plane. 

Complex  modes  are  likely  to  exist  in  a  helicopter  airframe  because  of  the 
nonproportional  distribution  of  damping.  In  the  AH-1G  airframe,  elastomeric 
isolators  and  nonstructural  access  panels  are  the  source  of  nonproportional 
damping  resulting  in  the  existence  of  complex  modes. 

MAJOR  MODES 


In  evaluating  the  overall  dynamic  response  of  a  system,  the  relative  contri¬ 
bution  of  each  of  the  modes  in  a  given  frequency  range  of  interest  is  of 
paramount  importance.  The  problem  of  identifying  the  significant  natural 
frequencies  of  a  structure  is  related  to  distinguishing  the  difference 
between  what  can  be  considered  to  be  a  major  mode  and  what  can  be  defined 
as  a  minor  or  local  mode.  At  a  particular  excitation  frequency,  it  is 
possible  to  excite  a  local  or  minor  elastic  mode  which  has  negligible 
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b.  Nyquist  diagram  (frequency  40-50  Hz) 
Figure  5.  Acceleration  Mobility  of  the  AH-1G. 


effect  on  the  overall  system  dynamic  response.  As  an  example,  consider 
the  flexural  vibrations  of  the  tall  skid  on  the  AH-1G  helicopter.  Accelera¬ 
tion  mobility  measurements  on  that  structure  will  show  the  natural  frequen¬ 
cies,  damping  coefficients  and  mode  shapes  of  all  the  elastic  modes  invol¬ 
ving  the  motion  of  the  tail  skid.  Obviously,  these  modes  are  not  relevant 
to  the  dynamic  behavior  of  the  airframe.  However,  with  only  normal  mode 
shapes,  natural  frequencies, and  structural  damping  coefficients  obtained 
from  the  test  data,  the  question  of  modal  significance  cannot  be  answered. 
If,  however,  the  normal  mode  is  combined  with  the  generalized  mass  the 
result  is  the  orthonormal  mode  which  has  physical  units.  Modes  of  physi¬ 
cal  units  are  absolutely  necessary  in  distinguishing  major  modes  from 
local  modes.  Additionally,  the  Euclidean  norm  of  each  mode  ik  is  a 
quantity  that  has  physical  significance  and  indicates  the  relative  contri¬ 
bution  of  each  mode  to  the  total  structural  dynamic  response. 

Local  or  Residual  Modes 


Certain  points  on  an  airframe  exhibit  response  characteristics,  as 
evidenced  by  their  modal  acceleration  properties,  which  may  be  signifi¬ 
cantly  higher  in  any  mode  than  those  of  other  points  on  the  structure. 

This  effect  defines  a  local  mode  since  only  geometrically  proximate  ortho¬ 
normal  modal  elements  are  extraordinarily  large  in  the  particular  mode. 

This  effect  is  most  critical  at  a  driving  point  although  it  may  also 
occur  in  transfer  mobilities.  The  result  is  that  high  frequency  local 
modes  can  have  a  pronounced  effect  on  mobilities  at  much  lower  frequen¬ 
cies.  During  dynamic  testing  of  the  AH-1G,  at  Kaman  Aerospace  Corporation, 
the  discovery  was  made  that  high  frequency  local  modes,  in  the  400-Hz 
range,  could  severely  influence  the  low  frequency  response  of  the  air¬ 
frame.  This  phenomenon  was  observed  when  force  excitation  was  applied 
at  low  frequency  for  driving  point  measurements.  It  occurred  during 
AH-1G  dynamic  tests  for  vertical  excitation  applied  proximate  to  the 
tail  skid  attachment  point  and  for  vertical  excitation  applied  on  the 
airframe  in  the  vicinity  of  the  horizontal  tail.  This  effect  was  not 
witnessed,  however,  for  hub  shaking  or  for  lateral  force  excitation 
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applied  at  the  tall  rotor  gearbox.  Although  a  common  happening,  this 
local  mode  effect  does  not  occur  at  all  driving  points. 

Figure  5  shows  the  real  and  Imaginary  parts  of  the  driving-point  accelera¬ 
tion  mobility  measured  at  the  tall  skid  attachment  point  of  the  AH-1G  air¬ 
frame  for  a  frequency  range  of  0  to  100  Hz.  It  Is  evident  that  the 
acceleration  mobility  Is  steadily  drifting  away  from  the  zero  reference 
mobility  axis.  Apparently,  the  bounds  of  a  local  mode  have  not  been  spanned. 
Indeed,  when  the  frequency  spectrum  Is  expanded  to  800  Hz,  as  shown  in 
Figure  7,  the  residual  modes  occurring  In  the  400-Hz  range  are  clearly 
evident.  The  modal  coefficients  associated  with  these  modes  are  large, 
and  If  not  properly  removed  analytically  can  contribute  significantly  to 
the  acceleration  mobility  measured  at  lower  frequencies. 

The  modal  parameter  extraction  techniques  that  were  employed  In  the  analy¬ 
sis  require  mobility  data  acquired  at  discrete  frequencies  on  either  side 
of  a  natural  frequency.  Since  the  methods  operate  on  a  mobility  differ¬ 
ence  proximate  to  the  natural  frequency,  the  local  mode  or  residual  effect 
is  minimized. 
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Figure  6.  Tail  vertical  driving  point  acceleration 
mobility  (frequency  0-100  Hz). 
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Figure  7.  Tail  vertical  driving  point  acceleration 
mobility  (frequency  0-800  Hz). 
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MODAL  TRANSFORMATION 


COMPLEX  MODES  TO  REAL  MODES 


The  AH-1G  airframe  Is  characterized  by  nonproportlonal  damping;  as  a 
result,  It  exhibits  complex  modes.  The  orthonormal  mode  shapes  extracted 
from  test  data  were  complex,  displaying  Imaginary  portions  of  considerable 
significance.  However,  the  NASTRAN  model  of  the  AH-1G  that  was  used  In 
the  analysis  considers  only  the  undamped  situation  In  the  extraction  of 
modal  parameters  yielding  the  classical  normal  mode  or  orthonormal  mode 
shapes.  Thus,  the  elements  of  each  mode  are  scalar  quantities  containing 
no  relative  phase  angle  Information  among  the  modal  elements.  Consequent¬ 
ly,  to  qualitatively  and  quantitatively  compare  the  complex  modes  developed 
from  test  data  with  the  real  modes  computed  from  the  NASTRAN  model  of  the 
AH-1G,  the  complex  modes  extracted  from  test  measurements  required  trans¬ 
formation  to  the  real  form.  This  process  of  converting  complex  modes  to 
real  modes  compromises  the  physical  realities  of  the  structure  but  is 
necessary  if  a  comparison  of  analytical  and  test  results  is  to  be  effect¬ 
ed.  The  technique  which  was  implemented  for  conversion  of  the  complex 
modes  to  real  modes  is  a  judgemental  process  based  on  a  reference  phase 
angle.  Let  4^  be  the  jth  element  of  the  kth  complex  mode,  y. 

The  complex  representation  of  a  modal  element  is 


'j = 


a  •  +  i 
J 


where  i  =  /T" 


The  modal  element  can  also  be  represented  by  a  magnitude  and  a  phase 
angle 


¥ 


j 


0  . 
J 


where  e.  =  tan 

J 


-1 


and  0°  <  9  i  360°. 
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Based  on  a  review  of  test  data  the  complex  axes  system  was  rotated  45° ,  as 
shown  on  Figure  8,  to  accommodate  the  transformation  from  complex  modes  to 
real  modes . 


Fiqure  8.  Coordinate  transformation. 

The  real  representation  of  the  complex  modes  was  based  on  following 
phase  angle  criteria 
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TRUNCATED  MODEL 


The  equation  of  motion  of  a  linear  structure  can  be  written  as 

[M]  {y}  ♦  Cd3  {y}  +  [K]  {y}  -  {f}  (13) 

assuming  a  steady  state  solution  of  the  form 
{y}  «  (y}eiu,t  and  {f}  -  {f}e1a)t 
and  substituting  into  equation  (13)  yields 

[  -  [M]u2  +  1W[D]  +  [K]j  (yl  *  {f}  (14) 

Define  the  displacement  Impedance  as 

[Z]  =  [ZR]  +  [ZX]=  [K]  -  u)2  [M]  +  i  uj[D]  (15) 

thus 

[Zj  (y)  =  m 

The  element  displacement  impedance  can  be  expressed  as 

The  displacement  impedance  matrix  defined  in  equation  (15)  can  be  pre¬ 
multiplied  by  [4>]~T  [$]T  and  postmultiplied  by  [$]  [$]'^  where  [<J>]  is 
the  matrix  of  modal  vectors 

[Z]  =  [*rT  [$]T[K]  [4>]  -  J  [$]T  [M]  [$]  +  i co  [$]T  [D]  [*]  [*]-1  (16) 
where  i  =  /T] 
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Using  the  orthogonality  characteristic  of  the  normal  modes,  the  diagonal 
generalized  mass  can  be  expressed  as 


[M*3  =  |>]T  [M]  [»]  (17) 

and  the  diagonal  generalized  stiffness  can  be  expressed  as 

DK*J  =  [t]T  [K]  [*]  (18) 

assuming  structural  damping  which  is  proportional  to  stiffness,  where  g  is 
the  structural  damping  coefficient,  yields 

Of  to]  M  *  EgK*J  (19) 

Substituting  equations(17) ,  (18)  and  (19)  into  equation  (16)  gives 

[Z]  =  0]"T  CK*  -  oj2  M*  +  i  g  K*J  (20) 

Define  the  modal  displacement  impedance  matrix  as 

£Z*>]  =  £K*  -  a)2  M*  +  i  g  K*J  (21) 

and  when  substituted  into  equation  (20)  yields 

[Z]  -  tCT  D-Z*  J  [♦]-’  (22) 

The  displacement  mobility  matrix  at  a  forcing  frequency  to  has  elements 
which  represent  the  ratio  of  the  displacement  phasor  along  coordinate  j 
to  the  applied  force  phasor  along  coordinate  k  when  the  applied  forces 
at  all  other  coordinates  are  zero.  The  full  displacement  mobility 
matrix  is  defined  as 

[V]  =[Z]_1  (23) 
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and  can  be  obtained  from  equation  (22),  thus 


[V]  -  M  of  =  M  fY*J  [»f 


(24) 


with  force  excitation  applied  at  a  particular  coordinate  k  the  response 
at  each  station  of  interest  j  can  be  obtained  yielding  a  column  of  the 
displacement  mobility  given  in  equation  (24),  thus 


As  indicated  in  equation  (25),  the  exact  displacement  mobility  of  a  physi¬ 
cal  structure  is  obtained  only  by  considering  an  infinite  number  of  nor¬ 
mal  modes.  Practical  considerations,  however,  restrict  the  number  of 
normal  modes  that  can  be  used  in  an  analysis. 


A  truncated  model  of  a  linear  system  can  be  formed  by  considering  fewer 
normal  modes  than  degrees  of  freedom  associated  wi'"i  the  physical  sys¬ 
tem.  If  the  number  of  normal  modes  is  limited, the  displacement  mobility 
becomes 

N 

(Y0(k)}  '  w  {Y**(k)i}+  ({,=  ^VVuoi'^i  +  {s)  (26) 

where  {6}  represents  the  difference  between  tne  actual  displacement 
mobility  and  the  approximation  using  only  N  normal  modes. 


A  truncated  model  of  a  linear  system  is  comprised  of  the  constant  coef¬ 
ficient  matrices  in  the  equation  of  motion  formed  using  fewer  normal 
modes  than  degrees  of  freedom  such  that  the  scalar  product  of  (6)  T  (6)  is 
minimized  with  respect  to  the  modal  coefficients  (Y*>.  Using  equa¬ 
tion  (26),  the  vector  representing  the  difference  between  the  exact  and 
the  approximation  to  the  displacement  mobility  can  be  written  as 

{Vj(k)}-  rt(k)iJ  «*}  (27) 
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and  the  transpose  of  {6}  becomes 


{ 6 } T  » 


(vj(k)}T-<v*>T^W  [*]T 


final ly 


(28) 


t{)T{{)  ■  {vj(k)}T  b(k)}  -  h(k)}T  W  c  V),*  <**> 

-  (Y*)T  E-  *(k)1J  WT  (vJ(k)} 

♦  ty*>T  M(k)1  -]  [$]^[4i]f‘ J  (v*) 

(29) 


Equation  (29)  is  a  scalar  product  and  the  derivative  of  a  scalar  with 
respect  to  a  vector  is  also  a  vector.  Thus,  the  modal  vector  {Y*}  must 
be  obtained  which  will  minimize  the  scalar  product  presented  in 
equation  (29).  This  can  be  accomplished  by  differentiating  the  scalar 
product,  equation  (29),  with  respect  to  the  modal  displacement  mobility 
vector  and  equating  the  result  to  zero  as  follows: 


{<S}T{6} 


3{Y*} 


0  ■  *(k)i J[t]T  {Yj(k)| +  ^(kM  -atrfc™  *(k)1 

(30) 


Solving  for  the  modal  displacement  vector  for  which  the  scalar  product  is 
a  minimum  yields 


(31) 


The  matrix  [i>]+  is  commonly  referred  to  as  the  generalized  inverse  or 
pseudoinverse  of  [<f>]  and  is  defined  by 


Of  if(»]T[t] 


-1 


[*] 


T 
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(32) 


Equation  (31)  may  be  rewritten  as 

f»(k)1JtY*)^v*g|»(k))J  =  W+|yj(k)}  (33) 

For  the  full  displacement  mobility  matrix  equation  (26)  becomes 

[Y]  =  [4]  fY*  J  [tf  +  [S]  (34) 

and  equation  (33)  can  be  extended  to 

IS  Y*  J  =  [4]+[Y]  [4]+T  (35) 

Equation  (35)  can  be  premultiplied  by  [4]  and  postmultiplied  by  [4]T  to 
define  the  truncated  mobility  matrix  as 

[Y] tr  s  [4]  C-Y*  i  [4]T  =  [4]  [4]+  [Y]  [4]+T  [4]T  (36) 

The  truncated  impedance  matrix  is  given  as 

[Z] TR  =  [4]+T  f  [4]+  =  W+T  IS  Z*  ;  [4]+  (37) 


The  truncated  matrices  are  of  order  J,  the  number  of  coordinates  of 
interest,  and  rank  N,  the  number  of  modes  considered,  and  are  therefore 
singular. 

If  the  product  of  [Y]TR  [Z]yR  is  premultiplied  by  [$]+  and  postmulti¬ 
plied  by  [*]  the  result  is 

M+  [Y]tr  [Ztr]  [*]  =  I  (38) 

and 

[4]T  [Z]tr  [Y]tr  [4]+T  =  I 
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or  using  equation  (43) 

T  «  7  (q}T  [4>]T  [H]  [*]  {q}  *  \  (q}T  C*  M*  J{q}  (45) 

Using  equation  (39)  the  kinetic  energy  becomes 

T  =  \  (q}T  [$]T  [H]tr  [$]  {q}  (46) 

Consider  a  structural  mass  and  stiffness  change  wherein  the  resulting 
displacements,  in  the  frequency  range  of  interest,  can  be  expressed 
approximately  by 

(y)  =  [$]  (r)  (47) 

The  kinetic  energy  for  the  modified  structure  can  be  expressed  approxi¬ 
mately  as 

T  =  \  (f)T  [<I>]T  ([M]  +  [AM]]  [<j>]  {f>  (48) 

which  upon  expansion  becomes 

T  =  \  (r}T  [<t]T  [M][$]  (f)  +  \  { r )T  [<1>]T[aM][<I>]  {f}  (49) 

Considering  only  the  portion  of  the  equation  which  contains  the  modifi¬ 
cation  in  the  mass 

[4>]T[AM][$]  =  [$]T[<i>]  T  [<I>]T[$]  "T  [$]T[AM][$]  [<D]T[4»]  -1  [<D]T[4>] 

.  J  4 

[$]T[aM][<J>]  =  [*]T[*]+T[$]T  [aM][*][*]+  [<J>]  (50) 

Defining  AMjR  =[$]+^[4>]T[AM][$][4>]+  the  kinetic  energy  becomes 

T  =  \  Cr}J  [$]Tf[M]TR  +  [AM]tr|  [*]  Cr]  (51) 

V  4 
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The  potential  energy  of  the  structure  Is  given  approximately  by 


where 


V  =  J  <r>T[*] 
=  \  (r)T[i] 

tAK]TR  5 


T([K]  +  [AK]J  MM 

T  [K]tr  +  [AK]TR]MM 


M+TMT[«]MM+ 


(52) 


The  truncated  influence  coefficient  matrix  for  the  modified  structure 
is  given  by 


^MOD^TR 


+  ^TR 


(53) 


and  the  eigenvalue  equation  which  reflects  the  mass  and  stiffness 
changes  is 


£CM0dWMTR  +  AMTR^$M0D^  "  ^MOD^ 


7 


MOD 


(54) 


Equations  Using  Orthonormal  Mode  Vectors 

It  is  more  convenient  to  use  the  orthonormal  mode  in  the  previous  equa¬ 
tions,  thereby  eliminating  the  effect  of  normalization  on  each  mode. 

The  orthonormal  mode  allows  a  comparison  of  the  relative  prominence  of 
each  mode.  The  orthonormal  mode  is  defined  as 

(55) 


(56) 


{4'}j  [M]m.  =  i 

For  N  modes  this  may  be  extended  to 
[T]T[H]  M  ■  [i] 
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The  generalized  mass  associated  with  each  mode  can  be  written  as 


M*  =  W]  [M]{4>>  1 


(57) 


If  equation  (55)  is  multiplied  by  the  generalized  mass  of  the  ith  mode 
and  compared  to  equation  (57),  the  elements  of  the  orthonormal  mode 
shape  become 


or 


(58) 


• 

Ji 


For  N  modes  and  J  points  of  interest  on  the  structure  equation  (58) 
may  be  expanded  to 

[*]  =  M  [V^  J  (59) 


Some  fundamental  matrix  operations  which  will  be  used  subsequently  are 

[t]T  ,  jo]t 


or  - 

M+  = 


mpj  c»]T  m  pmpjm 


where 


W 


[»]T[f]  -’of  =  — 

/tp 


M 


h]+  -  mtm 


r-'+T  -  W+T 


(60) 


If  the  above  relationships  are  substituted  into  equation  (36)  the 
truncated  mobility  matrix  becomes 
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(61) 


MTR  -  [¥][¥]+[Y][*]+T[¥]T 

=  [4\)D  Y*M*  ][f]T 

and  equation  (37)  becomes 

[Z]TR  =  tt]+TC  (62) 

The  truncated  mass,  stiffness  and  damping  matrices  presented  in 
equation  (39)  can  be  written  as  a  function  of  the  orthonormal  modal 
matrices  and  become 

(xrR  =■  [»]+Tcnw+ 

[K]tr  *  W+TC  j£  JW+  (63) 

•  W+Tr  a2  JM+ 

[3]tr  =  M+Tt  ga2  JM+ 

The  influence  coefficient  matrix,  using  the  orthonormal  mode  matrix 
is 

[c]TR  =  MC  $  dWT 

(64) 

The  classical  modal  eigenvalue  equation  for  the  truncated  model  can- 
be  written  as 

[c]TRrM]TRm  =  \  m  (65) 

ai 


47 


Considering  the  mass  and  stiffness  changes  to  the  structure 


AKtr  =  m+T[»]T[4K]t1lK*]<' 


The  modified  truncated  influence  coefficient  matrix,  in  terms  of  the 
orthonormal  mode  matrix, is  given  by 


tCM0D^TR  "  ^ 


CK]Tr  + 


[ak]tr  [V] 


-!mT 


and  the  eigenvalue  equation  which  reflects  the  mass  and  stiffness  changes 
is 


^MnrJ  £M]Td  +  [AM] 


TR 


Jtr 


]  [* 


Mnn- 


(68) 
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TRUNCATED  MODEL 


COMPUTER  EXPERIMENTS 


Computer  experiments  were  performed  on  mass  and  stiffness  changes  to 
determine  the  number  of  degrees  of  freedom  required  to  Insure  adequate 
representation  of  the  Truncated  Model  method  described  in  Reference  4. 

The  computer  experiments  were  conducted  for  mass  and  stiffness  changes 
individually  and  in  combination.  The  computer  experiments  were  per¬ 
formed  using  computer  codes  developed  for  the  IBM  360/50  D.O.S./V.S. 
Computer  System  using  the  FORTRAN  IV  Language. 

Truncated  parameters  of  a  linear  system  are  defined  herein  as  those  con¬ 
stant  coefficient  matrices  in  the  equations  of  motion  formed  from  fewer 
orthonormal  modes  than  degrees  of  freedom  such  that  the  Euclidean  Norm 
of  the  difference  between  the  exact  mobility  matrix  and  the  approxima¬ 
tion  to  the  mobility  matrix,  which  is  formed  using  the  orthonormal  modes, 
generalized  masses,  and  generalized  stiffnesses,  is  minimum  with  respect 
to  the  modal  or  generalized  mobility.  The  mobility  matrix  consists  of 
elements  which  are  defined  generically  as  the  partial  derivative  of  the 
response  at  coordinate  i  with  respect  to  force  excitation  applied  at 
coordinate  j  while  the  applied  force  at  all  other  coordinates  is  zero. 

Development  of  the  truncated  model  requires  the  structure's  natural 
frequencies,  orthonormal  modes  and  modal  damping  which  charac¬ 
terize  the  frequency  spectrum  of  interest.  The  natural  frequencies  and 
orthonormal  mode  shapes  necessary  to  effect  the  truncated  model  of  the 
AH- I G  airframe  were  generated  from  a  92  degree  of  freedom  NASTRAN  model. 
The  truncated  model  was  developed  using  92  degrees  of  freedom  and  16 
elastic  modes  which  encompassed  the  frequency  range  of  from  0  to  40  Hz. 


4.  Flannelly,  W.  G.,  and  Giansante,  N.,  EXPERIMENTAL  VERIFICATION  OF 
SYSTEM  IDENTIFICATION,  Kaman  Aerospace  Corporation,  Report  Number 
R- 1 275 ,  USAAMRDL  TR-74-64,  Eustis  Directorate,  U.  S.  Army  Air  Mobility 
Research  and  Development  Laboratory,  Fort  Eustis,  Virginia  23604, 
August  1974,  AD  784190. 


49 


The  procedure  used  In  developing  the  truncated  models  used  in  the  computer 
study  is  shown  schematically  in  Figure  9.  The  92  degree  of  freedom  NASTRAN 
model  presented  in  Reference  1  was  modified  to  be  compatible  with  the  mass 
configuration  of  the  baseline  AH-1G  helicopter  airframe  (S/N  57-15683), 
which  was  subjected  to  ground  vibration  tests.  Additionally,  the  number 
of  degrees  of  freedom  was  reduced  to  92  to  correspond  with  the  locations 
of  the  accelerometers  positioned  on  the  airframe.  A  92  degree  of  freedom 
truncated  model  was  subsequently  formed  using  the  normal  mode  shapes, 
natural  frequencies,  and  generalized  masses  extracted  from  the  reduced 
NASTRAN  model . 

In  generating  the  truncated  model,  the  normal  mode  shapes  obtained  from 
the  NASTRAN  model  were  converted  to  orthonormal  mode  shapes.  The  ortho¬ 
normal  mode  shape,  which  is  the  conventional  normal  mode  shape  normalized 
by  the  square  root  of  the  generalized  modal  mass,  has  physical  significance 
and  offers  a  convenient  means  of  comparing  the  relative  importance  of  the 
modal  responses  of  a  structure.  In  addition,  the  product  of  the  transpose 
of  the  orthonormal  mode  vector  and  the  orthonormal  mode  vector,  which 
yields  a  scalar,  provides  a  direct  measure  of  the  significance  of  a 
particular  mode,  since  the  numerically  greater  the  product  the  more  im¬ 
portant  the  contribution  of  that  modv*  to  the  overall  response  of  the  struc¬ 
ture.  On  the  contrary,  the  conventional  normal  mode  shape,  which  depends 
on  the  normalizing  station,  has  no  physical  significance  and  taken  inde¬ 
pendently  does  not  yield  any  information  regarding  the  dynamic  response 
characteristics  of  the  structure. 

As  indicated  in  Figure  9,  step  3,  mass  and/or  stiffness  changes  were  in¬ 
corporated  in  the  92-degree-of-freedom  truncated  model.  Each  of  the 
resulting  models,  which  included  particular  mass  changes  and  stiffness 
changes,  either  individually  or  in  combination,  became  a  baseline  model 
that  was  used  as  a  basis  of  comparison  for  the  reduced  truncated  models 
developed  subsequently  which  incorporated  the  same  mass  and/or  stiffness 
change  configuration.  The  reduced  truncated  models  contained  50,  30, 
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Figure  9.  Schematic  truncated  model  development. 
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and  20  degrees  of  freedom.  The  mass  changes,  stiffness  changes,  and  combi¬ 
nations  of  mass  and  stiffness  that  were  Incorporated  In  the  models  are 
presented  In  Tables  2,  3,  and  4,  respectively.  These  mass  and  stiffness 
change  configurations  required  the  formation  of  18  baseline  models  which 
were  used  to  evaluate  the  number  of  degrees  of  freedom  required  In  the 
reduced  models  to  substantiate  the  Truncated  Model  methodology. 

In  developing  the  reduced  truncated  models,  the  mass  and  stiffness  changes 
were  applied  at  grid  points  corresponding  to  the  grid  points  at  which 
the  same  changes  were  Incorporated  In  the  92  degree  of  freedom  models. 
Reduced  truncated  models  utilizing  50,  30,  and  20  degrees  of  freedom 
were  developed  for  the  18  baseline  truncated  models,  and  the  number  of 
modes  retained  varied  from  5  to  16  modes. 

Review  of  preliminary  vibration  shake  test  data  of  the  AH-1G  helicopter 
airframe  indicated  an  absence  of  dynamic  response  in  the  longitudinal 
(fore  and  aft)  direction  for  the  frequency  spectrum  of  interest.  Refer¬ 
ence  5  also  showed  a  lack  of  longitudinal  modes  except  for  the  pylon  fore 
and  aft  mode  and  the  main  rotor  mast  longitudinal  bending  mode.  There¬ 
fore,  the  modelling  strategy  in  developing  the  reduced  truncated  models 
was  to  omit  the  longitudinal  degrees  of  freedom,  except  for  those  associated 
with  the  main  rotor  mast.  Generally,  the  grid  points  at  which  concentrated 
masses  were  located  were  also  included  in  the  reduced  truncated  models. 
Additionally,  coordinates  at  which  mass  and/or  stiffness  changes  were 
effected  were  retained  in  the  reduced  models.  One  further  consideration 
in  generating  the  truncated  models  was  to  include  coordinates  which  were 
reasonably  distributed  throughout  the  structure  rather  than  to  include 
only  the  coordinates  concentrated  in  a  region  of  high  dynamic  response. 


Cronkhite,  James  D.,  and  Berry,  Victor  L.,  CORRELATION  OF  THE  AH-1G 
AIRFRAME  TEST  DATA  WITH  A  NASTRAN  MATHEMATICAL  MODEL.  NASA  CR-145119, 
February  1976. 
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TABLE  2.  MASS  CHANGES  FOR 

COMPUTER  EXPERIMENTS 


Outboard  wing  stores  removed 

Inboard  wing  stores  removed 

Inboard  and  outboard  wing  stores  removed 

p 

Mass  added  at  c.g.  location  (1.296  lb  sec  /in.) 

2 

Blade  and  hub  mass  added  at  hub  location  (1.86  lb  sec  /in.) 

Inboard  and  outboard  wing  stores  removed,  mass  added  at  c.g.  location 
blade  and  hub  mass  added  at  hub  location. 


TABLE  3.  STIFFNESS  CHANGES  FOR 
COMPUTER  EXPERIMENTS 

Add  lateral  stiffness  between  grid  points  4633T2,  4637T2 
(4.56  x  10*  lb/in.) 

Add  lateral  stiffness  between  grid  points  42245T2,  46445T2 
(1  x  106  Ib/in. ) 

Reduce  lift  link  stiffness  by  0.5  x  106  lb/ in. 

Add  stiffness  between  grid  points  4633T2,  4637T2  (4.56  x  105  lb/in.) 

Add  stiffness  between  grid  points  40145,  52045  (8.17  x  104  Ib/in.) 

Add  stiffness  between  grid  points  40145,  52045  (1.61  x  105  Ib/in.) 

Add  stiffness  between  grid  points  40145,  52045  (2.16  x  105  lb/in.) 

Add  stiffness  between  grid  points  42245T2,  46445T2  (.5  x  106  lb/ in.) 


TABLE  4.  MASS  AND  STIFFNESS  CHANGES 


Clean  wing  configuration 

Add  stiffness  between  grid  points  40145,  52045  (8.17  x  104  Ib/in.) 

Inboard  and  outboard  wing  stores  removed,  mass  added  at  c.g.  location, 
blade  and  hub  mass  added  at  hub  location 

Add  stiffness  between  grid  points  40145,  52045  (8.17  x  104  lb/in.) 

Inboard  wing  stores  removed 

Add  stiffness  between  grid  points  42245T2,  4644 5T2  (0.5  x  106  lb/ in.) 
Inboard  wing  stores  removed 

Add  stiffness  between  grid  points  40145,  52045  (8.17  x  104  lb/in.) 


This  procedure  yields  a  rectangular  modal  matrix  which  can  be  reliably 
operated  on  using  the  pseudoinverse  technique  and  which,  in  general,  yields 
numerically  stable  and  physically  meaningful  results.  Distribution  of 
coordinates  for  typical  reduced  truncated  models  having  50,  30,  and  20 
degrees  of  freedom  for  the  AH-1G  helicopter  airframe  in  the  clean  wing 
configuration  is  presented  in  Table  5.  In  performing  the  computer  experi¬ 
ments  for  a  particular  mass  and/or  stiffness  change  configuration,  no 
attempt  was  made  to  optimize  the  reduced  truncated  models  in  relation  to 
the  92  degree  of  freedom  model  incorporating  the  same  mass  and  stiffness 
changes.  However,  engineering  judgement  was  applied  in  selecting  the 
reduced  degrees  of  freedom  used  in  any  model. 

One  of  the  prime  objectives  of  the  Truncated  Model  methodology  is  to 
develop  reduced  models  which  can  successfully  predict  the  effects  of 
changes  in  mass  and  stiffness  of  the  structure.  The  reduced  truncated 
model  allows  the  rapid  determination  of  variations  in  configuration  and 
provides  a  means  for  evaluating  a  variety  of  structural  changes  in  the 
interest  of  improving  the  dynamic  response  of  the  vehicle.  In  the  computer 
experiments  which  were  conducted,  the  success  of  a  reduced  truncated  model 
in  adequately  predicting  mass  and  stiffness  changes  was  measured  by  com¬ 
paring  modal  parameters  obtained  from  the  reduced  truncated  model  with 
those  obtained  from  the  92  degree  of  freedom  baseline  for  the  same 
mass  and  stiffness  changes.  The  modal  data  which  were  compared  included 
natural  frequencies,  generalized  masses,  and  the  orthonormal  mode  shapes 
in  the  frequency  range  of  interest.  Using  the  orthonormal  mode  shape 
criterion,  a  single  numerical  value  was  determined,  for  each  mode  of 
interest,  which  provided  a  basis  of  comparison  of  the  effectiveness 
of  the  reduced  truncated  model.  For  a  specific  mass  and  stiffness 
change  configuration  this  numerical  value,  which  will  be  referred  to  as 
the  Modal  Correlation  Number,  was  calculated  for  each  mode  as  the  summa¬ 
tion  of  the  differences  between  the  orthonormal  modal  elements  of  the 
baseline  92  degree  of  freedom  truncated  model  and  the  respective  modal 
elements  of  the  reduced  truncated  model  divided  by  the  number  of  degrees 
of  freedom  in  the  reduced  truncated  model.  The  Model  Correlation  Number, 
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TABLE  5.  REDUCED  TRUNCATED  MODELS 

INBOARD,  OUTBOARD  WING  STORES  REMOVED 


92  D.O.F. 
Index 

Degree  of  Freedom 
Identification 

50  D.O.F. 

30  D.O.F. 

20  D.O.F. 

1 

40141T3 

★ 

X 

X 

X 

2 

40142T3 

X 

X 

4 

38045T2 

X 

X 

X 

5 

38045T3 

X 

X 

X 

6 

33845T2 

X 

7 

33845T3 

X 

9 

40145T3 

X 

10 

42245T2 

X 

X 

11 

42245T3 

X 

14 

46445T2 

X 

X 

X 

15 

46445T3 

X 

X 

X 

17 

48045T2 

X 

X 

18 

48045T3 

X 

X 

X 

19 

48845T2 

X 

X 

20 

49745T3 

X 

23 

52045T2 

X 

X 

X 

24 

52045T3 

X 

X 

X 

26 

2001 29T1 

X 

X 

X 

27 

2001 29T2 

X 

X 

X 

29 

2001 53T1 

X 

X 

X 

30 

2001 53T2 

X 

X 

X 

33 

1 24800T3 

X 

35 

25065T3 

X 

36 

19752T3 

X 

37 

19765T3 

X 

39 

25045T2 

X 

40 

26861 T3 

X 

41 

26865T2 

X 

42 

26869T3 

X 

*  "x"  Indicates  degree  of  freedom  used  in  the  model 
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TABLE  5.  CONTINUED 


92  D.O.F. 
Index 

Degree  of  Freedom 
Identification 

50  D 

45 

30045T3 

X 

47 

22703T2 

X 

48 

22725T3 

X 

50 

21325T2 

X 

52 

21341T3 

X 

54 

64229T3 

X 

57 

65929T3 

X 

58 

74229T3 

X 

61 

75929T3 

X 

63 

18635T2 

X 

64 

18625T3 

X 

66 

15625T2 

X 

67 

15629T3 

69 

15621T3 

X 

72 

13843T2 

X 

74 

14809T2 

X 

76 

11563T3 

X 

79 

9371T3 

X 

80 

9363T2 

X 

|  81 

9331T3 

83 

9333T2 

X 

84 

9339T3 

X 

86 

6167T2 

X 

87 

6169T3 

X 

89 

4661T3 

X 

90 

4633T2 

X 

91 

4637T2 

92 

3331T3 

X 

which  is  determined  for  each  mode  of  interest,  can  be  expressed  algebra! 
cally  as: 


MODAL  CORRELATION  NO.  * 


1 

r 


;/ 

j«A  j192  D.O.F. 


1/2 


In  the  above  expression, represents  the  orthonormal  modal  ele- 

J192  D.O.F. 

ment  for  the  jth  coordinate  and  the  ith  mode  for  the  92  degree  of  free¬ 
dom  truncated  model  for  a  particular  mass/ stiffness  change  configuration. 

The  term  4*^  represents  the  corresponding  orthonormal  modal  element 
J1RED 

for  the  reduced  degree  of  freedom  model.  Thus,  the  Modal  Correlation  Num 
ber  can  be  considered  as  representative  of  the  average  of  the  differences 
between  the  orthonormal  mode  elements  obtained  from  the  92  degree  of  free 
dom  baseline  model  and  the  corresponding  modal  elements  calculated  from 
a  reduced  degree  of  freedom  truncated  model. 


Approximately  90  separate  computer  experiments  were  performed  in  the  study. 
Table  6  presents  results  for  typical  reduced  degree  of  freedom  truncated 
models  for  the  mass  change  configuration  of  inboard  wing  stores  removed. 

A  comparison  of  natural  frequency  and  generalized  mass  is  shown  in  Table 
6a  and  6b,  respectively,  for  the  92  degree  of  freedom  truncated  model 
and  the  reduced  truncated  models  with  50,  30  and  20  degrees  of  freedom. 

Data  for  the  unmodified  structure  with  inboard  and  outboard  wing  stores 
attached  are  also  presented  in  the  table.  Removal  of  the  inboard  wing 
stores  does  not  significantly  change  the  fuselage  natural  frequencies. 

A  comparison  of  the  predicted  natural  frequencies  shows  excellent  agree¬ 
ment  among  the  models.  A  review  of  the  generalized  mass  predictions  indi¬ 
cates  good  agreement  with  the  92  degree  of  freedom  baseline  model  except 
for  the  second  mode  in  the  50  degree  of  freedom  model,  the  first,  second, 
and  twelfth  mode  for  the  30  and  20  degree  of  freedom  modes,  and  the 
thirteenth  mode  for  the  20  degree  of  freedom  model.  The  lack  of  correl¬ 
ation  for  the  aforementioned  modes  occurred  because  the  coordinates 
displaying  significant  response  for  the  modes  of  interest  were  omitted 
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TABLE  6.  MASS  CHANGE  CONFIGURATION  -  INBOARD  WING  STORES  REMOVED 
30  AND  20  DOF  MODELS  CONTAIN  NO  DOF  ON  ROTOR  MAST 


a.  Natural  frequency  comparison 


w/Stores 

BHMWBWSWfil  t  i  I’ll  if  11 \  ltf  i  Tr  mi  1 1 1  — 

Mode 

92  D.O.F. 
Unmodified 

92  D.O.F. 
Basel ine 

50  D.O.F. 

30  D.O.F. 

20  D.O.F. 

1 

4.18 

4.33 

4.33 

4.33 

4.33 

2 

5.12 

5.12 

5.12 

5.12 

5.12 

3 

6.87 

6.90 

6.90 

6.90 

6.90 

4 

8.02 

8.11 

8.11 

8.11 

8.11 

5 

12.39 

12.39 

12.39 

12.39 

12.39 

6 

15.75 

16.35 

16.35 

16.35 

16.35 

7 

16.37 

16.84 

16.84 

16.85 

16.84 

8 

18.85 

18.97 

18.96 

18.96 

18.97 

9 

20.60 

21.09 

21.09 

21.09 

21.10 

10 

23.93 

23.96 

23.96 

23.97 

23.97 

11 

24.70 

24.74 

24.74 

24.72 

24.72 

12 

27.62 

27.65 

27.65 

27.66 

27.65 

13 

29.95 

29.94 

29.95 

29.95 

29.96 

14 

31.06 

31.07 

31.06 

31.07 

31.02 

15 

33.13 

33.23 

33.23 

33.24 

33.23 

b.  Generalized  mass  comparison 


w/Stores 

Inboard  Winq 

Stores  Removed 

Mode 

92  D.O.F. 
Unmodified 

92  D.O.F. 
Baseline 

50  D.O.F. 

30  D.O.F. 

20  D.O.F. 

1 

1.089 

1.017 

1.017 

8.587 

8.583 

2 

0.842 

0.842 

1.812 

6.782 

6.825 

3 

0.404 

0.408 

0.408 

0.408 

0.408 

4 

0.732 

0.709 

0.709 

0.709 

0.709 

5 

7.524 

7.537 

7.534 

7.533 

7.577 

6 

3.242 

0.885 

0.886 

0.886 

0.884 

7 

0.832 

2.582 

2.583 

2.588 

2.595 

8 

2.831 

3.510 

3.507 

3.470 

4.829 

9 

2.754 

2.767 

2.767 

2.765 

2.764 

10 

4.575 

4.709 

4.714 

4.753 

4.708 

11 

2.751 

2.777 

2.779 

2.766 

2.793 

12 

1.506 

1.484 

1.484 

5.708 

5.724 

13 

2.369 

2.404 

2.378 

2.355 

3.601 

14 

1.796 

1.835 

3.44 

3.525 

3.236 

15 

1.172 

1.159 

1.159 

1.159 

1.158 
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from  the  model.  Modes  one  and  two  are  the  main  rotor  pylon  pitch  and 
pylon  roll  modes,  respectively;  therefore,  coordinates  on  the  rotor  mast  in 
the  longitudinal  and  lateral  directions  should  be  included  in  the  models 
to  satisfactorily  predict  the  associated  modal  parameters.  Modal  data  present¬ 
ed  in  Table  7  are  also  for  the  inboard  wing  stores  removed  configuration. 

The  reduced  truncated  models  for  which  the  data  were  developed  included 
coordinates  on  the  rotor  mast  resulting  in  improved  correlation  of  the 
modal  parameters.  However,  as  shown  in  Table  7b  the  generalized  masses 
predicted  for  the  eighth  and  thirteenth  mode  for  the  20  degree  of  freedom 
truncated  model  lack  agreement  with  the  generalized  masses  predicted  by 
other  models  for  the  corresponding  modes.  Further  refinement  of  the  coor¬ 
dinates  selected  for  the  model  would,  undoubtedly,  provide  improved  corre¬ 
lation  in  the  generalized  masses  for  the  aforementioned  modes.  The  pre¬ 
dicted  natural  frequencies  are  not  particularly  sensitive  to  the  specific 
coordinates  embodied  in  the  specific  models  as  long  as  the  coordinates 
are  reasonably  distributed  throughout  the  structure.  Table  8  presents  a 
review  of  the  orthonormal  mode  shapes,  in  terms  of  the  Modal  Correlation 
Number,  for  the  50,  30,  and  20  degree  of  freedom  models.  Generally,  the 
degree  of  correlation  for  each  mode  decreases,  as  witnessed  by  an  increase 
in  the  Modal  Correlation  Numbers,  as  the  number  of  degrees  of  freedom  in 
the  model  is  reduced.  Additionally,  the  ability  to  predict  the  higher  fre¬ 
quency  orthonormal  mode  shapes  decreases  as  the  model  size  decreases. 

Modal  data  for  the  configuration  with  inboard  and  outboard  wing  stores 
removed  are  presented  in  Table  9.  The  data  presented  in  Table  9a  indicate 
excellent  correlation  for  predicted  natural  frequency  regardless  of  size 
or  specific  degrees  of  freedom  contained  in  the  model.  Generalized  mass 
data  as  shown  in  Table  9b  are  sensitive  to  the  model  size  and,  more  parti¬ 
cularly,  dependent  upon  the  coordinates  included  in  the  model.  Specifi¬ 
cally,  the  50  degree  of  freedom  model  does  not  include  coordinates  required 
to  define  the  second  mode;  therefore,  there  is  no  correlation  of  generalized 
mass  compared  to  the  92  degree  of  freedom  baseline  model.  Similarly,  the 
30  and  20  degree  of  freedom  models  do  not  predict  satisfactorily  the  first 
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TABLE  7.  MASS  CHANGE  CONFIGURATION  -  INBOARD  WING  STORES  REMOVED 
MODELS  CONTAIN  DOF  ON  ROTOR  MAST 


a.  Natural  frequency  comparison 


w/Stores 

Inboard  Winq  Stores  Removea  j 

Mode 

92  D.O.F. 
Unmodified 

92  D.O.F. 
Baseline 

50  D.O.F. 

30  D.O.F. 

20  D.O.F. 

4.18 

4.33 

4.33 

4.33 

4.33 

5.12 

5.12 

5.12 

5.12 

5.12 

6.87 

6.90 

6.90 

6.90 

6.90 

8.02 

8.11 

8.11 

8.11 

8.11 

12.39 

12.39 

12.39 

12.39 

12.40 

15.75 

16.35 

16.35 

16.35 

16.35 

16.37 

16.84 

16.84 

16.84 

16.84 

18.85 

18.97 

18.97 

18.97 

18.96 

20.60 

21.09 

21.09 

21.09 

21.09 

10 

23.93 

23.96 

23.96 

23.97 

23.97 

11 

24.70 

24.74 

24.74 

24.74 

24.71 

12 

27.62 

27.65 

27.65 

27.65 

27.65 

13 

29.95 

29.94 

29.94 

29.95 

29.95 

14 

31.06 

31.07 

31.07 

31.07 

31.07 

15 

33.13 

33.23 

33.23 

33.23 

33.23 

b.  Generalized  mass  comparison 


“w/sTBKB - 

r  Inboard  Winq  Stores  Removed  i 

92  D.O.F.'  “ 

92  D.O.F. 

50  D.O.F. 

30  D.O.F. 

20  D.O.F. 

Mode 

Unmodified 

Basel ine 

1.089 

1.017 

1.017 

1.017 

1.017 

0.842 

0.842 

0.842 

0.842 

0.842 

0.404 

0.408 

0.408 

0.408 

0.408 

0.732 

0.709 

0.709 

0.709 

0.709 

7.524 

7.537 

7.528 

7.534 

7.539 

3.242 

0.885 

0.886 

0.885 

0.884 

7 

0.832 

2.582 

2.582 

2.588 

2.591 

8 

2.831 

3.510 

3.509 

3.469 

4.846 

9 

2.754 

2.767 

2.767 

2.766 

2.764 

KM 

4.575 

4.709 

4.714 

4.724 

4.749 

■  I  ■ 

2.751 

2.777 

2.778 

2.771 

2.786 

1.506 

1.484 

1.484 

1.484 

1.484 

2.369 

2.404 

2.375 

2.371 

3.393 

1.796 

1.835 

1.800 

1.794 

1.804 

1.172 

1.159 

1.159 

1.159 

1.159 
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TABLE  9. 


MASS  CHANGE  CONFIGURATION 

INBOARD,  OUTBOARD  WING  STORES  REMOVED 

30  AND  20  DOF  MODELS  CONTAIN  NO  DOF  ON  ROTOR  MAST 

a.  Natural  frequency  comparison 


w/stores 

IMIIflMi7ITO.IW.IWhiin^iTiTTgB!H 

jEElMiLfrl 

Removed 

Mode 

D.O.F'.1  ' 
Unmodified 

92  D.O.F. 
Baseline 

30  D.O.F. 

20  D.O.F.  1 

1 

4.18 

4.57 

4.57 

4.57 

4.57 

2 

5.12 

5.12 

5.12 

5.12 

5.12 

3 

6.87 

6.95 

6.95 

6.95 

6.95 

4 

8.02 

8.19 

8.19 

8.20 

8.19 

5 

12.39 

12.41 

12.41 

12.41 

12.41 

6 

15.75 

16.57 

16.57 

16.57 

16.57 

7 

16.37 

17.99 

17.99 

17.99 

17.98 

8 

18.85 

22.26 

22.25 

22.26 

22.28 

9 

20.60 

23.17 

23.17 

23.17 

23.16 

10 

23.93 

24.63 

24.63 

24.63 

24.62 

11 

24.70 

28.00 

27.86 

27.89 

27.87 

12 

27.62 

28.20 

28.18 

28.19 

28. 18 

13 

29.95 

30.41 

30.42 

30.51 

30.41 

14 

31.06 

31.21 

31.20 

31.23 

31.24 

15 

33.13 

33.30 

33.30 

33.31 

33.30 

b. 

Generalized 

mass  comparison 

Outboard  Wing  Stores 

Removed 

Mode 

HEESSflSfllBI 

92  D.O.F. 
Baseline 

50  D.O.F. 

30  D.O.F. 

'  20  D.O.F. 

1 

i.08: 

.922 

.922 

5.997 

5.994 

2 

0.842 

0.842 

6.857 

6.778 

6.880 

3 

0.404 

0.417 

0.417 

0.417 

0.417 

4 

0.732 

0.693 

0.693 

0.693 

0.692 

5 

7.524 

7.472 

7.470 

7.472 

7.520 

6 

3.242 

0.818 

0.818 

0.818 

0.816 

7 

0.832 

1.644 

1.645 

1.648 

1.645 

8 

2.831 

8.013 

7.986 

7.915 

7.863 

9 

2.754 

2.273 

2.271 

2.271 

2.263 

10 

4.575 

5.320 

5.316 

5.342 

5.340 

11 

2.751 

2.348 

2.779 

2.739 

2.775 

12 

1.506 

2.781 

1.264 

9.733 

9.891 

iKm 

2.369 

3.008 

4.887 

4.16 

10.80 

14 

1.796 

2.502 

2.490 

2.79 

2.47 

15 

1.172 

1.154 

1.153 

1.155 

1.150 
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and  second  modes  since  the  models  lack  modal  participation  of  these  modes 
as  a  result  of  the  coordinates  selected  for  the  models.  Table  10  presents 
data  for  models  wherein  the  coordinates  used  Include  the  effect  of  the 
modes  in  question.  Consequently,  the  generalized  masses  predicted  by 
the  reduced  models  agree  with  those  obtained  from  the  92  degree  of  freedom 
baseline  model. 

The  ability  of  a  model,  with  a  particular  coordinate  configuration,  to 
accurately  predict  the  generalized  mass  of  a  specific  mode  must  be  evalu¬ 
ated  with  regard  to  the  relative  significance  of  that  mode.  The  Importance 
of  any  mode  is  directly  proportional  to  the  product  of  the  transpose  of  the 
orthonormal  mode  and  the  orthonormal  mode.  Thus,  the  predictive  capability 
of  any  model  must  be  characterized  by  the  ability  to  successfully  deter¬ 
mine  the  modal  parameters  associated  with  the  significant  modes  in  the  fre¬ 
quency  spectrum  of  interest. 

Table  11  presents  data  for  the  first  fuselage  lateral  bending  orthonormal 
mode  for  the  50,  30,  and  20  degree  of  freedom  models  and  the  corresponding 
modal  elemental  data  which  has  been  extracted  from  the  92  degree  of  freedom 
baseline  model.  Comparison  of  the  modal  data  for  reduced  degree  of  freedom 
models  with  those  from  the  92  degree  of  freedom  baseline  model  indicates 
excellent  correlation  for  all  models.  Table  12  presents  data  similar  to 
that  given  in  Table  11,  except  the  data  is  for  the  second  fuselage  vertical 
bending.  Again,  there  is  excellent  agreement  among  the  orthonormal  elements 
of  the  reduced  models  compared  to  the  92  degree  of  freedom  model.  Thus, 
the  20  degree  of  freedom  truncated  model  is  a  satisfactory  representation 
of  the  actual  structure  with  excellent  modal  predictive  capability. 

Modal  data  are  presented  in  Table  13  for  a  stiffness  change  configuration 
realized  by  introducing  a  stiffening  element  having  a  spring  rate  of 
8.17  x  104  lb/in.  between  grid  points  40145  and  52045.  Grid  point  40145 

is  located  on  the  tail  boom  at  FS  401.33  in.,  BL  55.91  in.,  WL  0.0  in.  and 

grid  point  52045  is  located  on  the  vertical  tail  at  FS  520.67  in.,  BL  118.27 

in.,  arid  WL  0.0  in.  Table  13a  and  13b  present  natural  frequency  and 
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TABLE  10.  MASS  CHANGE  CONFIGURATION 

INBOARD,  OUTBOARD  WING  STORES  REMOVED 
MODELS  CONTAIN  DOF  ON  ROTOR  MAST 


a.  Natural  frequency  comparison 


| 

Mode 

92  D.O.F. 
Unmodified 

92  D.O.F. 
Baseline 

1 

4.18 

4.57 

4.57 

4.57 

4.57 

2 

5.12 

5.12 

5.12 

5.12 

5.12 

3 

6.87 

6.95 

6.95 

6.95 

6.95 

4 

8.02 

8.19 

8.19 

8.19 

8.19 

5 

12.39 

12.41 

12.41 

12.41 

12.41 

6 

15.75 

16.57 

16.57 

16.57 

16.57 

7 

16.37 

17.99 

17.99 

17.99 

17.99 

8 

18.85 

22.25 

22.26 

22.25 

22.24 

9 

20.60 

23.17 

23.17 

23.17 

23.16 

10 

23.93 

24.63 

24.63 

24.63 

24.62 

11 

24.70 

28.00 

27.86 

27.86 

27.87 

12 

27.62 

28.20 

28.18 

28.18 

28.18 

13 

29.92 

30.41 

30.42 

30.42 

30.41 

14 

31.06 

31.21 

31.21 

31.21 

31.21 

15 

33.13 

33.30 

33.30 

33.30 

33.30 

b.  Generalized  mass  comparison 


Inboard, 

Outboard  Wing 

Stores  Removed! 

■  T'FH'Hi—  I'I'H'HI 

Mode 

Unmodified 

Baseline 

1 

1.089 

0.922 

0.922 

0.922 

0.922 

2 

0.842 

0.842 

0.842 

0.842 

0.842 

3 

0.404 

0.417 

0.417 

0.417 

0.417 

4 

0.732 

0.693 

0.693 

0.693 

0.693 

5 

7.524 

7.472 

7.470 

7.471 

7.476 

6 

3.242 

0.818 

0.818 

0.818 

0.818 

7 

0.832 

1.644 

1.645 

1.645 

1.647 

8 

2.831 

0.013 

7.997 

7.986 

7.664 

9 

2.754 

2.273 

2.271 

2.270 

2.256 

10 

4.575 

5.320 

5.328 

5.323 

5.414 

11 

2.751 

2.340 

2.727 

2.698 

2.655 

12 

1.506 

2.701 

1.261 

1.259 

1.259 

13 

2.369 

3.008 

3.086 

3.070 

3.074 

14 

1.796 

2.502 

2.512 

2.503 

2.502 

15 

1.172 

1.154 

1.154 

1.153 

1.156 
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TABLE  11.  INBOARD,  OUTBOARD  WING  STORES  REMOVED 
FIRST  FUSELAGE  LATERAL  BENDING 

ORTHONORMAL  MODE 


92  D.O.F. 

30  D.O.F. 

20  D.O.F. 

1 

-0.0383 

-0.0383 

-0.0383 

2 

4 

-0.0324 

-0.0324 

-0.0324 

5 

0.0117 

0.0117 

0.0117 

6 

-0.1338 

7 

0.0078 

9 

0.0140 

10 

0.1491 

0.1492 

0.1491 

11 

14 

0.4150 

0.4150 

15 

.VKSilTiVlc «!' 

0.0232 

0.0232 

17 

0.6576 

0.6576 

-icisww; ; r,  H 

18 

0.0262 

0.0261 

0.0261 

19 

0.8163 

0.8163 

0.8163 

20 

0.0295 

0.0295 

23 

1.5495 

1.5495 

1.5495 

1.5495 

24 

IS 

0.0343 

0.0342 

0.0341 

26 

-  3  WKci 

-0.0046 

-0.0045 

-0.0045 

27 

HBr  M 

0.0780 

0.0780 

0.0780 

29 

miX.'. 

-0.0108 

-0.0108 

-0.0108 

30 

®  inSSi 

0.1980 

0.1980 

0.1980 

33 

-0.0017 

-0.0017 

35 

0.0029 

0.0029 

36 

0.0131 

0.0131 

37 

-0.0007 

-0.0007 

39 

-0.1439 

-0.1439 

40 

0.0189 

0.0189 

41 

-0.1729 

-0.1729 

42 

-0.0216 

-0.0216 

TABLE  11.  CONTINUED. 


§2  D.O.P.  ~ 
Element  No. 

$2  0.6. F. 

50  D.O.F. 

30  b.O.F 

45 

0.0047 

0.0047 

47 

-0.0775 

-0.0775 

48 

-0.0019 

-0.0019 

50 

-0.0729 

-0.0729 

52 

0.0262 

0.0262 

54 

0.0656 

0.0656 

0.0656 

57 

0.0935 

0.0935 

0.0935 

58 

-0.0673 

-0.0673 

-0.0673 

61 

-0.0956 

-0.0956 

-0.0956 

63 

-0.0256 

-0.0256 

64 

0.0006 

0.0006 

66 

0.0385 

0.0385 

0.0385 

67 

-0.0194 

-0.0194 

69 

0.0218 

0.0218 

72 

0.0527 

0.0527 

74 

0.0629 

76 

0.0115 

0.0629 

79 

0.0188 

0.0188 

80 

0.1350 

0.1350 

0.1350 

81 

0.0189 

0.0189 

83 

0.1601 

0.1601 

84 

-0.0137 

-0.0137 

86 

0.2173 

0.2173 

0.2173 

87 

-0.0094 

-0.0094 

-0.0095 

89 

0.0093 

0.0093 

90 

0.2600 

0.2600 

0.2600 

91 

0.2600 

92 

0.0060 

0.0061 

0.0060 

TABLE  12.  INBOARD  OUTBOARD  WING  STORES  REMOVED 
SECOND  FUSELAGE  VERTICAL  BENDING 
ORTHONORMAL  MODE 


D.O.F. 
merit  No. 


92  D.O.F. 

0.7799 

0.7798 

0.7798 

0.7173 

0.7172 

0.1540 

0.1540 

0.1540 

0.6386 

0.6385 

0.6384 

i 

0.0930 

0.0930 

0.5089 

0.5088 

0.6547 

0.6546 

0.2030 

0.2032 

0.2034 

0.6167 

0.6165 

0.2129 

0.2132 

0.2135 

0.2122 

0.3182 

0.3181 

0.3180 

0.3180 

0.1549 

0.1551 

0.1553 

0.1336 

0.1335 

0.1334 

0.1337 

0.1032 

0.1035 

0.1030 

-0.0883 

0.0885 

-0.1616 

-0.1618 

0.1620 

0.1608 

-0.4108 

-0.4108 

-0.4109 

-0.4104 

-0.0007 

-0.0007 

-0.0008 

-0.0009 

-0.0202 

-0.0202 

-0.0203 

-0.0202 

0.0096 

0.0096 

0.0096 

0.0098 

0.0077 

0.0076 

0.0007 

0.0077 

0.0932 

0.0932 

0.1009 

0.1009 

-0.1084 

-0.1084 

-0.1115 

-0.1115 

-0.0037 

-0.0037 

0.1909 

0.1908 

0.0074 

0.0075 

0.1818 

0.1818 

1 

2 

4 

5 

6 
7 
9 

10 

11 

14 

15 

17 

18 

19 

20 

23 

24 
26 
27 

29 

30 
33 

35 

36 

37 

39 

40 

41 

42 


TABLE  12.  CONTINUED. 


45 

0.3405 

0.3404 

47 

-0.0052 

-0.0052 

48 

-0.0080 

-0.0080 

50 

-0.0142 

-0.0142 

52 

-0.0449 

-0.0449 

54 

-0.1227 

-0.1226 

57 

-0.1413 

-0.1412 

58 

-0.1140 

-0.1139 

61 

-0.1160 

-0.1159 

63 

-0.0187 

-0.0187 

64 

-0.1220 

-0.1221 

66 

-0.0107 

-0.0107 

67 

-0.1077 

69 

-0.1088 

-0.1088 

72 

-0.0001 

-0.0001 

74 

-0.0099 

-0.0099 

76 

-0.0064 

-0.0064 

79 

0.1062 

0.1062 

80 

0.0472 

0.0472 

81 

0.1426 

83 

0.0360 

0.0360 

84 

0.1808 

0.1808 

86 

0.0999 

0.1000 

87 

0.4531 

0.4531 

89 

0.5420 

0.5420 

90 

0.0797 

0.0797 

91 

0.0798 

92 

0.6536 

0.6536 

26 

12 


-0.1139 

-0.1158 


-0.0108 

-0.1077 


-0.1140 

-0.1161 


0.0473 

0.1425 


0.1001 

0.4531 

0.0800 

0.6536 


0.0992 

0.4529 


TABLE  13.  STIFFNESS  CHANGE  CONFIGURATION  -  ADD  STIFFNESS 
8.17  x  104  lb/in.  Between  40145  AND  52045 
MODELS  CONTAIN  DOF  ON  ROTOR  MAST 


a.  Natural  frequency  comparison 


- ! 

With  Stiffness  Change  1 

MTBTTTI  iHf 

92  D.O.F. 

50  D.O.F. 

30  D.O.F. 

20  D.O.F. 

Mode 

Baseline 

1 

4.18 

4.18 

4.18 

4.18 

4.17 

2 

5.12 

5.27 

5.27 

5.27 

5.27 

3 

6.87 

6.87 

6.87 

6.87 

6.87 

4 

8.02 

12.36 

12.36 

12.36 

12.36 

5 

12.39 

12.86 

12.86 

12.86 

12.86 

6 

15.75 

16.31 

16.31 

16.31 

16.31 

7 

16.37 

17.62 

17.62 

17.62 

17.61 

8 

18.85 

20.60 

20.60 

20.60 

20.60 

9 

20.60 

23.57 

23.57 

23.57 

23.60 

10 

23.93 

24.38 

24.38 

24.38 

24.38 

11 

24.70 

27.60 

27.60 

27.60 

27.61 

12 

27.62 

28.88 

28.88 

28.89 

28.89 

13 

29.95 

30.87 

30.87 

30.87 

30.87 

14 

31.06 

32.86 

32.68 

32.68 

32.68 

15 

33.13 

41.34 

41.34 

41.34 

41.34 

b.  Generalized  mass  comparison 


_ With  Stiffness  Change 

Mode 

92  D.O.F. 
Unmodified 

92  D.O.F. 
Baseline 

50  D.O.F. 

30  D.O.F. 

20  D.O.F. 

1 

1.089 

1.089 

1.089 

1.089 

1.089 

2 

0.842 

0.8*  1 

0.820 

0.820 

0.821 

3 

0.404 

0.404 

0.404 

0.404 

0.404 

4 

0.732 

7.396 

7.394 

7.398 

7.408 

5 

7.524 

2.54D 

2.539 

2.551 

2.553 

6 

3.242 

0.789 

0.789 

0.789 

0.789 

7 

0.832 

6.80  ; 

6.796 

6.797 

6.823 

8 

2.831 

2.743 

2.743 

2.743 

14.350 

9 

2.754 

5.650 

5.649 

5.646 

11.350 

10 

4.575 

•1.474 

6.710 

8.493 

12.650 

11 

2.751 

1.511 

4.788 

4.788 

4.847 

12 

1.506 

1  834 

2.971 

7.849 

7.819 

13 

2.369 

1.  354 

6.415 

6.414 

6.401 

14 

1.796 

1 .214 

1.214 

1.214 

1.217 

15 

1.172 

2.678 

3.201 

3.389 

3.393 
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generalized  mass  data,  respectively,  for  truncated  models  having  92,  50, 

30,  and  20  degrees  of  freedom.  Additionally,  modal  data  are  given  for 
the  unmodified  92  degree  of  freedom  model. 

Natural  frequencies  derived  from  the  reduced  degree  of  freedom  models  are 
in  excellent  agreement  with  the  natural  frequencies  predicted  by  the  92 
degree  of  freedom  baseline  model.  However,  the  50  and  30  degree  of  freedom 
models  yield  generalized  masses  which  correlate  with  those  of  the  baseline 
model  for  the  first  eight  modes  and  for  the  fourteenth  mode  while  the 
20  degree  of  freedom  model  generalized  masses  correlate  with  the  first 
seven  modes  and  for  the  fourteenth  mode.  The  generalized  masses  associ¬ 
ated  with  the  remaining  modes  lack  agreement  with  the  baseline  model,  which 
indicates  that  the  coordinates  used  in  the  reduced  degree  of  freedom  models 
lack  adequate  representation  of  the  modes  in  question. 

Table  14  presents  first  fuselage  vertical  bending  modal  data  for  the  afore¬ 
mentioned  stiffness  change  configuration  for  the  50,  30,  and  20  degree  of 
freedom  models  and,  additionally,  the  corresponding  modal  elements  for  the 
92  degree  of  freedom  baseline  model.  The  Modal  Correlation  Number  for  the 
first  fuselage  vertical  mode  associated  with  each  reduced  model  is  also 
shown  in  Table  14  to  effect  a  comparison  between  modal  correlation  and  the 
value  of  the  Modal  Correlation  Number.  A  review  of  the  data  given  in  the 
table  indicates  excellent  correlation  of  the  modal  elements  obtained  from 
the  reduced  models  compared  to  the  baseline  92  degree  of  freedom  modal 
elements. 

Table  15  presents  data  for  a  combination  mass  and  stiffness  change.  The 
condition  that  is  described  is  for  the  clean  wing  configuration  with  a  strut 

4 

type  diagonal  element  having  stiffness  of  8.17  x  10  lb/in.  added  between 
grid  points  40145  and  52045.  Tables  15a  and  15b  show  natural  frequency  and 
generalized  mass  predictions,  respectively,  which  were  obtained  from  the  92, 
50,  30  and  20  degree  of  freedom  truncated  models  in  which  the  aforementioned 
mass  and  stiffness  change  combination  was  incorporated.  Each  of  the  re¬ 
duced  degree  of  freedom  models  predicted  extremely  accurately  the  natural 
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TABLE  14.  STIFFNESS  CHANGE  CONFIGURATION  -  ADD  STIFFNESS 
8.17  x  10^  Ib/in  between  40145  and  52045 
FIRST  FUSELAGE  VERTICAL  BFNDING 
ORTHONORMAL  MODE 


snnrr 

Element  No. 

4 

5 
8 
9 

10 

11 

14 

15 

17 

18 
19 
22 

23 

24 

25 

29 

30 

35 

36 

37 

39 

40 

41 

42 

45 

46 

47 


W OX 


-0.0141 

0.4178 

-0.0097 

0.4956 

0.0034 

0.5613 

0.0198 

0.6264 

0.0343 

0.6275 

0.0377 

0.0178 

0.0364 

0.5986 

-0.0890 

-0.0051 

0.0142 

-0.0207 

-0.0870 

-0.0886 

-0.0128 

0.0170 

-0.0129 

0.0304 

0.1184 

-0.0142 

-0.0144 


50  D.rt.F. 


-0.0141 

0.4178 

-0.0097 

0.4956 

0.0034 

0.5613 

0.0198 

0.6264 

0.0343 

0.6275 

0.0377 

0.0178 

0.0364 

0.5986 

-0.0890 

-0.0051 

0.0142 

-0.0207 

-0.0870 

-0.0886 

-0.0128 

0.0170 

-0.0129 

0.0304 

0.1184 

-0.0142 

-0.0144 


30  D.O.F. 


-0.0147 

0.4176 

-0.0097 

0.4953 

0.0028 

0.0191 

0.6260 


0.0178 

0.0372 

0.5983 

-0.0889 

-0.0051 

0.0143 


0.0886 


-0.0130 

0.0305 


20  D.O.F. 


-0.0152 

0.4177 

-0.0097 

0.4952 

0.0021 

0.0184 

0.6259 


0.0178 

0.0383 

0.5982 

-0.0049 

0.0139 


-0.0131 

0.0307 


-0.0145 


TABLE  14.  CONTINUED 


“snnrr 

Element  No. 

— m — 

“56  6.6.P. 

"  So  D.O.F. 

jTff  6.6.P. — 

48 

-0.0679 

-0.0679 

-0.0678 

50 

-0.0090 

-0.0090 

52 

-0.0857 

-0.0857 

54 

-0.1840 

-0.1840 

-0.1840 

-0.1841 

57 

-0.2671 

-0.2671 

-0.2670 

-0.2673 

58 

-0.1670 

-0.1670 

-0.1669 

-0.1670 

61 

-0.2260 

-0.2260 

-0.2258 

-0.2260 

63 

-0.0269 

-0.0269 

64 

-0.0266 

-0.0766 

66 

0.0057 

0.0057 

0.0057 

67 

-0.0285 

-0.0186 

69 

0.0083 

0.0083 

72 

0.0058 

0.0058 

74 

0.0165 

0.0165 

76 

0.1332 

0.1331 

79 

0.2319 

0.2319 

80 

0.0475 

0.0475 

0.0478 

81 

0.2603 

0.2602 

83 

0.0586 

0.0586 

84 

0.2471 

0.2471 

86 

0.0948 

0.0948 

0.0954 

87 

0.4375 

0.4375 

0.4372 

89 

0.5319 

0.5319 

90 

0.1083 

0.1083 

0.1090 

0.1089 

91 

0. 1084 

0.1084 

0.1091 

0.1090 

[  Modal  Correlati 

on 

Number 

6.168  x  10"6 

3.674  x  10"' 

*  7.131  x  10"4 

74 


TABLE  15. 


a. 


COMBINATION  MASS,  STIFFNESS  CHANGE 
MASS  CONFIGURATION  -  CLEAN  WING 
STIFFNESS  CONFIGURATION  -  ADD  STIFFNESS 

8.17  x  104  LB/IN.  BETWEEN  40145  AND  52045 
30  AND  20  DOF  MODELS  CONTAIN  NO  DOF  ON  ROTOR  MAST 

Natural  frequency  comparison 


Mass/Stiffness  Change  1 

92  D.O.F. 

92  D.O.F. 

50  D.O.F. 

30  D.O.F. 

20  D.O.r . 

Mode 

Unmodified 

Baseline 

■H| 

4.18 

4.57 

4.57 

4.57 

4.57 

5.12 

5.27 

5.27 

5.27 

5.27 

6.87 

6.95 

6.95 

6.95 

6.95 

8.02 

12.39 

12.39 

12.39 

12.39 

«881|pM 

12.39 

13.87 

13.87 

13.87 

14.03 

15.75 

16.58 

16.58 

16.58 

16.58 

16.37 

22.18 

22.18 

22.18 

22.74 

18.85 

23.17 

23.17 

23.17 

23.77 

20.60 

24.62 

24.62 

24.62 

25.83 

23.93 

26.64 

26.64 

26.64 

m3: 

24.70 

28.18 

28.18 

28.19 

H  Mg- 

•27.62 

30.13 

30.13 

30.13 

»» 

29.95 

30.91 

30.91 

30.91 

Wi 

31.06 

32.87 

32.87 

32.87 

*TB§| 

33.13 

41.97 

41.97 

41.97 

b.  Generalized  mass  comparison 


— 

_ Mass/Stiffness  Change 

Mode 

92  D,0.F, 
Unmodified 

92  D.O.F. 
Baseline 

50  D.O.F. 

30  D.O.F, 

20  D.O.F. 

1 

1.089 

0.921 

0.921 

5.985 

5.970 

2 

0.842 

0.820 

0.820 

38.950 

39.590 

3 

0.404 

0.416 

0.416 

0.416 

0.417 

4 

0.732 

7.223 

7.221 

7.222 

7.196 

5 

7.524 

1.920 

1.920 

1.921 

2.102 

6 

3.242 

0.802 

0..802 

0.802 

0.813 

0.832 

9.416 

9.420 

9.419 

6.311 

2.831 

2.275 

2.274 

2.275 

2,758 

2.754 

4.992 

4.982 

8.582 

2.119 

4.575 

3.067 

3.932 

3.930 

2.751 

1.293 

9.430 

9.399 

II 

1.506 

3.896 

4.650 

4.641 

|| 

2.369 

1.403 

6.692 

11.51 

1.796 

1.198 

1.199 

1.198 

H 

1.172 

_ 

2.797 

3.463 

3.462 
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frequencies  resultln  from  the  combination  of  mass  and  stiffness  change. 

The  generalized  masses,  which  are  dependent  upon  the  coordinates  Incor¬ 
porated  In  the  model,  show  a  variation  In  correlation  depending  upon 
the  mode  under  consideration.  For  example,  the  30  degree  of  freedom  model 
was  not  able  to  successfully  predict  the  generalized  masses  associated  with 
the  first  and  second  mode.  This  occurred  because  the  model  was  deficient  in 
coordinates  located  on  the  rotor  mast  and  the  first  and  second  modes  are 
the  pylon  pitch  and  pylon  roll  modes,  respectively.  Modal  element  data  Is 
presented  in  Table  16,  and  a  comparison  of  the  elements  for  the  50,  30,  and 
20  degree  of  freedom  models  with  the  data  from  the  92  degree  of  freedom 
model  reveals  excellent  correlation. 
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TABLE  16.  COMBINATION  MASS,  STIFFNESS  CHANGE 
MASS  CONFIGURATION  -  CLEAN  WING 
STIFFNESS  CONFIGURATION  -  ADD  STIFFNESS 
8.17  x  10^  Ib/in,  between  40145  and  52045 
FIRST  FUSELAGE  VERTICAL  BENDING 
ORTHONORMAL  MODE 


mwn 

i 

0.6382 

0.6382 

IHHI 

4 

0.0425 

0.0425 

.0625 

5 

0.4622 

0.4622 

Mil  |  '.H 

.4511 

6 

0.0119 

0.0118 

i  Wr 

7 

0.2636 

0.2636 

|g|  JR  X  9§ 

8 

-0.0278 

-0.0278 

-.0217 

9 

0.5575 

0.5575 

0.5575 

.5385 

0.0771 

0.0770 

0.0771 

11 

0.6386 

0.6386 

0.6386 

14 

0.1006 

0.1005 

0.1005 

.1002 

15 

0.7199 

0.7200 

0.7199 

.6815 

0.0875 

0.0874 

0.0875 

.0830 

18 

0.7216 

0.7216 

0.7216 

.6898 

22 

0.0058 

0.0058 

0.0058 

-.0007 

23 

-0.0485 

-0.0485 

-.0591 

24 

0.6861 

■ml 

.6599 

25 

-0.1392 

rusk 

29 

-0.0065 

-0.0065 

.30 

-0.0060 

-0.0060 

31 

-0.1402 

-0.1401 

33 

-0.0630 

-0.0630 

35 

-0.0622 

-0.0622 

37 

-0.1377 

-0.1377 

39 

-0.0111 

-0.0111 

40 

-0.0131 

-0.0131 

-0.0131 

-.0020 

41 

-0.0077 

-0.0077 

-0.0077 

.0047 

42 

-0.0057 

-0.0057 

-0.0130 

43 

-0.0070 

-0.0070 

46 

-0.0131 

-0.0131 

TABLE  16.  CONTINUED. 


— snnrr — 

Element  No. 

113991 

47 

-0.0130 

-0.0130 

-0.0130 

48 

-0.1096 

-0.1096 

-0.1096 

50 

-0.0116 

-0.0116 

52 

-0.1244 

-0.1244 

54 

-0.1511 

-0.1511 

-0.1511 

-.1498 

57 

-0.1661 

-0.1661 

-0.1661 

-.1708 

58 

-0.1306 

-0.1306 

-0.1306 

-.1344 

61 

-0.1278 

-0.1277 

-0.1278 

-.1380 

63 

-0.0168 

-0.0168 

0.0168 

64 

-0.1165 

-0.1165 

-0.1165 

66 

0.0003 

0.0003 

67 

-0.0389 

-0.0389 

72 

0.0065 

0.0065 

74 

0.0047 

0.0047 

76 

0.1249 

0.1249 

0.1249 

.1141 

77 

0.0206 

0.0206 

.0219 

30 

0.0354 

c ;  WHi 

81 

0.2785 

j 

.2917 

86 

Q.Gu71 

mBWS  .- *i§5 

0.0671 

87 

Q. 5214 

0.5213 

0.5214 

89 

0.6282 

0.6282 

90 

0.0667 

0.0668 

INCOMPLETE  MODEL 


The  Incomplete  Model  methodology  provides  a  technique  for  correcting  a 
dynamic  analytical  model  of  a  structure  based  on  natural  frequencies  and 
mode  shapes  developed  from  ground  vibration  tests  of  the  structure.  The 
analysis  depends  on  a  set  of  natural  frequencies  and  mode  shapes  which 
the  analytical  models  are  required  to  predict  and  mass  and  stiffness 
analytical  models  which  are  representative  of  the  structure  such  that 
small  changes  In  the  elements  of  each  will  yield  the  desired  characters 
tics. 

MASS  MATRIX  IMPROVEMENT 


The  approach  taken  in  the  Incomplete  Model  methodology  is  that  the  meas¬ 
ured  modes  of  a  structure  are  correct  and  the  analytical  mass  matrix  is 
modified  to  satisfy  the  orthogonality  relationship  among  the  measured 
modes  and  the  modified  mass  matrix. 


Consider  M^  an  analytical  mass  matrix  having  n  degrees  of  freedom  and  y  a 
matrix  containing  m  measured  orthonormal  modes  each  having  n  degrees  of 
freedom,  where  m<n.  The  degrees  of  freedom  must  correspond  to  the  points 
of  measurement  on  the  structure.  In  the  situation  considered  herein  MA 
corresponds  to  the  92  degree  of  freedom  mass  matrix  developed  from  the 
NASTRAN  model  of  the  AH-1G  airframe  and  the  modal  matrix  consists  of  modes 
developed  from  the  accelerations  measured  at  locations  on  the  airframe 
corresponding  to  the  same  points  of  interest  as  used  in  the  mass  matrix. 

If  AM  represents  the  changes  in  the  mass  matrix  required  to  satisfy  the 
orthogonality  relationship,  then 


fT  (M A  +  aM)v  =  I 
or 

>i'T  aMy  =  I  -  ma 

Q 

Where  m^  is  the  nondiagonal  ^M^f. 


(69) 

Equation  (69)  has  an  infinite  number 


of  solutions;  that  is,  there  are  an  infinite  number  of  aM  matrices  that 
will  satisfy  the  system.  It  is  possible  to  determine  a  aM  which  has  a 
minimum  Euclidean  norm  subject  to  the  constraints  established  in  equation 
(69).  Following  the  procedure  used  in  References  6  and  7,  the  function  to 
be  minimized  is  established  as 

e  *  | |Ma'1/2  am  Ma‘1/2||  (70) 

Minimization  of  equation  (70)  results  in  a  AM  matrix  which  is  a  minimum 
in  the  weighted  Euclidean  sense. 


Using  Lagrangian  multipliers,  A.,,  the  constraints  of  equation  (69)  may  be 

*  J 

applied  to  the  function  to  be  minimized  (equation  70).  Thus,  the  follow¬ 
ing  Lagrangian  function  may  be  written: 

m  m  T  1 

3  =  e  +  E  E  A. .  $'aM4>  -  I  +  m.  (71 ) 

i=l  J=1  1J  AJij 

If  equation  (71)  is  differentiated  with  respect  to  each  element  of  the  AM 
matrix  and  the  results  are  set  equal  to  zero  the  minimization  of  equation 
(70)  is  effected  within  the  constraints  of  equation  (69).  This  process 
results  in  a  matrix  of  partial  derivatives  given  as 


dam. 


=  2M 


-1 


AM  M 


-1 


+  ¥  A 


YT  =  0 


1J 


(72) 


Solving  equation  (72)  for  the  aM  matrix  yields 
aM  =  -  1/2  Ma  nT  4,T  Ma 


(73) 


6.  Berman,  A. ,  et  al.,  IMPROVEMENT  OF  ANALYTICAL  DYNAMIC  MODELS  USING 
MODAL  TEST  DATA.  AAIAA/ASME/ASCE/AHS  21st  Structures,  Structural 
Dynamics  and  Materials  Conference,  May  12-14,  1980,  Seattle,  Washington. 

7.  Baruch,  M. ,  and  Bar  Itzhock,  I.Y.,  OPTIMAL  WEIGHTED  0RTHAG0NALIZATI0N 
OF  MEASURED  MODES.  AIAA  Journal,  vol.  16,  No.  4,  1978,  pp.  346-351. 
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where  a  is  a  square  matrix  of  order  m  whose  elements  a^j  are  the  Lagran- 
gian  multipliers.  Substituting  equation  (73)  into  equation  (69)  yields 
the  solution  for  a 

A  =  -2  mA_1  (I  -  mA)mA"1  (74) 

When  equation  (74)  is  substituted  into  equation  (73),  the  corrected  mass 
matrix  is  given  by 

M  =  Ma  +  aM  =  Ma  +  Ma  v  mA_1  (I  -  mA)mA_1  yT  MA  (75) 

The  corrected  mass  matrix,  which  is  symmetrical,  is  consistent  with  the 
measured  modes.  The  only  matrix  inversion  required  in  the  process  is 
that  of  mA,  which  is  of  the  order  of  the  measured  modes. 

STIFFNESS  MATRIX  IMPROVEMENT 


In  a  similar  manner,  it  is  possible  to  correct  the  stiffness  matrix 
using  the  natural  frequencies  and  mode  shapes  developed  from  test  re 
suits,  the  improved  mass  matrix,  and  the  eigenvalue  equation  for  the 
ith  mode 


[K  -  ^2M]  'F.j  =  0 


(76) 


In  addition,  three  constraint  equations  are  imposed  on  the  improved 
stiffness  matrix,  including  the  eigenvalue  equation  for  all  the  measured 
modes,  the  orthogonality  relationship  for  the  stiffness  matrix,  and  the 
requirement  for  a  symmetrical  stiffness  matrix.  The  constraint  equations 
are  given  by 


Kh*  -  Mtsj  =  0 


T„  2 

f  K<F  -  SJ 


K  -  K1  =  0 


(77) 

(78) 

(79) 
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The  function  to  be  minimized  is 


(80) 


Lagrangian  multipliers  a,.  ,  a.  and  ac  may  be  used  in  conjunction 

Kij  °ij  sij 

with  the  stiffness  matrix  constraints  of  equations  (77),  (78)  and  (79), 
respectively.  The  resulting  Lagrangian  function  is  given  by 


n 

m 

6  =  e  +  E 

i  =  l 

E 

j  =  l 

Ak 

n 

m 

+  E 

E 

A 

i  =  l 

j  =  l 

0 

ij 


(81) 


■I  hi  T  o  n  n  j 

+  E  EX  ( ¥  KW  -  “  )  ^  ■;  +  E  E  X.  (K-K1)., 

i=l  j=l  °ij  1J  i=l  j=l  sij  1J 


Differentiating  equation  (81)  with  respect  to  the  elements  of  the  K  matrix 
and  setting  the  result  to  0  yields 

2M~1aKM'"1  +  a/  +  4'X0  Y7  +  A$  -  XST  =0  (82) 

Solving  equation  (82)  for  the  required  change  in  the  stiffness  matrix 
gives 


aK  =  K  -  Ka  =  1/2  M  [xK  fT  +  nQ  y1  +  Xs  -  x$T]  (83) 

Substituting  for  K  in  Equation  (79)  yields  a  solution  for  the  quantity 
(xs  -  xj)  which,  when  substituted  into  equation  (83),  eliminates  xs  arid 
results  in 


K  Ka  -  1/4  M  [nKT  +  A^1)  M  -  1/4  M^[A0  +  AQThM  (84) 

Inserting  equation  (84)  into  the  constraint  equation  (78)  and  using 
=  I  gives 

(aq  +  xT)  =  4(4,TKah'  -  n2)  -  (x^M'y  +  yTMxK)  (85) 
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Equation  (85)  may  be  substituted  into  equation  (84)  and  the  result  is 
K  =  K A  -  ( 4'TKa»f  -  Q2)^  -  l/4MnKT(I  -  MnT)M 

-  1/4M(I  -  (86) 

This  equation  is  substituted  into  the  constraint  equation  (77),  giving 

M(I  -  yyTM)xK  =  (I  -  MwT)Kay  (87) 

and  its  transpose  is  given  by 

akT(I  -  MyvT)M  =  yTKA ( yyTM  -  I)  (88) 

The  quantities  on  the  left  side  of  equations  (87)  and  (88)  appear  in 
equation  (86);  these  substitutions  result  in  the  corrected  stiffness 
matrix  expressed  as 

K  =  Ka  -  (KA4,yTM  +  MwTKa)  +  MT(*V  +  (89) 
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MASS  AND  STIFFNESS  CHANGES 


An  important  criterion  for  evaluating  the  effectiveness  of  an  analytical 
model  developed  from  dynamic  test  data  is  Its  ability  to  predict  the 
effect  of  mass  and  stiffness  changes  to  the  structure.  To  judge  the  use¬ 
fulness  of  the  Truncated  and  Incomplete  Model  methodologies,  analytical 
models  were  developed  from  dynamic  test  data  for  the  mass  and  stiffness 
configurations  that  were  tested.  A  truncated  model  and  an  incomplete 
model  were  developed  for  the  high  gross  weight,  forward  center-of-gravity 
configuration  and  were  used  as  the  baseline  analytical  models.  Mass  and 
stiffness  changes  were  made  to  the  analytical  models  and  the  results,  in 
the  form  of  modal  data,  were  compared  to  similar  data  extracted  from  test 
results  for  the  corresponding  airframe  configuration.  Figure  10  presents 
a  flow  diagram  representing  the  typical  process  which  was  used  in  apply¬ 
ing  the  Truncated  Model  methodology  regarding  mass  and  stiffness  changes. 

A  flow  diagram  depicting  implementation  of  the  Incomplete  Model  method¬ 
ology  for  mass  and  stiffness  modifications  to  the  AH-1G  airframe  is  shown 
in  Figure  11 . 

MASS  CHANGES 

Three  gross  weight  configurations  were  tested  during  the  program.  However, 
since  the  high  gross  weight,  forward  center-of-gravity  condition  was  used 
as  the  baseline,  this  allowed  two  mass  change  conditions  to  be  implemented 
analytically.  The  two  mass  changes  were: 

1.  High  gross  weight,  forward  center  of  gravity  to  mean  gross 

weight: 

a.  remove  inboard  wing  stores 

b.  add  forward  fuel,  3S7  lb 

c.  add  aft  fuel ,  276  lb 

2.  High  gross  weight,  forward  center  of  gravity  to  low  gross 

weight: 

a.  remove  inboard,  outboard  wing  stores 

b.  add  forward  fuel,  429  lb 

c.  add  aft  fuel ,  354  lb 
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Figure  10.  Truncated  model  configuration  change. 
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Figure  11.  Incomplete  model  configuration  change. 
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The  wing  stores  were  considered  to  be  concentrated  masses  rigidly  attached 
to  the  wings  at  degrees  of  freedom  corresponding  to  the  NASTRAN  model  of 
the  AH-1G  and  at  which  points  accelerometers  were  located  on  the  AH-1G  air¬ 
frame.  The  addition  of  fuel  to  the  forward  and  aft  fuel  tanks  in  the  test 
article  was  also  treated  as  the  addition  of  concentrated  masses  at  grid 
points  in  the  NASTRAN  model  which  corresponded  to  accelerometers  positioned 
geometrically  proximate  to  the  fuel  tank  locations  on  the  AH-1G  airframe. 
Since  the  fuel  is  a  continuum,  consideration  of  fuel  changes  as  discrete 
impedance  changes  is  extremely  approximate.  To  investigate  this  effect 
many  computer  experiments  were  performed  related  to  the  manner  of  fuel 
and  wing  stores  discretization.  Tables  17,  18,  and  19  represent  typical 
data  used  for  computer  runs  for  mass  changes  in  developing  the  high  gross 
weight  to  low  gross  weight  configuration.  Similarly,  Tables  20,  21,  and 
22  contain  data  used  for  generating  the  mean  gross  weight  configuration 
from  the  high  gross  weight  forward  center  of  gravity  baseline  model. 


TABLE  17.  MASS  CHANGE  -  HIGH  GROSS 
WEIGHT  TO  LOW  GROSS  WEIGHT 


DOF 

NASTRAN 

GRID  IDENT. 

MASS  CHANGE 
LB-SEC2/IN. 

DESCRIPTION 

47 

22703T2 

.917 

Add  AFT  Fuel 

48 

22725T3 

.917 

Add  AFT  Fuel 

54 

64229T3 

-1.370 

Remove  Inboard  Wing  Store,  R.S 

57 

65929T3 

-  .997 

Remove  Outboard  Wing  Store.  R.S 

58 

74229T3 

-1.370 

Remove  Inboard  Wing  Store,  L.S 

61 

75929T3 

-  .997 

Remove  Outboard  Wing  Store,  L.S 

66 

15625T2 

1.110 

Add  Forward  Fuel 

67 

15629T3 

.555 

Add  Forward  Fuel 

69 

15621T3 

.555 

Add  Forward  Fuel 

TABLE  18.  MASS  CHANGE  -  HIGH  GROSS 
WEIGHT  TO  LOW  GROSS  WEIGHT 


DOF 

NASTRAN 

GRID  IDENT. 

MASS  CHANGE 
LB-SEC2/IN- 

DESCRIPTION 

47 

22703T2 

.917 

Add  AFT  Fuel 

48 

22725T3 

.917 

Add  Aft  Fuel 

54 

64229T3 

-1.370 

Remove  Inboard  Wing  Store,  R.S 

57 

65929T3 

-  .997 

Remove  Outboard  Wing  Store,  R.S 

58 

74229T3 

-1.370 

Remove  Inboard  Wing  Store,  L.S 

61 

75929T3 

-  .997 

Remove  Outboard  Wing  Store,  L.S 

62 

18645T1 

1.110 

Add  Forward  Fuel 

63 

18635T2 

.4575 

Add  Forward  Fuel 

64 

18625T3 

.4575 

Add  Forward  Fuel 

66 

15625T2 

.6535 

Add  Forward  Fuel 

67 

15629T3 

.3268 

Add  Forward  Fuel 

69 

15621T3 

.3268 

Add  Forward  Fuel 

TABLE  19.  MASS  CHANGE  -  HIGH  GROSS 
WEIGHT  TO  LOW  GROSS  WEIGHT 


DOF 

NASTRAN 

GRID  IDENT. 

MASS  CHANGE 

lb-sec2/in. 

DESCRIPTION 

47 

22703T2 

.917 

Add  AFT  Fuel 

48 

22725T3 

.917 

Add  Aft  Fuel 

55 

65011T1 

-2.367 

Remove  Inboard,  Outboard  Wing 
Stores,  L.S 

56 

65011T3 

-2.367 

Remove  Inboard,  Outboard  Wing 
Stores,  L.S 

59 

75011T1 

-2.367 

Remove  Inboard,  Outboard  Wing 
Stores,  R.S 

60 

7501 1T3 

-2.367 

Remove  Inboard,  Outboard  Wing 
Stores,  R.S 

62 

18645T1 

1.110 

Add  Forward  Fuel 

63 

18635T2 

.4575 

Add  Forward  Fuel 

64 

18625T3 

.4575 

Add  Forward  Fuel 

66 

15625T2 

.6535 

Add  Forward  Fuel 

67 

15629T3 

.3268 

Add  Forward  Fuel 

69 

15621T3 

.3268 

Add  Forward  Fuel 
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TABLE  20.  MASS  CHANGE  -  HIGH  GROSS 

WEIGHT  TO  MEAN  GROSS  WEIGHT 


DOF 

NASTRAN 

GRID  IDENT. 

MASS  CHANGE 
LB  SEC2/IN. 

DESCRIPTION 

47 

22703T2 

.715 

Add  AFT  Fuel 

48 

22725T3 

.715 

Add  Aft  Fuel 

54 

64229T3 

-1.370 

Remove  Inboard  Wing  Store,  R.S 

58 

74229T3 

-1.370 

Remove  Inboard  Wing  Store,  L.S 

66 

15625T2 

.925 

Add  Forward  Fuel 

67 

15629T3 

.4625 

Add  Forward  Fuel 

69 

15621T3 

.4625 

Add  Forward  Fuel 

TABLE  21.  MASS  CHANGE  -  HIGH  GROSS 

WEIGHT  TO  MEAN  GROSS  WEIGHT 


NASTRAN 

GRID  I DENT. 

MASS  CHANGE 
LB  SEC2/IN. 

DESCRIPTION 

22703T2 

.715 

Add  AFT  Fuel 

22725T3 

.715 

Add  AFT  Fuel 

64229T3 

-1.370 

Remove  Inboard  Wing  Store,  R.S 

74229T3 

-1.370 

Remove  Inboard  Wing  Store,  L.S 

18645T1 

.925 

Add  Forward  Fuel 

18635T2 

.381 

Add  Forward  Fuel 

18635T3 

.381 

Add  Forward  Fuel 

15625T2 

.544 

Add  Forward  Fuel 

15629T3 

.272 

Add  Forward  Fuel 

15621T3 

.272 

Add  Forward  Fuel 

TABLE  22.  MASS  CHANGE  “  HIGH  GROSS 

WEIGHT  TO  MEAN  GROSS  WRIGHT 


DOF 

NASTRAN 

GRID  IDENT. 

MASS  CHANGE 
LB  SEC2/IN. 

47 

22703T2 

.715 

48 

22725T3 

.715 

55 

65011T1 

-1.370 

56 

65011T3 

-1.370 

59 

75011T1 

-1.370 

60 

75011T3 

-1.370 

62 

18645T1 

.925 

63 

18635T2 

.381 

64 

18635T3 

.381 

66 

15625T2 

.544 

67 

15629T3 

.272 

69 

15621T3 

.272 

DESCRIPTION 


Add  AFT  Fuel 
Add  AFT  Fuel 

Remove  Inboard  Wing  Store,  L.S 

Remove  Inboard  Wing  Store,  L.S 

Remove  Inboard  Wing  Store,  R.S 

Remove  Inboard  Wing  Store,  R.S 

Add  Forward  Fuel 
Add  Forward  Fuel 
Add  Forward  Fuel 
Add  Forward  Fuel 
Add  Forward  Fuel 
Add  Forward  Fuel 


COMPARISON  OF  ANALYTICAL  AND  TEST  RESULTS 


INCOMPLETE  MODEL  -  FREQUENCY  RESPONSE 


One  measure  of  the  success  of  an  analytical  dynamic  model  is  its  ability 
to  predict  the  dynamic  response  of  a  structure  compared  to  test  measure¬ 
ments.  The  acceleration  mobility,  or  accelerance,  reflects  the  adequacy 
of  the  analytical  mass  and  stiffness  matrices  in  duplicating  test  data 
over  a  frequency  band  of  interest.  The  acceleration  mobility  matrix  for 
a  linear  system  is  given  as 


[M  -  K]  -  i  a  [K] 

C  (l) 

U) 


-1 


(90) 


The  individual  elements  of  the  acceleration  mobility  matrix  characterize 
the  response  at  the  ith  degree  of  freedom  for  force  excitation  at  the 
jth  degree  of  freedom. 


The  Incomplete  Model  methodology  was  applied  to  the  92  degree  of  freedom 
NASTRAN  model  to  develop  improved  mass  and  stiffness  matrices  representa¬ 
tive  of  the  AH- 16  airframe.  Subsequently,  equation  (90)  was  used  to 
determine  the  acceleration  mobility  matrices  of  the  system.  In  evaluating 
the  acceleration  mobility  of  the  structure  at  a  particular  forcing  fre¬ 
quency  the  value  of  the  structural  damping  coefficient,  g,  must  be  speci¬ 
fied.  For  frequencies  of  excitation  proximate  to  the  natural  frequencies 
the  modal  damping  coefficient  obtained  from  test  measurements  was  used  in 
the  analysis.  However,  use  of  the  modal  damping  coefficients  are  theor¬ 
etically  valid  only  at  the  system  natural  frequencies.  In  the  current 
program  at  forcing  frequencies  off-resonance  in  the  range  of  9-13  Hz,  a 
linear  distribution  of  structural  damping  was  applied  based  on  the  modal 
damping  coefficients  associated  with  the  natural  frequencies  at  7.32  Hz 
and  14.35  Hz.  Additionally,  for  comparative  purposes  values  of  structural 
damping  of  .05  and  .10  were  used  for  the  same  frequency  band.  Equation  (90) 
was  used  to  develop  acceleration  mobility  matrices  for  a  frequency  spectrum 
of  5  Hz  to  22  Hz.  Typical  frequency  response  curves  are  presented  in 
Figures  12,  13,  14,  and  15  showing  comparisons  between  the  analytical 
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Figure  12.  Acceleration  mobility  Y^q  ^485 
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Figure  14.  Acceleration  mobility  Yy^i  y521 
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model  predictions  incorporating  the  Improved  mass  and  stiffness  matrices 
developed  via  the  Incomplete  Model  theory  and  test  measurements.  Since  the 
mobility  acceleration  Is  a  complex  quantity,  data  Is  presented  in  terms 
of  magnitude  (g's/lb)  and  phase  angle  (degrees).  The  values  of  structural 
damping  which  were  used  in  the  analysis  are  shown  In  Table  23.  Figures 
12a  and  12b  present  analytical  and  test  results  for  the  magnitude  and  phase 
angle,  respectively,  for  the  acceleration  mobility  YZ5Q  2435*  which  repre¬ 
sents  the  acceleration  in  the  vertical  direction  at  fuselage  station  50  In. 
for  force  excitation  applied  in  the  vertical  direction  at  fuselage  station 
485  in.  Figures  13a  and  13b  show  data  in  the  same  format  for  the  acceler¬ 
ation  mobility  Y2202R,Z485’  which  represents  vertical  acceleration  mobility 
for  a  coordinate  on  the  right  wing  of  the  airframe  for  excitation  applies 
at  the  tail  skid  location.  Similarly,  Figures  14  and  15  present  accelera¬ 
tion  mobility  data  for  the  response  stations  indicated  for  force  excitation 
applied  laterally  at  the  tail  rotor  gearbox  location.  A  review  of  the  afore¬ 
mentioned  figures  indicates  good  agreement  between  the  theoretical  predic¬ 
tions  and  the  test  measurements,  particularly  at  frequencies  of  excitation 
proximate  to  the  natural  frequencies  of  the  airframe.  The  predicted  vibra¬ 
tion  response  obtained  from  the  92  degree  of  freedom  NASTRAN  model  is  also 
presented  on  the  aforementioned  figures.  The  data  indicates  that  applica¬ 
tion  of  the  Incomplete  Model  methodology  has  significantly  improved  the 
analytical  model  vibration  response  prediction  capability.  The  natural 
frequencies  predicted  using  the  reduced  NA5TRAN  model  did  not  agree  with 
test  data;  however,  the  incomplete  model  results  showed  a  high  degree  of 
correlation  with  test  data  in  the  vicinity  of  the  natural  frequencies. 
Additionally,  the  incomplete  model  improved  the  vibration  response  predic¬ 
tions  for  the  off  resonance  frequency  bands.  One  factor  which  must  be  con¬ 
sidered  is  that  the  test  data  presented  in  the  aforementioned  figures  re¬ 
flects  the  complex  modes  which  characterize  the  AH-1G  airframe,  while  the 
analytical  predictions  are  based  on  a  mathematical  model  developed  from 
orthonormal  modal  data  consisting  of  scalar  elements.  Introduction  of 
complex  orthonormal  modes  into  the  analysis  would,  undoubtedly,  improve 
the  correlation. 
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TABLE  23.  STRUCTURAL  DAMPING  USED  FOR 

ACCELERATION  MOBILITY  PREDICTION. 


Natural 

Frequency 

Hz 

D 

9 

Frequency 

Ranqe 

Hz 

6.28 

.16 

5.8  -  6.4 

7.32 

.066 

7.0  -  9.0 

14.35 

.18 

13.8  -  14.6 

15.99 

.103 

15.6  -  17.0 

17.60 

.083 

17.2  -  18.4 

22.06 

.148 

_ 

20.0  -  22.4 

TRUNCATED  MODEL  MASS  AND STIFFNESS  CHANGES 


Ms Inc  the  Truncated  Model  methodology,  the  effects  of  mass  changes  on  the 
AH- 1 G  airframe  were  predicted.  The  baseline  model  used  In  the  analysis 
was  based  on  six  modes  which  were  extracted  from  test  data  and  considered 
to  be  significant  and  linearly  independent.  The  analytical  predictions  of 
natural  frequencies  and  orthonormal  mode  shape  elements  resulting  from  the 
mass  changes  were  compared  to  test  results  for  the  respective  mass  config¬ 
urations. 

The  AH-1G  airframe  is  characterized  by  high  modal  response  associated  with 
the  tail  boom,  particularly  in  the  area  of  the  vertical  fin,  and  to  a  lesser 
degree  in  the  nose  section  of  the  fuselage.  Generally,  the  center  section 
of  the  airframe  displays  low  modal  activity. 

MASS  CHANGES 

The  mass  changes  presented  in  Table  17  were  used  in  the  truncated  model  for 
predicting  the  low  gross  weight  natural  frequencies  and  orthonormal  mode 
shapes  from  the  high  gross  weight  baseline  configuration.  Table  24  presents 
a  comparison  of  predicted  natural  frequencies  and  the  results  obtained  from 
shake  tests  of  the  low  gross  weight  configuration. 


TABLE  24.  NATURAL  FREQUENCIES  FOR  TRUNCATED  MODEL 
MASS  CHANGE  -  HIGH  TO  LOW  GROSS  WEIGHT 


Natural  Frequencies 

Test 

Analysis 

6.29 

6.38 

7.19 

7.37 

14.78 

14.52 

17.71 

16.12 

19.14 

18.30 
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A  comparison  of  the  orthonormal  modes  associated  with  the  first  three 
natural  frequencies  shown  In  the  preceding  table  showed  excellent  agree¬ 
ment  between  the  modal  elements  obtained  from  test  data  and  those  developed 
analytically.  However,  a  comparison  of  the  modal  data  developed  from  shake 
tests  of  the  high  gross  weight  configuration  and  the  low  gross  weight  condi¬ 
tion  displayed  minor  changes  in  the  respective  modal  elements.  Table  25 
presents  results  for  the  first  mode  which  Is  the  fuselage  first  lateral 
bending  mode.  Similar  results  were  obtained  for  the  natural  frequencies 
at  7.19  Hz  and  14.78  Hz.  Table  26  shows  a  comparison  of  test  and  analyti¬ 
cal  results  for  the  mode  associated  with  17.71  Hz.  Results  for  the  19.14  Hz 
mode  are  shown  in  Table  27, indicating  the  predictive  capability  of  the 
truncated  model.  Qualitatively,  the  results  are  in  agreement  and  the 
truncated  model  satisfactorily  predicts  the  direction  of  change  in  the 
modal  elements  in  developing  the  low  gross  weight  data  from  the  high  gross 
weight  configuration. 


TABLE  25.  TRUNCATED  MODEL  MASS  CHANGE  -  HIGH 
GROSS  WEIGHT  TO  LOW  GROSS  WEIGHT 
(Orthonormal  Mode  6.29  Hz) 


Natural 

Frequency 

_ 

High  Gross  Weight 

Test 

Low  Gross  Weight 

Test 

Analysis 

6.28 

6.29 

6.38 

D.O 

F. 

Y464L 

14 

.354 

.292 

.379 

Y480 

17 

.669 

.600 

.697 

Y490 

19 

.501 

.539 

.524 

Y517 

21 

1.423 

1.669 

1.441 

Y521 

23 

1 .883 

1.891 

1.902 

YHUB 

30 

.7530 

.650 

.767 

Y95R 

80 

.220 

.149 

.228 

Y50 

90 

.354 

.371 

.368 
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TABLE  26.  TRUNCATED  MODEL  MASS  CHANGE  -  HIGH 
GROSS  WEIGHT  TO  LOW  GROSS  WEIGHT 
(Orthonormal  Mode  17.71  Hz) 


Natural 

Frequency 

High  Gross  Weight 

Test 

15.99 

Low  Gross  Weight 

Test 

17.71 

Analysi s 

16.12 

D.O.F. 

Z396R  2 

.576 

.573 

.511 

Y380  4 

.342 

.307 

.571 

Z34C  7 

.728 

-.198 

.392 

Y480  17 

.430 

.610 

.772 

X540  22 

.904 

.681 

.807 

Z540  24 

-.663 

-.358 

-.561 

Y60L  86 

-.651 

.168 

.070 

Z50  92 

.334 

-.070 

.139 

TABLE  27.  TRUNCATED  MODEL  MASS  CHANGE  -  HIGH 
GROSS  WEIGHT  TO  LOW  GROSS  WEIGHT 
(Orthonormal  Mode  Comparison) 


Natural 

Frequency 

High  Gross  Weight 

Test 

17.60 

Low  Gross  Wei  gilt 

Test 

19.14 

Analysi s 

18.30 

D.O. 

F. 

Z396R 

2 

.190 

.619 

.533 

Z400 

9 

.204 

.580 

.539 

Z464 

15 

.145 

.739 

.480 

Z460 

18 

.125 

.866 

.531 

X540 

22 

.158 

.986 

.706 

Z540 

24 

-.141 

-.362 

1 

ro 

00 

00 

Z93R 

81 

.186 

-.0280 

.201 

Z46R 

89 

.297 

-.124 

.337 
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Similarly,  the  mass  changes  indicated  in  Table  17  were  applied  to  the 
high  gross  weight  mass  configuration  to  develop  the  natural  frequencies 
and  orthonormal  mode  shapes  associated  with  the  mean  gross  weight  condi¬ 
tion.  Table  28  shows  the  predicted  natural  frequencies  and  test  results 
for  the  mean  gross  weight  vehicle.  The  analytical  data  were  achieved  by 
using  the  high  gross  weight  truncated  model  and  incorporating  analytically 
the  mass  changes  to  effect  the  mean  gross  weight  condition. 


TABLE  28.  NATURAL  FREQUENCIES  FOR  TRUNCATED  MODEL  MASS  CHANGE  - 
HIGH  TO  MEAN  GROSS  WEIGHT 


Natural 

Frequencies 

Test 

Analysis 

6.17 

6.30 

7.15 

7.34 

14.39 

14.28 

17.88 

16.03 

24.12 

22.45 

For  the  first  two  modes  given  in  Table  28  the  mean  gross  weight  modal  data 
generated  analytically  from  the  high  gross  weight  baseline  truncated  model 
were  in  agreement  with  the  similar  data  developed  from  dynamic  tests  of 
the  mean  gross  weight  configuration.  However,  the  orthonormal  modes 
developed  from  the  high  gross  weight  test  data  did  not  vary  significantly 
from  the  respective  modes  obtained  from  the  mean  gross  weight  test  results. 
Table  29  presents  a  comparison  of  the  modal  elements  for  the  14.39  Hz  mode. 
The  analysis  results  display  the  correct  trend  although,  generally,  the 
changes  from  the  high  gross  weight  baseline  configuration  are  small. 

Table  30  shows  elemental  mode  data  for  the  24.12  Hz  mode. 

Comparison  of  the  mode  shapes  developed  analytically  for  the  low  and  mean 
gross  weight  configurations  by  introducing  concentrated  mass  changes  to 
the  high  gross  weight  configuration  truncated  model  generally  showed  good 
agreement  with  natural  frequencies  and  the  respective  modes  generated  from 
test  data. 
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TABLE  29.  TRUNCATED  MODEL  MASS  CHANGE  -  HIGH 
GROSS  WEIGHT  TO  MEAN  GROSS  WEIGHT 
(Orthonormal  Mode  14.39  Hz) 


Natural 

Frequency 

High  Gross  Weight 

Test 

14.35 

Mean  Gross  Weight 

Test 

14.39 

Analysis 

14.28 

D.O 

F. 

Z396L 

3 

.159 

-.211 

.151 

Y422L 

10 

>.817 

-.543 

-.780 

Y440 

13 

-1.148 

-1.064 

-1 .093 

Y480 

17 

-1.021 

-.624 

-.960 

Y521 

23 

.376 

.669 

.380 

Y95R 

80 

-.349 

.198 

-.345 

Y60L 

86 

-.545 

-.351 

-.547 

Y50 

90 

-.676 

-.423 

-.665 

TABLE  30.  TRUNCATED  MODEL  MASS  CHANGE  -  HIGH 


GROSS  WEIGHT  TO  MEAN  GROSS  WEIGHT 
(Orthonormal  Mode  24.12  Hz) 


Natural 

Frequency 

High  Gross  Weight 

Test 

22.06 

Mean  Gross  Weight 

Test 

24.12 

Analysis 

22.45 

D.O. 

F. 

Z396R 

2 

.482 

.309 

.438 

Z400 

9 

.543 

.328 

.509 

Z423 

11 

.649 

.432 

.617 

Z464 

15 

.709 

.378 

.679 

Y490 

19 

.226 

.189 

.210 

Y517 

21 

.090 

.195 

.131 

Z61R 

88 

-.170 

-.093 

-.197 

Z50 

92 

-.247 

-.188 

-.286 
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STIFFNESS  CHANGES 


The  Truncated  Model  theory  was  used  to  predict  the  effect  of  the  three 
stiffness  changes  which  were  implemented  during  the  ground  vibration  tests 
of  the  AH-1G  airframe.  Three  strut  type  stiffness  changes  were  applieJ. 

The  Truncated  Model  developed  from  the  high  gross  weight,  forward  center- 
of-gravity  configuration  was  used  as  the  baseline  and  the  stiffness  changes 
were  analytically  incorporated  in  that  model. 

During  the  ground  vibration  tests  of  the  stiffness  modifications  force 
excitation  was  applied  in  the  vertical  direction  only.  The  results  of 
the  shake  tests  indicated  that  removal  of  the  lift  link  did  not  signifi¬ 
cantly  alter  the  natural  frequencies  or  response  of  the  airframe  vertical 
bending  modes.  Removal  of  the  lift  link  effectively  separates  the  upper 
body,  consisting  of  the  transmission  and  rotor  pylon,  from  the  basic  airframe; 
consequently,  the  airframe  modal  characteristics  are  not  affected.  A  rotor 
pylon  mode  in  the  longitudinal  direction  was  observed  in  the  test  results 
but  did  not  appear  in  the  analytical  results.  This  is  probably  due  to 
inadequate  representation  of  the  pylon  structure  in  the  original  analytical 
model . 

The  first  stiffness  change  involved  analytically  deleting  the  lift  link 
calculated  stiffness  (2.717  x  10^  Ib/in.)  between  coordinates  Z195T  and 
Z195B.  Table  31  lists  the  natural  frequencies  developed  theoretically  and 
those  obtained  from  test  data. 
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TABLE  31.  NATURAL  FREQUENCIES  FOR  TRUNCATED  MODEL 


STIFFNESS  CHANGE  -  LIFT  LINK  REMOVED 


Natural 

Frequencies 

Test 

Analysis 

7.16 

7.32 

14.65 

14.35 

20.84 

22.06 

Table  32  presents  data  for  the  first  vertical  bending  mode  at  7.16  Hz. 
Similarly,  the  data  given  in  Table  33  is  for  the  14.65  Hz  mode.  Test 
results  for  the  lift  link  removed  configuration  revealed  a  mode  at  20.84  Hz 
with  primary  response  in  the  vertical  direction.  Analysis  indicated  a  mode 
at  22.06  Hz,  also  in  the  vertical  direction.  However,  there  was  little 
correlation  between  the  respective  modal  elements. 

The  second  stiffness  change  that  was  implemented  in  the  airframe  involved 
substitution  of  a  DAVI  lift  link,  with  stiffness  of  7500  lb/iri.  in  the 
vertical  direction,  for  the  conventional  lift  link.  Regarding  the  analy¬ 
tical  model,  this  change  was  reflected  by  deleting  stiffness  of  2.709  x 
107  Ib/in.  between  coordinates  Z195T  and  Z195B, which  are  the  approximate 
lift  link  attachment  points. 
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TABLE  32.  TRUNCATED  MODEL  STIFFNESS  CHANGE  - 
LIFT  LINK  REMOVED 
(Orthonormal  Mode  7.16  Hz) 


Lift  Link  In 
Test 


Lift  Link  Out 


Test 


Analysis 


Natural 

Frequency 


D.O.F. 


7.32 


7.16 


7.32 


Z404R  1 
Z464  15 
Z460  18 
Z485  20 
Z540  24 
ZMHUB  29 
X60L  87 


.139 

.407 

.439 

.610 


.145 

.477 

.512 

.722 

1.177 

.316 

.228 


.130 

.399 

.430 

.602 

.967 

.305 

.177 


.974 

.314 

.174 


TABLE  33.  TRUNCATED  MODEL  STIFFNESS  CHANGE  - 
LIFT  LINK  REMOVED 
(Orthonormal  Mode  14.65  Hz) 


Lift  Link  In 

Lift  Link  Out 

Natural 

Frequency 

Test 

Test 

Analysis 

14.35 

14.65 

14.35 

D.O.F. 

Y380 

4 

.478 

.343 

.471 

Y339L 

6 

.246 

.249 

.247 

Y422L 

10 

.817 

.487 

.810 

Y464L 

14 

1.073 

.667 

1.061 

Y490 

19 

.917 

.495 

.909 

Y520 

23 

-.376 

-.358 

-.375 

Y60L 

86 

.545 

.120 

.549 

Y50 

90 

.676 

.154 

.675 

This  value  is  the  difference  in  spring  rate  between  the  original  lift  link 
(2.717  x  107  Ib/in.)  and  the  spring  rate  of  the  DAVI  lift  link  (7500  Ib/in.) 
Table  34  presents  the  predicted  natural  frequencies  and  test  values  for  the 
DAVI  lift  link  configuration. 


TABLE  34.  NATURAL  FREQUENCIES  FOR  TRUNCATED  MODEL  STIFFNESS 
CHANGE  -  DAVI  LIFT  LINK  CONFIGURATION 


Natural  Frequencies 

Test 

Analysis 

7.17 

7.28 

14.89 

14.35 

15.62 

15.95 

20.90 

22.04 
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During  the  test,  the  DAVI  lift  link  condition  exhibited  a  rotor  pylon  mode 
in  the  fore-and-aft  direction  at  a  frequency  of  5.80  Hz,  which  was  not  evi¬ 
dent  in  the  truncated  model  analysis.  Typical  orthonormal  mode  elements  for 
both  theory  and  test  results  are  presented  in  Tables  35  and  36.  For  the  DAVI 
lift  link  installed  configuration  the  modal  elements  derived  analytically 
are  in  agreement  with  the  respective  test  resilts.  However,  there  is  prac¬ 
tically  no  difference  in  modal  data  for  the  t*o  configurations  tested. 


TABLE  35.  TRUNCATED  MODEL  STIFFNESS  CHANGE  - 
DAVI  LIFT  LINK  INSTALLED 
(Orthonormal  Mode  7.17  Hz) 


Lift  Link  In 

DAVI  Lift  Link  Installed 

Natural 

Test 

Test 

Analysis 

Frequency 

7.32 

7.17 

7.28 

D.O. 

F. 

Z396R 

2 

.104 

.105 

.109 

Z423 

11 

.193 

.212 

.189 

Z464 

15 

.407 

.471 

.399 

Z460 

18 

.439 

.502 

.430 

X540 

22 

-.541 

O 

CM 

10 

• 

1 

-.540 

Z540 

24 

.974 

1.123 

.967 

Z46R 

89 

.201 

.252 

.200 

Z50 

92 

.232 

.288 

.235 
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TABLE  36.  TRUNCATED  MODEL  STIFFNESS  CHANGE  - 
DAVI  LIFT  LINK  INSTALLED 
(Orthonormal  Mode  15.62  Hz) 


Lift  Lin!;  in 

DAVI  Lift  Link  Installed 

Natural 

Frequency 

Test 

Test 

Analysis 

15.99 

15.62 

15.95 

D.O. 

F. 

Z404R 

1 

.676 

.782 

.662 

Z381 

5 

.436 

.540 

.424 

Z423 

11 

.431 

.592 

.418 

Z464 

15 

.30? 

.465 

.295 

Z460 

18 

.361 

.522 

.357 

Z540 

24 

-.663 

-.549 

-.644 

Z61R 

88 

.236 

.285 

.213 

Z50 

92 

.334 

.443 

.300 

i 
! 
1 

The  third  stiffness  condition  that  was  shake  tested  introduced  a  stiffen¬ 
ing  element  between  the  tail  boom  and  the  vertical  tail.  The  spring  rate 

4 

of  the  element  was  calculated  as  5.53  x  10  lb/in.  in  a  direction  coinci- 

i 

dent  with  the  longitudinal  axis  of  the  stiffener.  The  natural  frequencies  < 

obtained  via  test  and  those  analytically  generated  are  detailed  in  Table  j 
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TABLE  37.  NATURAL  FREQUENCIES  FOR  TRUNCATED  MODEL  STIFFNESS 
CHANGE  -  TAIL  STIFFENER  CONFIGURATION 


Natural 

Frequencies 

Test 

Analysis 

7.13 

10.83 

14.66 

14.47 

20.94 

23.20 

The  orthonormal  mode  data  associated  with  the  test  natural  frequencies  of 
7.13  Hz  and  20.94  Hz  are  presented  In  Tables  38  and  39.  For  the  tail 
stiffener  configuration  theoretical  and  test  results  did  not  compare  as 
favorably  as  for  the  other  mass  and  stiffness  conditions  investigated.  A 
possible  explanation  is  that  the  attachment  points  of  the  stiffener  ele¬ 
ment  are  soft',  consequently  the  effective  stiffness  of  the  element  is  modi¬ 
fied  dv  the  loca'  flexibility  of  the  attachment  points.  Thus,  the  effec¬ 
tive  spring  rate  introduced  by  the  stiffener  is  insignificant  and  differs 
from  the  value  used  in  the  analysis,  causing  a  lack  of  correlation  between 
test  data  and  theoretical  results. 


TABLE  38.  TRUNCATED  MODEL  STIFFNESS  CHANGE  - 
TAX'.  STIFFti'ER  INSTALLED 
(Orthononal  Mode  7.13  Hz) 


No  Tail  Stiffener 

Tail  Stiffener  Installed 

Natural 

Frequency 

Test 

Test 

Analysis 

7.32 

7.13 

10.83 

D.O.F. 

Z404R  1 

.139 

.131 

.510 

Z423  11 

.193 

.193 

.416 

Z464  lb 

.407 

.425 

.531 

Z460  18 

.439 

.462 

.588 

Z485  20 

.610 

.644 

.603 

Z540  24 

.974 

.709 

.493 

Z50L  87 

.174 

.205 

.237 

Z50  92 

.232 

.266 

_ 

.399 

- -  - .  ■  .  — i 

TABLE  39.  TRUNCATED  MODEL  STIFFNESS  CHANGE  - 
TAIL  STIFFENER  INSTALLED 
(Orthorormal  Mode  20.94  Hz) 


Natural 

Frequency 

No  Tull  Stiffener 

Tail  Stiffener  Installed 

Test 

Test 

Analysis 

22.06 

20.94 

23.20 

D.O.F. 

Z404R 

1 

.59? 

.556 

.232 

Z381 

5 

.436 

.294 

.195 

Z340 

7 

-.515 

-1.048 

-.871 

Z464 

15 

.709 

.363 

.567 

Z485 

20 

.459 

.174 

.497 

X540 

22 

.792 

.695 

.196 

Z61R 

88 

-.175 

-.168 

-.266 

Z50 

92 

-.247 

-.235 

-.371 

Z50 


INCOMPLETE  MODEL  MASS  AND  STIFFNESS  CHANGES 


MASS  CHANGES 


Using  the  Incomplete  Model  development,  the  effect  of  mass  changes  on  the 
AH-1G  airframe  were  predicted.  The  procedure  described  previously  for  the 
Truncated  Model  approach  was  also  used  in  the  Incomplete  Model  application 
for  predicting  the  effects  of  the  mass  and  stiffness  modifications  to  the 
airframe  which  were  shake  tested.  Table  40  shows  the  theoretically  de¬ 
rived  natural  frequencies  and  those  extracted  from  test  data  for  the  low 
gross  weight  condition. 

Table  41  presents  orthonormal  modal  data  at  a  frequency  of  6.29  Hz,  which 
is  the  first  lateral  bending  mode  of  the  .\1rframe.  Modal  data  for  the 
17.71  Hz  natural  frequency  are  delineated  in  Table  42.  The  tabular  data 
present  evidence  that  the  incomplete  model,  generally,  is  capable  of  pre¬ 
dicting  the  effect  of  mass  change  for  the  low  gross  weight  configuration. 

In  particular,  the  changes  in  the  modal  elements  are  usually  in  the  correct 
direction. 

Table  43  shows  predicted  natural  frequencies  developed  from  the  Incomplete 
Model  theory  and  the  values  obtained  from  shake  tests  of  the  mean  gross 
weight  configuration.  The  theoretical  results  were  obtained  by  applying 
the  appropriate  mass  changes  to  the  NASTRAN  model,  which  was  originally 
improved  on  the  basis  of  the  high  gross  weight  orthonormal  modes  developed 
from  test  data. 
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TABLE  40.  NATURAL  FREQUENCIES  FOR  INCOMPLETE  MODEL 
MASS  CHANGE  -  HIGH  TO  LOW  GROSS  WEIGHT 


Natural  Frequencies 

Test 

Analysis 

6.29 

6.29 

7.19 

7.32 

14.78 

14.34 

17.71 

15.97 

19.14 

20.52 

TABLE  41 .  INCOMPLETE  MODEL  MASS  CHANGE  -  HIGH 
GROSS  WEIGHT  TO  LOW  GROSS  WEIGHT 
(Orthonormal  MoHc  C.2C  Hz) 


Natural 

Frequency 

High  Gross  Weight 

Low  Gross  Weight 

Test 

Test 

Analysis 

6.28 

6.29 

6.28 

D.O, F. 

Y464L 

14 

.354 

.292 

.357 

Y480 

17 

.669 

.600 

.673 

Y4l“'vi 

19 

.501 

.539 

.505 

Y517 

21 

1.423 

i  .669 

1.427 

Y521 

23 

1.883 

1.891 

1.887 

YHUB 

30 

.753 

,650 

.751 

Y95R 

80 

.220 

.149 

.220 

Y50 

90 

.354 

.371 

.354 
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TABLE  42.  INCOMPLETE  MODEL  MASS  CHANGE  -  HIGH 
GROSS  WEIGHT  TO  LOW  GROSS  WEIGHT 
(Orthonormal  Mode  17.71  Hz) 


Natural 

Frequency 

High  Gross  Weight 

Low  Gross  Weight 

Test 

Test 

Analysis 

15.99 

17.71 

15.97 

D.O.F. 

Z396R 

2 

.576 

.573 

.412 

Y360 

4 

.342 

.307 

.573 

Z340 

7 

.728 

-.198 

.341 

Y480 

17 

.430 

.610 

.773 

XS40 

22 

.904 

.682 

.725 

Z540 

24 

-.663 

-.358 

-.527 

Y60L 

80 

-.651 

.168 

-.473 

Z50 

92 

.334 

-.07 

.127 

TABLE  43.  NATURAL  FREQUENCIES  FOR  INCOMPLETE  MODEL  MASS  CHANGE 


MASS  CHANGE  -  HIGH  TO  MEAN  GROSS  WEIGHT 


Natural 

Frequencies 

Test 

Analysis 

6.17 

6.32 

7.15 

7.34 

14.39 

14.29 

17.88 

16.00 

24.12 

19.94 

Data  shown  in  Tables  44  and  45  are  for  the  14.39  Hz  and  24.12  Hz  modes 
which  are  natural  frequencies  obtained  from  dynamic  tests  of  the  mean 
gross  weight  configuration.  From  a  review  of  the  mean  gross  weight  data 
presented  it  is  apparent  that  the  mass  changes  that  were  effected  did 
not  significantly  alter  the  airframe  natural  frequencies  or  associated 
mode  shapes. 

TABLE  44.  INCOMPLETE  MODEL  MASS  CHANGE  -  HIGH 
GROSS  WEIGHT  TO  MEAN  GROSS  WEIGHT 
(Orthonormal  Mode  14.39  Hz) 


Natural 

Frequency 

High  Gross  Weight 

Mean  Gross  Weight 

Test 

Test 

Analysis 

14.35 

14.39 

14.29 

D.O.F. 

Z396L 

3 

.159 

-  .211 

.135 

Y422L 

10 

-  .817 

-  .542 

-  .774 

Y440 

13 

-1.148 

LO 

o 

1 

-1.089 

Y480 

17 

-1.021 

-  .624 

-  .964 

Y521 

23 

.376 

.669 

.356 

Y95R 

80 

-  .349 

.198 

-  .343 

Y60L 

86 

-  .545 

-  .351 

-  .541 

Y50 

90 

-  .676 

CO 

<NJ 

1 

-  .660 
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TABLE  45.  INCOMPLETE  MODEL  MASS  CHANGE  -  HIGH 


GROSS  WEIGHT  TO  MEAN  GROSS  WEIGHT 
(Orthonormal  Mode  24.12  Hz) 


Natural 

Frequency 

High  Gross  Weight 

Mean  Gross  Weight 

Test 

Test 

Analysis 

22.06 

24.12 

19.94 

D.O. 

F. 

Z396R 

2 

.482 

.309 

.266 

Z400 

9 

.543 

.328 

.338 

Z423 

11 

.649 

.432 

.386 

Z464 

15 

.709 

.378 

.365 

Y490 

19 

.226 

.189 

-.076 

Y51  7 

21 

.090 

.195 

-.036 

Z61R 

88 

-.170 

.093 

-.281 

Z50 

92 

-.247 

-.188 

-.283 

STIFFNESS  CHANGES 


Incomplete  Model  system  identification  was  also  used  to  predict  the  effect 
on  natural  frequencies  and  mode  shapes  of  the  stiffness  changes  that  were 
incorporated  in  the  AH-1G  airframe  during  the  ground  vibration  phase  of 
the  program.  Table  46  presents  the  natural  frequencies  for  the  lift  link 
removed  stiffness  configuration. 
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TABLE  46.  NATURAL  FREQUENCIES  FOR  INCOMPLETE 
MODEL  STIFFNESS  CHANGE  - 
LIFT  LINK  REMOVED 


Natural 

Frequencies 

Test 

Analysis 

7.27 

14.34 

21.04 


Tables  47  and  48  present  orthonormal  modal  data  obtained  analytically  and 
similar  data  extracted  from  test  results  for  the  14.65  Hz  mode  and  the 
20.84  Hz  test  mode. 

TABLE  47.  INCOMPLETE  MODEL  STIFFNESS  CHANGE  - 
LIF l  LINK  REMOVED 
(Orthonormal  Mode  14.65  Hz) 


Lift  Link  In 


Test 


14.35 


Natural 

Frequency 


Lift  Link  Removed 

Test 

Analysis 

14.65 

14.35 

.343 

.483 

.249 

.244 

.487 

.822 

.667 

1.082 

.495 

.924 

-  .358 

-  .377 

.120 

.542 

.154 

-  .090 

TABLE  48.  INCOMPLETE  MODEL  STIFFNESS  CHANGE 
LIFT  LINK  REMOVED 
(Orthonormal  Mode  20.84  Hz) 


Lift  Link  In 

Lift  Lin 

k  Removed 

Natural 

Frequency 

Test 

Test 

Analysis 

22.06 

20.84 

21.02 

D.O. 

F. 

Z404R 

1 

.597 

-.223 

-.331 

Z423 

11 

.649 

-.225 

-.109 

Z464 

15 

.709 

-.258 

-.005 

X540 

22 

.792 

.109 

-.204 

Z268L 

42 

-.208 

-.083 

.172 

Z200L 

62 

-.052 

1 

o 

o 

CTi 

.172 

Y95R 

80 

.014 

-.015 

.215 

Z50 

92 

-.247 

-.103 

.159 

Natural  frequency  data  for  the  condition  wherein  the  DAVI  lift  link  was 
substituted  for  the  conventional  lift  link  are  given  in  Table  49. 

TABLE  49.  NATURAL  FREQUENCIES  FOR  INCOMPLETE 
MODEL  STIFFNESS  CHANGE  - 
DAVI  LIFT  LINK  INSTALLED 


Z 


Mode  shape  data  associated  with  the  natural  frequencies  of  7.17  and  15.62 
are  presented  in  Tables  50  and  51  respectively.  In  both  the  lift  link 
removed  and  the  DAVI  lift  link  installed  configurations  the  test  results 
do  not  differ  significantly  from  the  test  data  for  the  conventional  lift 
link  condition.  Although  the  analytical  predictions  of  the  stiffness  changes 
agree  reasonably  well  with  the  respective  test  results, there  is  no  basis  for 
judging  the  effectiveness  of  the  Incomplete  Model  methodology  in  this  situa¬ 
tion. 

TABLE  50.  INCOMPLETE  MODEL  STIFFNESS  CHANGE  - 
DAVI  LIFT  LINK  INSTALLED 
(Orthonormal  Mode  7.17  Hz) 


Natural 

Frequency 

D.O.F. 


Lift  Link  In 
Test 
7.32 


DAVI  Lift  Link  In 


Analysis 

7.27 


Z396R 

2 

.104 

.105 

.075 

Z423 

11 

.193 

.212 

.158 

Z464 

15 

.407 

.471 

.323 

Z460 

18 

.439 

.502 

.334 

X540 

22 

-.541 

-.620 

-.409 

Z540 

24 

.974 

1.123 

.785 

Z46R 

89 

.201 

.252 

.148 

X50 

92 

.232 

.288 

.181 

Natural 

Lift  Link  In 

DAVI  Lift  L 

ink  In 

Test 

Test 

Analysis 

Frequency 

15.99 

15.62 

16.02 

D.O.F. 

Z404R  1 

.676 

.782 

.677 

Z381  5 

.436 

.540 

.437 

Z423  11 

.431 

.592 

.431 

Z464  1 5 

.302 

.465 

.295 

Z460  18 

.361 

.522 

.351 

Z540  24 

-.663 

-.549 

-.686 

Z61R  88 

.236 

.285 

.249 

Z50  92 

.334 

.443 

.354 

Natural  frequency  data  for  the  tail  stiffener  configuration  are  shown  in 
Table  52.  Orthonormal  mode  shape  data  are  presented  in  Tables  53  and  54 
for  the  7.13  Hz  and  20.94  Hz  modes,  respectively.  Analytical  results 


TABLE  52.  NATURAL  FREQUENCIES  FOR  INCOMPLETE 
MODEL  STIFFNESS  CHANGl  - 


TAIL 

STIFFENER  CONFIGURATION 

Natural 

Frequencies 

Test 

Analysis 

7.13 

10.92 

14.66 

14.55 

20.94 

22.24 

TABLE  53.  INCOMPLETE  MODEL  STIFFNESS  CHANGE  - 
TAIL  STIFFENER  INSTALLED 
(Orthornormal  Mode  Comparison) 


No  Tail  Stiffener  Tail  Stiffener  Installed 


Natural 

Frequency 

Test 

Test 

Analysis 

7.32 

7.13 

10.92 

D.O.F. 

Z404R 

1 

.139 

.131 

.590 

Z423 

11 

.193 

.193 

.510 

Z464 

15 

.407 

.425 

.617 

Z460 

18 

.439 

.462 

.678 

Z485 

20 

.610 

.644 

.690 

Z540 

24 

.974 

.707 

.352 

Z60L 

87 

.174 

.205 

.235 

Z50 

92 

.232 

.266 

.373 

TABLE  54.  INCOMPLETE  MODEL  STIFFNESS  CHANGE  - 
TAIL  STIFFENER  INSTALLED 

(Orthonormal  Mode  Comparison) 


— 

No  Tail  Stiffener 

Tail  Stiffener  Installed 

Natural 

Frequency 

Test 

Test 

Analysis 

22.06 

20.94 

22.24 

D.O.F. 

Z404R  1 

.597 

.556 

.349 

Z381  5 

.436 

.294 

.267 

Z340  7 

-.515 

-1.048 

-.785 

Z464  1 5 

.709 

.363 

.477 

Z485  20 

.459 

.174 

.380 

X540  22 

.792 

.695 

.123 

Z61R  88 

-.175 

-  .168 

-.237 

Z50  90 

r*^ 

C\J 

1 

-  .235 

-.302 
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showed  the  first  vertical  bending  mode  at  10.92  Hz,  but  the  test  results 
with  the  tall  stiffener  installed  did  not  reflect  the  change  In  system 
natural  frequency.  It  appears  the  tall  stiffener  was  ineffective  because 
of  the  softness  associated  with  the  mounting  system  of  the  stiffener. 

Generally,  both  the  incomplete  model  and  the  truncated  model  predicted  the 
effect  of  mass  and  stiffness  changes  to  the  airframe.  The  mass  and  stiff¬ 
ness  changes  that  were  incorporated  in  the  airframe  did  not  alter  signifi¬ 
cantly  the  natural  frequencies  and  mode  shapes  associated  with  the  unmodi¬ 
fied  airframe.  However,  the  incomplete  model  and  truncated  model  did  pre¬ 
dict  the  correct  trend  resulting  from  the  changes  to  the  airframe. 
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DONE  AND  HUGHES  EQUATION 


The  technique  of  Done  and  Hughes,  which  is  described  in  Reference  8,  was 
applied  to  determine  the  effect  of  the  mass  and  stiffness  changes  that 
were  implemented  during  the  program.  The  method  is  applicable  to  concen¬ 
trated  mass  changes  and  for  strut  type  stiffness  modifications. 


The  acceleration  mobility  between  coordinates  i  and  j,  resulting  from  the 
nth  discrete  impedance  change,  is  related  to  the  acceleration  mobility  prior 
to  the  change  via  the  following  relationship: 


,.n 

vij 


,.n-1 

Yij 


f..n-l  ..n-1  ..n-1  ..n-1  ) 

Y.  -  Y.  Y-  -  Y. 
irn  1sn  Jrn  ^sn 


..n-1  ..n-1  ..n-1 

+  Y  -  2Y  +  Y 
Ars  r  r  rs  s  s 
nn  nn  nn  nn 


(91) 


where  Yn-1  is  the  mobility  (with  respect  to  the  coordinates  specified  by 

the  subscripts)  prior  to  the  nth  discrete  impedance  change  and  Yn  is  the 
mobility  after  the  nth  discrete  impedance  change 


Ars  = 
n  n 


384  K 


r  s 
n  n 


for  a  stiffness  change  of  Kr 

between  coordinates  r  and  sn  n 

n  n 


Ib/in. 


; —  for  a  mass  change  of  nv.  lb  at  coordin- 

l  r 

rn  ate  rn.  n 


8.  Done,  G.  T.  S.,  and  Hughes,  A.  D.,  The  Response  of  a  Vibrating  Struc¬ 
ture  as  a  Function  of  Structural  Parameters.  Journal  of  Sound  and 
Vibration,  Vol.  38,  No.  2,  1975. 
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when  the  Impedance  change  Is  a  lumped  mass  at  coordinate  r  ,  the  moblli 
ties  containing  the  subscripts  sp  in  equation  (91)  are  all  defaulted  to 
zero 


i.e.,  for  a 

n  n 


yn-1  _  yn-1  _  gn-1 

Yisn  “  Yjs  "  Yrnsn 
n  n  n  n 


?r;  -  o 

n  n 


This  is  a  computational  artifice  tc  make  the  form  of  equation  (91)  the 
same  for  both  stiffness  and  mass  changes. 


The  subscripts  in  equation  (91)  are  defined  as  follows: 

i  and  j  refer  to  the  coordinate  pair  for  which  the  changed 
mobility  is  being  calculated 


r  and  s  refer  to  the  connection  coordinates  of  the  nth 
n  n 

impedance  change  (subject  to  the  artifice  mentioned  above). 


To  compute  equation  (91)  for  the  nth  change,  all  the  mobilities  on  the 
right  must  be  known  after  the  (n-l)th  change.  In  other  words,  if  there 
are  N  discrete  changes  to  be  effected,  one  must  start  with  the  following 
mobil ities: 


Y° 


ir_ 


Y? 

is. 


n 


Y? 

or 


n 


Y° 


Jsn 

v° 

Vn 

Y° 

Vn 

v° 

snsn 


for  n 


1  to  N 


130 


MASS  CHANGES 


The  configuration  changes  of  the  AH-1G  from  high  gross  weight  to  mean 
gross  weight,  and  from  high  gross  weight  to  low  gross  weight  are  repre¬ 
sented  as  simultaneous  discrete  mass  changes.  Table  20  shows  the  coordina- 
ates  and  the  mass  changes  associated  with  the  change  in  configuration  from 
high  gross  weight  to  mean  gross  weight.  Table  17  shows  the  coordinates 
and  mass  changes  associated  with  the  change  in  configuration  from  high 
gross  weight  to  low  gross  weight. 

For  two  pairs  of  coordinates  (Z90R,  Z485  gunner  seat  with  respect  to 
tail  vertical  shaking  location  and  Z140R,  Z485  pilot  seat  with  respect 
to  tail  vertical  shaking  location)  between  6  Hz  and  10  Hz,  acceleration 
mobilities  were  calculated.  The  mobilities  of  the  original  structure  were 
derived  by  the  modal  series  method  using  the  modes  measured  from  the  high 
gross  weight  configuration.  New  mobilities  resulting  from  the  discrete 
mass  changes  of  Tables  17  and  20  were  calculated  using  the  derived  original 
mobilities  and  the  Done  and  Hughes  equation.  Finally,  mobilities  based  on 
the  mode  shapes  of  the  new  configurations  (Mean  gross  weight  and  Low  gross 
weight)  were  also  derived  for  the  same  pairs  of  coordinates.  The  results 
are  presented  graphically  in  Figures  16,  17,  18  and  19. 

The  discretization  of  the  mass  changes  due  to  the  fuel  redistribution  asso¬ 
ciated  with  the  helicopter  configuration  changes  is  approximate  at  best; 
the  continuum  character  of  these  changes  can  be  a  major  inaccuracy  in  the 
application  of  the  Done  and  Hughes  equation. 

STIFFNESS  CHANGES 

Equation  91  was  used  to  predict  the  effect  on  acceleration  mobility  of  the 
stiffness  changes  which  were  incorporated  in  the  airframe  during  the  test 
phase  of  the  program  and  which  have  been  described  previously.  Figures  20 
and  21,  respectively,  present  the  acceleration  mobilities  Yz-j g5T  Z4g5  and 
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MOBILITIES  : 
DERIVED  FROM 
PRE-CHANGE  MOOES 
PREDICTED  USING 
DONE  &  HUGHES 
EQUATION 
/  DERIVED  FROM 
I  POST-CHANGE  MODE: 


FREQUENCY 


Figure  17.  YZ140R#Z485  High  to  mea 
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NOBILITIES: 
DERIVED  FROM  PRE 
CHANGE  MODES 
1  PREDICTED  USING 
■'  DONE  &  HUGHES 
EQUATION 

.  DERIVED  FROM  POS 
CHANGE  MODES 


J 


) 


gross  weight  change 


6  8  10 


FREQUENCY  (HZ) 


Figure  19.  yz140r,z485  High  t0  1ow  gross  weight  chan9e- 
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ACCELERATION  MOBILITY 


v'zi 95T  Zi95j  Tor  the  lift  link  removed  configurations.  Similarly,  data 
for  the  stiffness  condition  wherein  the  DAVI  lift  link  replaced  the  conven¬ 
tional  lift  link  are  shown  in  Figures  22  and  23.  The  acceleration  mobil¬ 
ity  data  predicted  using  the  Done  and  Hughes  equation  are  in  reasonable 
agreement  with  the  synthesized  data.  Additionally,  the  shift  in  natural 
frequency  resulting  from  the  stiffness  modification  is  predicted  correctly. 
The  predicted  acceleration  mobility  data  indicate  that  the  natural  frequency 
of  the  system  occurs  at  a  lower  frequency  which  is  consistent  with  the 
acceleration  mobility  data  derived  using  the  measured  modes  of  the  changed 
configuration. 

Data  for  the  tail  stiffener  installed  configuration  are  shown  in  Figure  24. 
The  results  indicate  that  the  effective  impedance  associated  with  the 
installation  of  the  stiffening  element  is  much  lower  than  the  value  which 
was  calculated  for  the  stiffener.  The  acceleration  mobilities  predicted 
from  the  Done  and  Hughes  expression  are  very  small  compared  to  those  derived 
using  the  orthonormal  modes  of  the  changed  stiffness  configuration. 

A  logical  explanation  for  these  results  is  that  the  stiffening  member 
attachment  points  to  the  airframe  are  soft  and  behave  as  series  connected 
impedances  to  the  stiffening  element.  This  is  a  manifestation  of  the  pro¬ 
blem  of  local  flexibility  occurring  in  the  testing  of  structures.  The 
resultant  stiffness  of  the  member  and  the  local  flexibility  is  lower 
than  either  of  the  individual  impedances. 
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ACCELERATION  MOBILITY 


8  10 


FREQUENCY  (HZ) 


X540.Z485 


Figure  24.  Y 


Installation  of  tail  stiffener 


ERROR  ANALYSIS 


In  order  to  evaluate  the  effect  of  possible  error  in  the  test  measurements 
on  the  Incomplete  Model  technique  for  model  improvement,  it  is  essential  to 
simulate  realistic  errors  in  the  test  data.  In  the  present  analysis,  error 
was  applied  to  the  orthornormal  mode  shapes  to  reflect  possible  errors  in 
the  measurement  of  airframe  accelerations. 

The  types  of  errors  that  were  incorporated  in  the  study  included  uniformly 
distributed  random  percentage  error  between  specified  limits  and  a  con¬ 
stant  specified  percentage  bias  error.  Both  types  of  error  were  applied 
to  the  orthonormal  mode  elements. 

The  simulated  measurement  errors  representative  of  the  current  state  of  the 
art  and  used  in  the  analysis  ranged  to  +  10%  random  error  and  10%  bias  error. 
The  purpose  of  the  error  analysis  was  to  determine  the  sensitivity  of  the 
Incomplete  Model  methodology  to  possible  errors  in  test  measurements.  The 
procedure  that  was  followed  involved  applying  specified  levels  of  random 
and  bias  error  to  the  orthonormal  modes  derived  from  test  measurements. 

Using  the  orthonormal  modes  contaminated  with  the  simulated  error,  the 
Incomplete  Model  theory  was  applied  to  the  92  degree  of  freedom  NASTRAN 
model  yielding  modified  mass  and  stiffness  matrices.  These  matrices  were 
subseq  itly  used  in  developing  acceleration  mobilities  at  particular 
forcing  frequencies.  Frequency  spectrums  for  particular  points  of  interest 
on  the  airframe  were  developed  incorporating  the  simulated  measurement 
errors  and  the  results  compared  to  the  zero  measurement  error  condition. 

Figure  25  presents  typical  acceleration  mobility  data,  YZ5q  2485’  describ~ 
ing  the  vertical  response  at  degree  of  freedom  Z50  for  force  excitation 
applied  rrti<  '  at  the  tail  skid  location  Z485.  Data  for  YY52i  y521 
are  snown  in  Figure  26.  In  each  situation  the  frequency  response  data  are 
shown  for  the  zero  error  condition  and  for  several  levels  and  a  combination 
of  random  and  bias  error  applied  to  the  orthonormal  mode  shape  data. 
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Generally,  the  acceleration  mobility  frequency  response  predicted  using 
the  Incomplete  Model  development  Is  relatively  insensitive  to  the  level 
of  measurement  error  applied  In  the  study.  The  data  Indicate  that  even 
with  the  high  level  of  error  used  in  the  analysis  the  frequency  response 
is  well  defined  and  approximates  the  no  error  results.  The  areas  of  the 
frequency  response  proximate  to  the  natural  frequencies  are  more  sensitive 
to  measurement  error  than  are  the  off-resonance  responses  However,  even 
with  the  high  levels  of  measurement  error  incorporated  in  the  study,  the 
system  natural  frequencies  can  be  easily  identified. 

Introducing  the  aforementioned  levels  of  error  into  the  mode  shapes  which 
were  extracted  from  test  measurements  affected  both  the  mass  and  stiffness 
matrices  developed  from  the  Incomplete  Model  system  identification.  The 
respective  changes  in  the  mass  and  stiffness  matrices  and  their  sensitivity 
to  measurement  error  are  reflected  by  the  root  mean  square  of  the  mass  and 
stiffness  changes  and,  more  significantly,  by  the  ratio  of  the  root  mean 
square  of  the  changes  to  the  root  mean  square  of  the  modified  mass  and 
stiffness  matrices.  The  ratio  of  the  root  mean  square  of  the  mass  changes 
to  the  root  mean  square  of  the  elements  of  the  mass  matrix  is  representa¬ 
tive  of  the  average  change  in  the  individual  elements  to  satisfy  the  incom¬ 
plete  model  constraints.  Similarly,  the  parameter  referenced  to  the  stiff¬ 
ness  matrix  reflects  the  average  change  in  the  elements  of  the  stiffness 
matrix.  Tables  55  and  56  present  data  in  root  mean  square  form  for  the 
mass  and  stiffness  matrices,  respectively.  Additionally,  the  effects  of 
measurement  error  on  selected  elements  of  the  mass  and  stiffness  matrices 
are  given  in  Tables  57  and  58,  respectively.  Also  shown  in  the  tables  are 
the  elemental  mess  and  stiffness  values  for  the  92  degree  of  freedom 
NASTRAN  model.  The  data  indicate  that  the  Incomplete  Model  methodology 
yields  mass  and  stiffness  matrices  that  are  insensitive  to  measurement 
error,  at  least  for  the  levels  of  error  considered  in  the  analysis. 
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TABLE  55.  EFFECT  OF  MEASUREMENT  ERROR  ON 
INCOMPLETE  MODEL  MASS  MATRIX 


%  Error 

Random  Bias 

RMS  Of  Mass 
Change 

RMS  Of  Mass 

RMS  Of  Mass 
Change/Mass 

0 

0 

.02490 

.2665 

.0934 

5 

5 

.02063 

.2665 

.0774 

10 

0 

.02317 

.2665 

.0869 

10 

10 

.02027 

.2665 

.0761 

TABLE  56.  EFFECT  OF  MEASUREMENT  ERROR  ON 

INCOMPLETE  MODEL  STIFFNESS  MATRIX 


%  Error 

Random  Bias 

RMS  Of  Stiff¬ 
ness  Change 

RMS  Of 
Stiffness 

RMS  Of  Stiffness 
Change/Stiffness 

0  0 

4.2678  +  04 

2.7919  +  05 

.1529 

5  5 

4.3006  +  04 

2.7919  +  05 

.1540 

10  0 

4.4255  +  04 

2.7919  +  05 

.1585 

10  10 

4.1885  +  04 

2.7919  +  05 

.1500 
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TABLE  57.  EFFECT  OF  ERROR  ON  HASS  ELEMENTS 


Mass 

Element 

Coordinates 

92  D.O.F. 

NASTRAN 

Model 

Incomple 

te  Model 

0%  Random 
0%  Bias 

5%  Random 
5%  Bias 

10%  Random 
0%  Bias 

10%  Random 
10%  Bias 

(23.23) 

.9744 

1.1647 

1.0621 

1.0771 

1.0808 

(25.25) 

3.6783 

3.6761 

3.6743 

3.6657 

3.6813 

(35.35) 

2.3522 

2.3023 

2.2901 

2.2959 

2.2955 

(56.56) 

8.0584 

8.1556 

7.8237 

7.7976 

8.1717 

(60.60) 

8.1160 

8.1505 

8.0465 

8.0618 

8.1170 

TABLE  58.  EFFECT  OF  ERROR  ON  STIFFNESS  ELEMENTS 


Stiffness 

Element 

Coordinates 

92  D.O.F. 

NASTRAN 

Model 

Incomplete  Model 

0%  Random 
0%  Bias 

5%  Random 
5%  Bias 

10%  Random 
0%  Bias 

10%  Random 
10%  Bias 

(23,23) 

6.562+04 

1 . 226+05 

1 . 2676+05 

1 .349+05 

1.328+05 

(25.25) 

4.036+06 

4.083+06 

4.0828+06 

4.091+06 

4.082+06 

(35.35) 

2.107+05 

2.267+05 

2.3417+05 

2.243+05 

2.468+05 

(56.56) 

1.34] ;05 

1.196+06 

1.1681+06 

1.147+06 

1 . 1 58+06 

(60.60) 

1.327+05 

3.288+05 

3.0633+05 

2.962+05 

2.888+05 
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APPLICATION  OF  FORCE  DETERMINATION  FORCuS 

The  Force  Determination  methodology  is  described  In  detail  In  Reference  2; 
a  brief  description  of  the  technique  is  presented  here. 

Using  Force  Determination,  the  magnitudes  and  phasings  of  the  vibratory 
forces  and  moments  acting  on  a  helicopter  can  be  obtained  for  the  flight 
conditions  of  the  mission  profile  and  the  calculated  loads  can  be  reproduced 
on  a  helicopter  during  ground  vibration  tests  in  a  test  bay.  "Ground  flying" 
allows  fatigue  data  to  be  accumulated  at  a  rapid  rate,  and  because  the 
airframe  is  tested  on  the  ground  conditions  pertinent  to  fatigue  damage 
are  observable.  Further,  "ground  flying"  provides  the  opportunity  to 
conduct  vibration  surveys,  to  monitor  acoustic  emission,  and  to  evaluate 
vibration  reduction  devices  and  structural  design  changes. 


Force  Determination  methodology  was  recently  applied  to  a  Bell  AH-1G  air¬ 
frame;  to  implement  the  technique,  34  strain  gage  accelerometers  were 
positioned  on  the  airframe  to  measure  vertical,  lateral,  and  longitudinal 
accelerations.  The  aircraft  was  flown  through  typical  mission  profiles 
for  each  of  four  gross  weight  configurations,  and  the  airframe  acceler¬ 
ometer  signals  were  recorded  yielding  a  vector  of  accelerations  y  for 
each  flight  condition.  Subsequently,  the  aircraft  was  subjected  to  ground 
vibration  tests  using  both  random  and  sinusoidal  force  excitation.  The 
acceleration  and  force  signals  were  processed  using  a  Fast  Fourier  Trans¬ 
form  aralyzer  to  yield  a  calibration  matrix  relating  the  acceleration 
at  point  i  for  excitation  applied  at  point  j.  The  acceleration  mobility 
calibration  matrix  is  given  by 


3V8Fi 

8VSF2  •  ■ 

SV9F1 

8V3F2  •  • 

3V9F1 


.  3yi/3F 


3yj/3F 


9V9Fm 
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The  calibration  matrix  consists  of  N  rows,  one  for  each  accelerometer 
location  and  M  columns,  one  for  each  position  on  the  airframe  at  ,<hich 
a  force  is  acting.  In  the  subject  program,  the  main  rotor  vertical, 
lateral  and  longitudinal  forces  and  a  lateral  force  at  the  tail  rotor 
location  were  considered. 

The  acceleration  response  of  a  structure  at  a  particular  frequency  of 
excitation  is  expressed  by  the  product  of  the  acceleration  mobility  matrix 
and  the  external  forces  acting  on  the  structure 

(y>  =  [Y]  m  (93) 


The  accelerations  and  forces  are,  in  general,  complex  quantities  consist¬ 
ing  of  magnitude  and  relative  phase  angle.  Therefore,  the  acceleration 
mobility  matrix  consists  of  complex  elements,  also.  Theoretically,  if 
there  are  M  forces  to  be  determined  it  is  necessary  to  measure  only  M 
acceleration  on  the  airframe.  However,  the  accelerations  and  the  acceler¬ 
ation  mobility  matrix  in  equation  (93)  are  measured  parameters  and  are 
subject  to  error.  Consequently,  the  solution  of  equation  (93)  for  forces 
could  yield  meaningless  results. 


To  eliminate  this  effect  the  property  of  redundancy  is  introduced  wherein 
many  more  accelerations  are  measured  on  the  airframe  compared  to  the 
forces  which  are  to  be  determined.  As  a  result,  the  calibration  matrix 
of  equation  (93)  becomes  rectangular,  having  many  more  rows  representing 
the  measured  airframe  accelerations  than  columns  which  are  associated 


with  the  forces  to  be  evaluated.  The  resulting  system  of  equations  can 
be  solved  by  using  the  pseudoinverse,  which  effects  a  multivariate  least 
squares  solution  to  the  problem.  Thus,  the  solution  to  equation  (93) 


becomes 

or 


IF)  =  [Y]+  y 
{F }  =  (CY]T[Y]]*1[Y]T  (y> 


(94) 
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In  equation  (94),  [Y]  is  the  calibration  acceleration  mobility  matrix 
that  is  obtained  from  analysis  of  ground  vibration  test  data,  (y }  is  the 
airframe  acceleration  data  acquired  from  flight  tests,  and  is  the 
forces  acting  on  the  airframe  that  are  to  be  determined.  Figure  27  is 
a  pictorial  representation  summarizing  the  Force  Determination  process. 

In  the  current  program  the  forces  that  were  determined  using  the  Force 
Determination  methodology  of  Reference  2  were  applied  to  the  improved 
92  degree  of  freedom  NASTRAN  model  of  the  AH-1G  airframe,  which  was 
improved  via  the  Incomplete  Model  methodology.  Forces  were  applied 
to  the  improved  model  at  a  frequency  of  10.8  Hz,  which  is  the  main  rotor 
2/rev  source  of  excitation.  The  analysis  considered  main  rotor  vertical, 
lateral  and  longitudinal  forces  in  addition  to  a  lateral  force  at  the 
tail  rotor  gearbox  to  be  acting  on  the  airframe.  Table  59  presents  a 
summary  of  the  magnitudes  and  phases  of  the  vibratory  forces  acting  on 
the  airframe  for  the  six  flight  conditions  investigated. 


The  steady-state  equations  of  motion  in  matrix  form  of  a  linear  system 
with  proportional  structural  damping  may  be  expressed  as 


CM] 


UtM 

2 


[K] 


{y}  =  (F) 


(95) 


At  any  forcing  frequency,  w,  the  acceleration  response  and  the  force 
excitation  are  related  by  the  system  impedance  matrix 


[Z]  (y>  =  (F) 


(96) 


From  equation  (95)  the  impedance  matrix  is  given  by 

[Z]  =  [M]  -  Oiial  [K]  (97) 

OJ 

The  acceleration  response  at  a  particular  frequency  of  excitation  is 
given  by 

{ y  }  =  [Z]"1  ^F>  (98) 
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AN  AIRCRAFT  IS  FLOWN  THROUGH  A  TYPICAL  MISSION  PROFILE. 
FIXED  SYSTEM  MOUNTED  ACCELEROMETER  SIGNALS  ARE  RECORDED 
NO  OTHER  INSTRUMENTATION  IS  REQUIRED. 

FROM  THIS  IN-FLIGHT  ACCELERATIONS  ARE  OBTAINED,  OR  {V} 

//////// 


SIGNAL 

PROCESSING 


THE  AIRCRAFT  IS  CALIBRATED  BY  MEANS  OF  A  SHAKE  TEST. 
ACCELEROMETER  AND  FORCE  SIGNALS  ARE  PROCESSED  TO  YIELD 
MEASURED  ACCELERATION  MOBILITY  MATRIX  OF  THE  AIRCRAFT. 


CALCULATE  MAGNITUDE  AND  PHASES  OF  FORCES  AND  MOMENTS  ACTING 
ON  THE  AIRCRAFT  IN  FLIGHT  FROM  {y}  AND  [Y]  + 


Figure  27.  Force  determination  process. 
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Equation  (98)  was  used  in  applying  the  forces  presented  in  Table  59  to 
the  improved  92  degree  of  freedom  NASTRAN  model.  The  impedance  matrix  in 
equation  (98)  consists  of  the  improved  mass  and  stiffness  matrices  developed 
using  the  Incomplete  Model  methodology. 

At  a  particular  frequency  of  excitation  the  impedance  matrix  and,  therefore, 
the  acceleration  response  of  the  system  depends  on  the  structural  damping 
coefficient. 

In  calculating  the  forced  response  of  the  AH-1G  airframe  at  the  2/ rev, 
10.8-Hz,  excitation  frequency  of  the  main  rotor,  a  structural  damping 
coefficient  of  .1224  was  used  which  results  from  the  linear  variation  of 
damping  between  the  second  and  third  natural  frequencies  of  the  system  as 
given  in  Table  23.  The  results  of  applying  the  main  rotor  and  tail  rotor 
forces  obtained  from  Reference  2  for  the  six  flight  conditions  considered 
herein  are  presented  in  Tables  60  through  65.  The  data  are  shown  for  the 
37  accelerometers  used  in  the  program  and  are  identified  with  respect  to 
the  coordinate  definition  shown  in  Table  2.  Additionally,  the  magnitude 
and  phase  angle  of  the  applied  forces  and  their  respective  coordinates  on 
the  airframe  are  presented. 

Acceleration  data  for  both  analysis  and  flight  are  also  presented  in 
Figures  28  through  33.  Generally,  the  analysis  correctly  predicts  the 
distribution  of  acceleration  magnitude  and  phase  angle  throughout  the 
airframe.  The  accelerations  predicted  for  the  coordinate  on  the  vertical 
tail  pylon  do  not  agree  with  the  test  results.  A  contributing  factor  is 
that  the  location  of  the  accelerometers  on  the  top  of  the  pylon  do  not 
coincide  with  the  geometrical  location,  in  terms  of  fuselage,  water  line, 
and  butt  line  stations  of  the  respective  degrees  of  freedom  in  the  NASTRAN 
model . 


D  NASTRAN  MODEL 


TABLE  60.  ACCELERATION  RESPONSE  FOR  MAXIMUM  STRAIGHT  AND 
LEVEL  FLIGHT  AT  9500  LB  GROSS  WEIGHT 


Frequency  =  10.8  Hz 

Structural  Damping  =0.122  I 

Applied 

Force  (lb) 

_ 

Direction 

Real 

Imaginary  1 

Lateral  at  Tail  Rotor  Gearbox 

-73. 

-36. 

Longitudinal  at  hub 

439. 

-85. 

Lateral  at  hub 

-278. 

-187. 

Vertical  at  hub 

319. 

1113. 

Accelerometer 

Magnitude  (q’s)  i 

Phase  (deg)  1 

Identification 

Analysis 

Flight 

Analysis 

Flight 

Z50 

.236 

0.200 

177. 

166. 

Z100T 

.102 

0.139 

192. 

154. 

Z210T 

.086 

0.191 

133. 

161 

Z340 

.105 

0.309 

91. 

40. 

Z400 

.168 

0.210 

60. 

29. 

Z460 

.143 

0.105 

150. 

304. 

Z540 

.246 

0.494 

240. 

235. 

Z90R 

.114 

0.078 

166. 

145. 

Z90L 

.105 

0.120 

118. 

136. 

Z140R 

.042 

0.056 

111. 

104. 

Z140L 

.131 

0.101 

66. 

110. 

Z200R 

.163 

0.171 

59. 

4. 

Z200L 

.394 

0.283 

50. 

81. 

Z260R 

.177 

0.182 

73. 

63. 

Z260L 

.197 

0.222 

64. 

53. 

Z396R 

.373 

0.122 

50. 

77. 

Z396L 

.062 

0.260 

171. 

23. 

ZLONG 

.092 

0.094 

63. 

65. 

ZLATR 

.122 

0.107 

60. 

71. 

ZCOLL 

.160 

0.121 

68. 

80. 

Y50 

.248 

0.037 

233. 

246. 

Y30 

.149 

0.022 

232. 

162. 

Y140 

.100 

0.037 

225. 

144. 

Y220B 

.087 

0.036 

35. 

111. 

Y220T 

.523 

0.389 

214. 

232. 

Y300 

.246 

0.206 

48. 

120. 

Y380 

.308 

0.363 

60. 

155. 

Y440 

.310 

0.628 

75. 

172. 

!  Y490 

.370 

0.474 

185. 

200. 

Y517 

1.891 

1.123 

'  214. 

221. 

X140 

.068 

0.027 

23. 

86. 

X180T 

1.244 

0.272 

172. 

141. 

X540 

.280 

0.444 

245. 

227. 

X200R 

.141 

0.142 

27. 

48. 

X200L 

.092 

0.114 

217. 

250. 

XI 90 

.081 

0.115 

36. 

70. 

X200L 

_ 

.047 

0.027 

23. 

333. 

1R/1 


TABLE  61.  ACCELERATION  RESPONSE  FOR  A  45-DEGREE 
TURN  AT  9500  LB  GROSS  WEIGHT 


1  frequency  =  10.8  Hz 

Structural  Damping 

=  6.122 

!  Applied  Force  (lb)  1 

1  Direction 

Real 

Imaginary  I 

1  Lateral  at  Tail  Rotor  Gearbox 

-  50. 

-  14. 

Longitudinal  at  hub 

139. 

485. 

Lateral  at  hub 

-221. 

-143. 

Vertical  at  hub 

443. 

385. 

Accelerometer 

i  Magnitude  (g's) 

i  Phase  (deg)  1 

Identification 

Analysis 

Flight 

Analysis 

Flight  i 

Z50 

.239 

0.226 

116. 

175. 

Z100T 

.126 

0.120 

103. 

165. 

Z210T 

.062 

0.204 

96. 

154. 

Z340 

.038 

0.236 

292. 

85. 

Z400 

.102 

0.172 

310. 

357. 

Z460 

.044 

0.031 

134. 

234. 

Z540 

.168 

0.449 

252. 

184. 

Z90R 

.089 

0.088 

123. 

184. 

Z90L 

.125 

0.120 

90. 

149. 

Z140R 

.012 

0.029 

71. 

172. 

Z140L 

.079 

0.079 

43. 

127. 

Z200R 

.092 

0.171 

323. 

319. 

Z200L 

.226 

0.195 

25. 

86. 

Z260R 

.072 

0.116 

10. 

27. 

Z260L 

.099 

0.157 

22. 

32. 

Z396R 

.180 

0.129 

31. 

351. 

Z396L 

.167 

0.184 

254. 

351. 

ZLONG 

.045 

0.052 

355. 

40. 

ZLATR 

.063 

0.072 

16. 

57. 

ZCOLL 

.094 

0.076 

35. 

67. 

Y50 

.248 

0.019 

232. 

183. 

Y30 

.142 

0.020 

231 . 

180. 

Y140 

.096 

0.024 

230. 

147. 

Y220B 

.092 

0.046 

45. 

94. 

Y220T 

.396 

0.328 

211. 

224. 

Y300 

.195 

0.156 

42. 

112. 

Y380 

.145 

0.269 

19. 

144. 

Y440 

.093 

0.474 

359. 

161. 

Y490 

.282 

0.483 

211. 

174. 

Y51 7 

1.425 

0.595 

213. 

208. 

XI 40 

.057 

0.029 

312. 

350. 

X180T 

1.386 

0.333 

109. 

144. 

X540 

.101 

0.418 

231 . 

176. 

X200R 

.120 

0.143 

36. 

50. 

X200L 

.138 

0.095 

245. 

277. 

XI 90 

.050 

0  105 

343. 

90. 

X200L 

.040 

0.028 

311. 

325. 

155 


TABLE  62.  ACCELERATION  RESPONSE  FOR  SIDEWARD 
FLIGHT  AT  9500  LB  GROSS  WEIGHT 


Structural  Damping  =  0.122 


led  Force 


Direction 


Lateral  at  Tail  Rotor  Gearbox 
Longitudinal  at  hub 
Lateral  at  hub 
Vertical  at  hub 


Accelerometer 

Identification 


Z50 

Z100T 

Z210T 

Z340 

Z400 

Z460 

Z540 

Z90R 

Z90L 

Z140R 

Z140L 

Z200R 

Z200L 

Z260R 

Z260L 

Z396R 

Z396L 

ZLONG 

ZLATR 

ZCOLL 

Y50 

Y30 

Y140 

Y220B 

Y220T 

Y300 

Y380 

Y440 

Y490 

Y517 

XI 40 

X180T 


Magnitude  (g's 


Analysis  Flight 


.154 

.114 

.129 

.267 

0.162 

0.003 

0.295 

0.103 

0.113 

0.093 

0.111 

0.111 

0.234 

0.198 

0.220 

0.151 

0.145 

0.112 

0.118 

0.138 

0.101 

0.059 

0.047 

0.028 

0.144 

0.056 

0.079 

0.135 

0.081 

0.061 

0.029 

0.191 

0.298 

0.064 

0.034 

0.100 

0.026 


Phase 


Analysis 


TABLE  63.  ACCELERATION  RESPONSE  FOR  APPROACH 
AND  LANDING  AT  9500  LB  GROSS  WEIGHT 


Structural  Damoina 

1 1*111 1  BBWB 

1  Applied  Force  (lb)  1 

Direction 

Real 

Imaginary  ii 

Lateral  at  Tail  Rotor  Gearbox 
Longitudinal  at  hub 

Lateral  at  hub 

Vertical  at  hub 

-  6. 

-  51. 

-137. 

1051. 

2. 

249. 

107. 

382. 

Accelerometer 

Identification 

Mi  li  1 1  ■trn— —  1 

Analysis 

■nunn 

■nrms 

Z50 

.097 

0.085 

100. 

122. 

Z100T 

.065 

78. 

75. 

Z210T 

0.078 

43. 

133. 

Z340 

.091 

0.277 

0. 

14. 

Z400 

.157 

0.187 

350. 

9. 

Z460 

.122 

0.023 

5. 

200. 

Z540 

.151 

0.504 

236. 

188. 

Z90R 

.047 

0.021 

58. 

71. 

Z90L 

.099 

0.074 

49. 

84. 

Z140R 

.045 

0.035 

23. 

35. 

Z140L 

.101 

0.073 

26. 

60. 

Z200R 

.172 

0.111 

0. 

315. 

Z200L 

.276 

0.195 

19. 

45. 

Z260R 

.175 

0.162 

13. 

16. 

Z260L 

.173 

0.198 

17. 

23. 

Z396R 

.145 

0.173 

25. 

354. 

Z396L 

.114 

0.227 

319. 

347. 

ZLONG 

.080 

0.083 

15. 

20. 

ZLATR 

.097 

0.095 

21. 

38. 

ZCOLL 

.141 

0.105 

23. 

37. 

Y50 

.107 

0.058 

256. 

2. 

Y30 

.068 

0.026 

242. 

358. 

Y140 

.043 

0.013 

240. 

18. 

Y220B 

.036 

0.024 

98. 

92. 

Y220T 

.335 

0.230 

214. 

221 . 

Y300 

.131 

0.079 

40. 

114. 

Y380 

.309 

0.120 

13. 

132. 

Y440 

Y490 

Y517 

.472 

.298 

.514 

0.183 

0.141 

0.074 

9. 

0. 

228. 

162. 

136. 

44. 

XI 40 

.027 

0.019 

315. 

354. 

X180T 

.700 

0.177 

80. 

130. 

X540 

.182 

0.483 

203. 

184. 

X200R 

.038 

0.088 

100. 

52. 

X200L 

.060 

0.057 

261. 

280. 

XI 90 

X220L 

.036 

.018 

0.076 

0.02C 

356. 

314. 

61. 

323. 
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TABLE  64.  ACCELERATION  RESPONSE  FOR  A  LEFT 

ROLLING  PULLOUT  AT  9500  LB  GROSS  WEIGHT 


Structural 

Applied  Force 

imwmm 

Direction 

Real 

Imaginary 

Lateral  at  Tail  Rotor  Gearbox 

-  19. 

-  97. 

Longitudinal  at  hub 

320. 

-985. 

Lateral  at  hub 

-635. 

-572. 

Vertical  at  hub 

729. 

3157. 

Accelerometer 

Identification 

PflEP8 

IWii  imtJJWMHil 

GuPl xS 

ora 

Z50 

.468 

0.632 

126. 

137. 

Z100T 

.329 

0.495 

114. 

120. 

Z210T 

.263 

0.460 

105. 

169. 

Z340 

.167 

0.943 

115. 

58. 

Z400 

.216 

0.686 

71. 

53. 

Z460 

.275 

0.100 

105. 

38. 

Z540 

.615 

1.299 

261. 

233. 

Z90R 

.260 

0.320 

128. 

133. 

Z90L 

.398 

0.480 

92. 

118. 

Z140R 

0.221 

109. 

118. 

Z140L 

.344 

0.408 

75. 

107. 

Z200R 

.299 

0.206 

79. 

38. 

Z200L 

0.895 

59. 

87. 

Z260R 

.420 

0.515 

80. 

82. 

Z260L 

.482 

0.692 

72. 

70. 

Z396R 

.593 

0.450 

64. 

74. 

Z396L 

.191 

0.778 

234. 

37. 

ZLONG 

.187 

0.268 

80. 

82. 

ZLATR 

.273 

0.335 

73. 

85. 

ZCOLL 

0.413 

80. 

93. 

Y50 

.678 

0.422 

269. 

308. 

Y30 

.412 

0.208 

264. 

291. 

Y 1 40 

.273 

0.110 

258. 

269. 

Y220B 

.179 

0.105 

72. 

166. 

Y220T 

1.540 

0.908 

214. 

231. 

Y300 

.584 

0.363 

65. 

154. 

Y380 

.782 

0.732 

53. 

178. 

Y440 

.962 

1.339 

43. 

195. 

Y490 

.312 

1.005 

325. 

332. 

Y517 

3.484 

0.724 

243. 

264. 

XI 40 

0.089 

333. 

60. 

X180T 

2.836 

0.748 

111. 

153. 

.622 

1.311 

262. 

227. 

X200R 

.244 

0.267 

59. 

55. 

•  336 

0.097 

243. 

279. 

XI 90 

WmmM 

0.310 

46. 

92. 

X220L 

■■ 

0.070 

327. 

28. 

158 


TABLE  65.  ACCELERATION  RESPONSE  FOR  A  RIGHT  ROLLING 
PULLOUT  AT  9500  LB  GROSS  WEIGHT 


I  Frequency  3  10.8  Hz 

Structural  Damping 

3  0.122  i 

i  Applied  Force  (lb)  1 

I  Direction 

Real 

Imaginary  ! 

1  Lateral  at  Tail  Rotor  Gearbox 

-  63. 

-  75. 

Longitudinal  at  hub 

287. 

-938. 

Lateral  at  hub 

-392. 

-467. 

Vertical  at  hub 

1532. 

1586. 

Accelerometer 

Magnitude 

(g's) 

Phase  (deg)  i 

Identification 

Analysis 

iTqht 

Analysi 

s 

Flight 

Z50 

.607 

0.651 

142. 

156. 

Z1007 

.319 

0.440 

127. 

136, 

Z210T 

.171 

0.517 

96. 

186. 

Z340 

.174 

0.869 

2. 

40. 

Z400 

.390 

0.624 

354. 

32. 

Z460 

.108 

0.107 

67. 

329. 

Z540 

.391 

1.405 

268. 

218. 

Z90R 

.214 

0.239 

137. 

158. 

Z90L 

.298 

0.440 

105. 

131. 

Z140R 

.068 

0.096 

60. 

123. 

Z140L 

.241 

0.333 

54. 

113. 

Z200R 

.384 

0.307 

5. 

338. 

Z200L 

.700 

0.818 

38. 

85. 

Z260R 

.346 

0.418 

29. 

64. 

Z260L 

.381 

0.641 

35. 

54. 

Z396R 

.466 

0.333 

49. 

49. 

Z396L 

.381 

0.679 

305. 

27. 

ZLONG 

.185 

0.224 

24. 

56. 

ZLATR 

.232 

0.282 

34. 

73. 

ZCOLL 

.321 

0.323 

48. 

81. 

Y50 

.576 

0.454 

266. 

328. 

Y30 

.330 

0.218 

261. 

316. 

Y140 

.215 

0.093 

260. 

303. 

Y220B 

.193 

0.077 

84. 

156. 

Y220T 

.999 
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Figure  29.  Acceleration  response  for  a  45-degree  turn 
at  a  gross  weight  of  9500  lbs. 


Figure  32.  Acceleration  response  for  a  left  rolling 
pullout  at  a  gross  weight  of  9500  lbs. 
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CONCLUSIONS 


1.  Techniques  were  developed  and  implemented  during  dynamic  tests  of  the 
AH-1G  airframe  which  yielded  accurate  acceleration  mobility  frequency 
response  data. 

2.  Analytical  methods  which  were  developed  at  Kaman  Aerospace  Corporation 
were  applied  to  dynamic  test  measurements  of  the  AH-1G  airframe  and  accur¬ 
ately  defined  the  natural  frequencies,  damping  coefficients  and  complex 
orthonormal  mode  shapes  in  the  frequency  range  0-25  Hz. 

3.  Truncated  Model  system  identification  was  performed  successfully  on  the 
dynamic  test  data  of  the  AH- 1 G  airframe.  Truncated  models  can  predict  the 
effect  of  concentrated  mass  changes  and  strut  type  stiffness  changes  on 

the  natural  frequencies  and  orthonormal  modes  of  the  structure. 

4.  The  Incomplete  Model  methodology  was  applied  to  a  92  D.O.F.  NASTRAN 
model  of  the  AH-1G  airframe  resulting  in  an  improved  model  which  predicted 
successfully  the  natural  frequencies  and  frequency  response  of  the  struc¬ 
ture.  Analytical  predictions  of  the  vibration  characteristics  of  the 
airframe  using  the  improved  model  compared  favorably  with  test  data  and 
exhibited  significant  improvement  compared  to  the  predictive  capability  of 
the  original  92  D.O.F.  NASTRAN  model. 

5.  The  main  rotor  and  tail  rotor  forces  obtained  from  Reference  2  were 
applied  to  the  improved  92  D.O.F.  NASTRAN  model  and  the  predicted  acceler¬ 
ations  were  in  agreement  with  the  accelerations  measured  in  flight. 

6.  The  Incomplete  Model  method  proved  to  be  insensitive  to  levels  of  bias 
and  random  error  applied  to  the  measured  orthonormal  modes  used  in  the 
analysis. 

7.  The  method  of  Done  and  Hughes  was  able  to  predict  the  effect  on  natural 
frequencies  and  complex  orthonormal  mode  elements  resulting  from  concen¬ 
trated  mass  changes  and  strut  type  stiffness  changes  made  to  the  airframe. 


RECOMMENDATIONS 


Based  on  the  results  of  this  effort,  it  is  recommended  that: 

1.  The  Incomplete  Model  and  Truncated  Model  methodologies,  which  currently 
operate  on  real  orthonormal  modal  data,  be  extended  to  incorporate  complex 
modes,  thereby  improving  the  vibration  prediction  capability  of  the  methods. 

2.  The  experimental  and  analytical  procedures  that  were  developed  and 
applied  in  the  current  program  be  applied  to  a  composite  structure,  such 
as  a  helicopter  tail  cone,  to  validate  and  promote  acceptance  of  the 
methods  for  all  types  of  structures. 
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APPENDIX  A 

DYNAMIC  TESTING  OF  THE  AH-16  AIRFRAME 


GROUND  VI3RATION  TESTS 


The  dynamic  testing  of  the  AH-1G  airframe  was  conducted  in  Kaman  Aerospace 
Corporation* s  full-size  dynamic  test  facility.  During  shake  tests  the 
airframe  was  suspended  using  a  bungee  system  to  simulate  free  flight.  A 
photograph  of  the  airframe  suspension  system  is  presented  in  Figure  A-l. 

A  full  view  of  the  airframe  in  the  test  facility  is  shown  in  Figure  A-2. 

The  electromagnetic  shaker  system  for  excitation  applied  laterally  at  the 
tail  rotor  gearbox  location  is  presented  in  Figure  A-3.  Force  excitation 
was  also  applied  vertically  at  the  tail  skid  location  and  at  the  hub  loca¬ 
tion  in  the  vertical  direction. 

INSTRUMENTATION 

Thirty  seven  STATHAM  Model  A69TL  strain  gage  accelerometers  were  located 
on  the  airframe  as  described  in  Reference  2.  The  accelerometers  are 
designed  for  flight  and  general  purpose  use  with  a  measurement  range  of 
+5  g's  and  a  natural  frequency  of  approximately  375  Hz.  In  addition, 

55  ENDEVCO  Model  2262-25  piezoresistive  accelerometers  were  placed  on  the 
airframe  during  the  present  study.  The  ENDEVCO  Model  2262-25  accelerometers 
have  a  rated  range  of  +  25g  and  a  natural  frequency  of  approximately  750  Hz. 
The  accelerometers  were  rigidly  attached  to  the  structure,  via  mounting 
blocks,  using  existing  bolt  or  connector  holes  or  with  appropriate  bonding 
agents  to  minimize  extraneous  localized  response  effects.  Typical  acceler¬ 
ometer  installations  are  shown  in  Figures  A-4  and  A-5.  The  accelerometers 
shown  in  the  foreground  on  the  wing  of  the  airframe  are  STATHAM  Model  A69TL 
accelerometers  with  direction  of  major  sensitivity  in  the  vertical  and 
longitudinal  directions.  Also  shown  in  the  figure  are  the  ENDEVCO  Model 
2262-25  with  direction  of  major  sensitivity  in  the  vertical  direction. 

The  accelerometers  shown  in  Figure  A-5  are  ENDEVCO  Model  2262-25.  The 
accelerometer  shown  in  the  lower  part  of  the  figure  is  placed  with  direction 


Figure  A-4.  Typical  accelerometer  installation  on  AH-1G  wing. 


of  major  sensitivity  in  the  lateral  direction,  and  the  accelerometer  in  the 
upper  part  of  the  figure  Is  positioned  with  direction  of  major  sensitivity 
In  the  vertical  direction.  Force  was  measured  using  an  INTERFACE  Load  Cell 
Model  1010-AF. 

The  accelerometers  used  in  the  program  were  wired  to  B&I  Instrument  Model 
24-200  balance  boxes  and  bridge  power  supplies.  The  balance  boxes,  used 
as  signal  conditioning  units,  had  their  output  signals  directed  to  a 
Pulse  Code  Modulation  encoder  (PCM).  Data  signals,  at  excitation  frequen¬ 
cies  of  interest,  were  transmitted  from  the  test  site,  stored  on  magnetic 
tape,  and  processed  on  the  IBM  360/50  to  yield  the  acceleration  mobilities 
of  the  airframe.  Prior  to  data  acquisition  and  reduction  selective  accelero¬ 
meter  and  force  data  were  processed  using  the  HP5420A  Digital  Signal  Analy¬ 
zer.  Both  random  and  sinusoidal  force  excitation  were  used  in  this  phase 
of  the  program  to  check  reciprocity  and  linearity  of  the  mobility  accelera¬ 
tion  frequency  response  of  the  airframe.  Figure  A-6  presents  instrumenta¬ 
tion  used  in  the  dynamic  testing  of  the  AH-1G  airframe. 


TECHNIQUES  AND  PROCEDURES  FOR  VIBRATION 
TESTING  OF  THE  AH-1G  HELICOPTER _ 


The  development  of  effective  techniques  for  shaking  the  AH-1G  and  analyz¬ 
ing  the  acquired  vibration  data  constituted  a  major  part  of  the  research 
work  associated  with  testing.  The  validity  of  the  methods  employed  rests 
heavily  upon  the  consistency  between  the  measured  structural  mobilities 
and  the  theoretical  models  for  which  these  mobilities  are  derived.  This 
consistency  is  critical  since  the  measured  mobilities  are  used  not  only 
for  obtaining  the  global  modal  parameters  of  the  test  vehicle,  but  also 
to  derive  mobilities  which  were  not  measured  directly. 

Digital  signal  analyzers  have  made  it  possible  to  measure  the  response  of 
structures  to  any  physically  realizable  excitation.  However,  the  inter¬ 
pretation  of  measured  structure  to  specified  excitation  forces  is  subject 
to  the  mathematical  model  used  in  the  process  of  analyzing  the  data.  The 
model  may  be  more  or  less  sophisticated,  depending  on  the  test  data.  The 
dynamic  testing  of  a  structure  like  the  helicopter  poses  a  number  of  speci¬ 
fic  problems.  These  problems  are  associated  with:  (1)  the  size  and  com¬ 
plexity  of  the  structure;  (2)  nonuniform  distribution  of  mass,  stiffness 
and  damping;  and  (3)  the  correct  application  of  linear  vibration  theory 
to  the  process  of  data  acquisition  and  analysis. 

It  has  been  implied  that  the  techniques  adopted  for  the  structural  dynamic 
testing  are  closely  related  with  the  theory  underlying  the  vibration  analy¬ 
sis.  Therefore,  a  discussion  of  the  specific  test  procedures  must  necess¬ 
arily  be  preceded  by  a  brief  summary  of  the  theoretical  considerations  in¬ 
cluding:  (1)  the  theory  of  the  generalized  linear  structure;  (2)  the  prin¬ 
cipal  characteristics  of  acceleration  mobility  data^  (3)  testing  procedures 
for  global  parameters  and  the  estimation  of  these  parameters;  (4)  testing 
procedures  for  obtaining  mo(  shapes  and  the  method  of  calculating  mode 
shapes;  and  (5)  methods  for  -ing  mobilities  from  modal  data. 
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THEORY  OF  THE  GENERALIZED  LINEAR  STRUCTURE 


The  dynamic  properties  of  any  structure  can  always  be  characterized  by  a 
relationship  between  a  selected  set  of  motion  coordinates  and  the  set  of 
externally  applied  forces,  i.e.. 


Motion' 

Character  of 

'Force 

Vector 

Structure 

X 

Vector 

The  character  of  the  structure  implied  in  equation  (A-l)  will  be  termed 
mobility.  '  If  the  motion  vector  is  a  vector  of  displacements/velocities/ 
accelerations,  the  character  of  the  structure  is  termed  displacement/ 
veloity/acceleration  mobility,  rospecti vely. 

The  central  phenomenon  of  vibration  theory  is  cyclic  oscillation,  hence 
the  quantities  that  go  into  equation  (A-l)  are  generally  sought  in  the 
frequency-domain;  for  example,  in  acceleration  measurements 

{y(«)>  =  [Y(a»)]{f  («)}  (A-2) 

where  {y(u)}  is  the  Fourier  transform  of  the  accelerations;  [Y(«)]  is  the 
acceleration  mobility  matrix;  and  {f (oj) }  is  the  Fourier  transform  of  the 
vector  of  generalized  forces,  compatible  with  the  selected  set  of  coor¬ 
dinates.' 

From  a  measurement  standpoint,  the  jkth  element  of  the  matrix  [YU)] 
relates  the  acceleration  measured  along  the  jth  coordinate  when  the  only 
force  acting  on  the  structure  is  that  applied  along  the  kth  coordinate; 
i.e., 

iy  ')  =  YjkU)  fk(o>)  When  fi  ^  k  =0  (A-3) 

Linear  vibration  response  of  a  structure  may  be  characterized  by  the  fol¬ 
lowing  conditions:  (1)  the  response  of  the  structure  to  random  forcing 
is  stationary  in  time  (i.e.,  forced  vibrations  are  steady);  (2)  the  ele¬ 
ments  of  the  matrix  [Y(w)]  are  functions  of  frequency  only,  and  depend  on 


178 


neither  the  motion  coordinates  nor  the  forcing  vector;  and  (3)  the 
mobility  matrix  [Y(co)]  is  symmetric;  i.e.,  =  Y^ . 

The  foregoing  conditions  have  specific  practical  implications  in  vibration 
testing  and  analysis.  The  first  condition  is  necessary  for  any  structure 
to  survive  continuous  operation  under  arbitrary  dynamic  excitation.  The 
second  condition  more  or  less  stipulates  the  type  of  shake  test  data  that 
is  adequate  for  analysis  based  on  a  linear  model  of  the  structure.  If 
the  mobility  functions  measured  for  different  force  levels  are  not  the 
same,  the  assumption  of  linearity  is  not  satisfied.  This  is  usually  the 
case  when  only  part  of  the  structure  may  be  participating  in  the  response. 

As  the  force  level  is  increased,  more  and  more  of  the  relevant  motion  coor¬ 
dinates  of  the  structure  start  to  participate  in  the  response.  The  range 
of  linear  response  is  reached  only  when  the  measured  mobility  remains  un¬ 
changed  with  changing  force  levels.  The  third  requirement  is  that  of 
reciprocity.  If  the  shaking  and  measurement  stations  are  interchanged, 
the  same  mobility  should  be  recorded,  otherwise  the  [Y]  matrix  will  not 
be  symmetric,  as  required  by  the  linear  model. 

In  the  foregoing  characterization  of  a  linear  system,  no  assumptions  are 
made  about  the  nature  of  the  damping  mechanisms  occurring  in  the  structures. 
All  the  conditions  required  for  linear  modeling  can  be  verified  in  the 
process  of  the  actual  shake  test  of  the  structure. 


The  relationship  between  the  Fourier  transform  of  the  force  vector  and 
that  of  the  displacement  vector  of  a  steadily  vibrating  undamped  multi¬ 
ple  degree  of  freedom  system  can  be  written  as: 
r 


2 

O) 


[M]  +  [K] 


{y(«)j  = 


f  (to) 


(A-4) 


where  [M]  and  [K]  are  real,  symmetric  mass  and  stiffness  matrices,  respec¬ 
tively.  Thus,  the  displacement  mobility  matrix  Tor  an  undamped  system  is 
simply 


[YU)] 


U 


a)  [M]  +  [K]J 


-1 


(A-5) 
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The  presence  of  damping  in  its  most  general  form  can  be  modeled  by  intro¬ 
ducing  a  frequency-dependent  complex  damping  matrix  into  equation  (A-5), 
i  .e. , 


[Y(«)]„ 


u2  [M]  +  [K]  +  [DR(u)]  +  i [D1 (w)] 


(A-6) 


This  analytical  development  has  meaning  only  in  the  frequency  domain  for 
the  general  case  of  damping.  This  is  mainly  because  the  physical  quanti¬ 
ties  that  can  be  used  to  characterize  the  arbitrary  damping  of  a  structure 
are  related  to  the  energy  dissipated  per  cycle  of  oscillation.  In  cases 
where  the  time  domain,  force/motion  relationship,  representing  the  damping 
mechanism  is  known,  the  damped  equations  of  motion  can  be  developed  in  the 
time  domain  and  then  Fourier  transformed  into  the  frequency-domain. 
However,  as  described  in  Reference  9,  the  inverse  Fourier  transform  of  the 
frequency-domain  equations  that  may  adequately  describe  an  arbitrarily 

I 

damped  system  may  not  yield  a  time-domain  system  of  equations  that  makes 
physical  sense.  In  other  words,  arbitrary  damping  mechanisms  may  not  be 
susceptible  to  a  time-domain  description.  Mathematical  models,  developed 
from  time-domain  equations  of  motion,  usually  fail  to  identify  global 
characteristics  of  structures  with  significant  damping. 


D 

In  general,  the  elements  of  [DU)]  are  small,  compared  to  those  of  the 
[K]  matrix.  Also,  in  order  for  reciprocity  conditions  to  be  met  and 
for  energy  to  be  dissipated,  the  damping  matrix  must  be  symmetric  and 
nonnegative  definite  over  the  entire  frequency  range. 


For  a  damped  system,  then 


|[K]  +  i [DU)]  -  a) 2 


{yUb  =  {fUb 


(A-7) 


9.  Crandall,  S.  H.,  The  Roll  of  Damping  in  Vibration  Theory.  Journal 
of  Sound  and  Vibration,  Vol .  11,  No.  1,  1970. 
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Consider  the  complex  frequency-dependent  characteristic  value  problem 


f[K]  +  1[0(»)]1 


{()>}  =  A(uj)[M]  { } 


(A-8) 


R  I 

where  (4>}  =  (4>  }  +  i {<J>  }  is  the  complex  characteristic  vector  which  can  be 

R  I 

assumed  to  be  frequency  independent;  A(w)  =  A  ( w)  +  i A  (w;  is  the  frequen¬ 
cy  dependent  complex  eigenvalue. 


If  combinations  of 
equation  (A-8),  then 


Aj((jo), 


and 


( 


\(w)»  (4>) 


k 


exist,  which  satisfy 


{4.} 


T 

k 


[K]  +  i [D(w )] 


(4>}j  =  Aj((ij){4>}J 


(A-9) 


and 


(4>}T 


[K]  +  1[D(u)] 


IKK^k 


(A- 1 0 ) 


{g>} T  denotes  the  transpose  of  (4>).  By  virtue  of  the  symmetry  of  the  [K], 
[M]  and  [D(co)]  matrices,  equations  (A-9)  and  (A-10)  lead  to  the  follow¬ 
ing  orthogonality  relationships: 


and 


where 


<+>J  [M]  {*>k  -  mj 


[K]  +  i  [D(co)  ] 


<*>k  = 


k.  * 


5jk 


(A-ll) 

(A— 1 2 ) 


nij  =  (4>}J  [M]  {cp}^ 
kj  =  w]  [K]  (4>}j 
dj(«)  =  (4»}]  [0(a))].{4>}j 


(A- 13) 
(A— 1 4 ) 
(A-15) 
(A-l 6 ) 
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It  follows  that 


w] 


[K]  -  w  [M]  +  i[D(to)] 


<«k  = 


kj  -  to  idj(w) 


"jk 


A.(u)  -  of  m  . 


6.. 

Jk 


(A-17) 


If  the  vectors  {<]>}.  exist,  it  can  easily  be  verified  that  only  the 

*J 

imaginary  parts  of  \Aw)  need  be  frequency  dependent.  So  that 

J 


Aj(w)  =  Aj  +  i  Xj  (to) 


(A-18) 


By  post-multiplying  the  transpose  of  equation  (A-8)  by  which  is 

the  complex  conjugate  of  (<R,  the  following  equation  is  obtained: 

J 


{«I 


[K]  -  i [D(w)] 


T 


{*}j  =  Xj(w) {<«,}]  [M]^ 


( A  - 1 9 ) 


Similarly,  the  complex  conjugate  of  equation  (A-8)  can  be  premultiplied 
by  {oil  to  get 


{*}j 


[K]  -  iCO(co)] 


{<M  .•  =  A  (U))(d>}  ‘  [M]{(j>}  . 

J  v*  J  J 


(A-20) 


From  equations  (A-19)  and  (A-20) 


A.  (to)  +  A* (to)  =  2  Wl  [K]  [M]  4)* 

J  J  J  J  J  J 


(A-21) 


and 


A  (to)  -  A* (to)  =  2iW7  [[)(<*>)]  (❖)*/ (<5)}T[M]{<)))* 
J  J  J  j  J  J 


(A-22) 


The  right-hand  side  of  equation  (A-21)  is  a  frequency  independent  quantity. 
However,  the  right-hand  side  of  equation  (A-22)  is  frequency  dependent, 
establishing  the  validity  of  the  claim  made  in  equation  (A-18). 
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A  complex  L  x  N  modal  matrix  [4>]  can  be  defined  such  that  its  jth  column 
is  the  L  x  1  vector  {<}>}  .;  j  =  1,  2,  .  ..N,  where  L  is  the  number  of  coor- 

J 

dinates  chosen  to  describe  the  system  and  N  is  the  number  of  modes  of 
the  system.  In  principle,  N  is  infinite;  however,  in  practice,  over  a 
given  frequency  range,  only  a  finite  number  of  system  modes  are  necessary. 


Equation  (A-7)  can  be  rewritten  to  give 


{ y  (oo ) }  =  [$] 


[<D]Tf[K]  -  u2  [M]  +  i[DM]l 


[*] 


-1 


[*]V(u>)}  (A-23) 


and,  using  the  orthogonality  relationships,  leads  to  the  results 


{:*(«)}  -  01 


1 

<>?  -  “2)  +  1>!(W 

l  . 

[*]'{f(u)} 


(A-24) 


By  definition,  {  y (co )}  =  [Y(oj)]{ f (to)} ;  hence, 


N 

CY(«)3  =  I 

n=l 


{4,}  {<J>}' 

n  n 


m„ 


1 


(AR  -  w2)  +  i  A*  (to) 
x  n  '  nv  ' 


(A-25) 


R  I  2 

An  and  Xp (to)  have  units  of  (frequency)  ,  and  from  pnvsical  considerations, 

both  a||  and  An(co)  are  positive.  Therefore,  it  is  possible  to  define 


An  -  “n 

n  n 


(A-26) 


and 


aJ(u)  =  gn(w)«n 


(A-27) 


The  matrix  of  modal  acceleration  coefficients  of  the  nth  mode  is  defined 
as 


[A] 


n 


_  1_ 

mn 


(A-28) 


•-  -  -  .J>.  rV  V.  •*. 
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The  acceleration  mobility  matrix  and  the  displacement  mobility  matrix  are 
related  by 


[¥(«■)]  =  -  »2[yM] 


(A-29) 


Using  equations  (A-26),  (A-27),  (A-28),  and  (A-29),  the  jkth  acceler¬ 
ation  mobility  can  be  written  as 


N 


(w)  =  •  l  Aj kn 
n=l 


w2/ft  2 

_  '  n _ 

(1  -  oZ/ft^)  +  ignH 


R  I 

and  Ajkn  =  Ajkn  +  1  Ajkn’  t*ie  comPlex  modal  acceleration  for 

response  at  coordinate  j  and  force  excitation  k  for  mode  n. 


(A-30) 


ftn  =  the  natural  frequency  of  the  nth  mode 
gn  =  damping  coefficient  of  the  nth  mode 

to  -  frequency  of  excitation 


In  the  most  general  case,  the  dependence  of  gn(<o)  on  frequency  may  not  be 
known.  However,  it  is  expedient  to  take  advantage  of  the  fact  that  the 
ig  (to)  term  in  equation  (A-30)  is  dominant  only  in  the  frequency  range 

ii  2  2 

near  the  natural  frequency  of  the  nth  mode  where  w  /ftn  =*  1.  Thus,  any 
suitable  representation  of  gn(to)  which  matches  the  correct  value  in 
the  neighborhood  of  u  =  fin.  may  be  assumed. 


The  general  form  of  the  jkth  element  of  the  acceleration  mobility  matrix 
can  be  written  as 


N 

Y„.  -  ?J„  ♦  1yL  ■  E$k  ♦  IEjk  -  l  Aj 

n=l 


to2 /ft  2 

n 


jk  jk 


‘  jk 


kn  (1  -  to2/ft2n)  +  ign(03) 


(A-31) 


R  I 

where  E^  +  iE^  represents  the  rigid  body  acceleration  coefficients.  In 

most  cases,  is  very  smell  compared  to  the  rest  of  the  terms  in  the 
series  and  is  usually  neglected. 
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(A-32) 

(A-33) 

(A-34) 

(A-35 ) 

(A-36) 

( A-37 ) 


Equations  (A-32),  (A-33),  and  (A-34)  represent  the  measured  acceleration 
mobility  as  a  linear  combination  of  the  mode  functions. 
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Plots  of  F  (u)  and  F  (u)  as  functions  of  frequency  ratio  for  three 
values  of  the  damping  coefficient  are  shown  in  Figure  A-7.  The  Nyquist 
plots  of  the  complex  F(w)  functions  are  shown  in  Figure  A-8. 

,,R 

The  F  (w)  function  is  characterized  by  two  peaks  at 


“in  *  “n/1*  9n2  '  \J  1  +  <A-38> 

and 

“2n  fin/  +  9„2  +  9n  ]\  +  9„*  (A-39) 


while  the  F*(w)  function  has  only  one  peak  at  /  l  +  gp2 

Note  that 


w2n 


2 


(A-40) 


which  increases  with  increased  damping. 

"R 

From  the  plots  in  Figure  A-7,  it  is  seen  that  linear  combinations  of  F 
and  F*  vary  rapidly  in  the  vicinity  of  the  natural  frequency,  and  are 
either  negligible  or  slowly  varying  with  frequency  in  the  regions  away 
from  the  natural  frequency. 
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Separated  Modes  -  Equation  (A-31 )  carries  the  basic  implication  that  the 

effects  of  the  structure's  modes  occurring  at  different  frequencies  on 
the  measured  mobility  are  additive  in  the  frequency  domain.  If  a  mode 
occurs  at  a  frequency,  in  the  neighborhood  of  which  the  contributions  from 
the  other  modes  of  the  structure  are  either  negligible  or  are  weakly  vary¬ 
ing  with  frequency,  such  a  mode  is  said  to  be  well  separated.  The  nature 
of  the  measured  mobility  in  this  frequency  range  will  be  dominated  by  that 
particular  mode. 


Classical  Modes  -  In  the  case  of  a  classical  mode,  i.e.,  when  the  system 
mode  shape  is  the  same  for  the  damped  system  as  it  would  be  for  the  un¬ 
damped  system  the  A^n  is  a  real  number  ^A^  =  0  j  »  and  real 

part  of  the  measured  acceleration  mobility  will  show  two  turning  points 
for  each  separated  mode  and  the  imaginary  part  will  show  a  single  turning 
point  only.  For  a  classical  mode,  equation  (A-40)  can  be  approximated 
to  give  an  estimate  of  the  damping  coefficient:  gn  -  (co2n  -  wln)/^n  . 

Complex  Modes  -  For  the  general  case  of  nonclassical  or  complex  modes, 

R  I 

both  Ajj^  and  Ajkn  are  significant.  The  measured  real  and  imaginary  mobili- 

**  R 

ties  of  a  well-separated  mode  contain  linear  combinations  of  both  F  (w)  and 

F1 (oi)  in  proportions  given  in  equations  (A-32)  and  (A-33).  In  particular, 

if  Ajkn  »  Ajkn,  it  is  the  imaginary  part  of  the  acceleration  mobility 

which  will  show  two  turning  points,  while  the  real  part  will  show  a  single 
turning  point  only. 
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Coupled  Modes  -  System  modes,  occurring  In  frequency  ranges  such  that 
their  mutual  contributions  to  the  measured  mobility  in  this  frequency 
range  are  rapidly  varying  functions  of  frequency,  are  said  to  be  coupled. 

Mode  Clusters  -  A  mode  cluster  (or  a  cluster  of  modes)  is  characterized 
by  a  group  of  system  modes  which  are  coupled  together  by  virtue  of  the 
proximity  of  their  resonances.  Mode  clusters  are  usually  separated  by 
regions  of  negligible  or  slowly  varying  mobility  values  in  the  frequency 
domain. 

Mode  clusters  generally  have  the  appearance  of  single  modes  in  wideband, 
low  frequency  resolution  mobility  measurements.  Higher  resolution  data 
usually  helps  to  reveal  the  modal  content  of  a  particular  mode  cluster. 
Figure  A-9  shows  broad  band  (0  -  200  Hz)  mobility  of  a  helicopter  vertical 
tail  shake,  measuring  vertical  acceleration  at  the  nose.  Between  mode 

clusters,  measured  mobility  varies  slowly  close  to  the  zero  value. 

In  fact,  what  appears  to  be  a  single  mode  in  the  0  -  10  Hz  frequency  range 
is  actually  a  cluster  of  two  modes  as  shown  in  Figures  4a  and  4b  ,  which 
are  higher  resolution  measurement  of  the  same  mobility  in  the  5  -  10  Hz 
range. 

Identification  of  mode  clusters  is  useful  in  determining  which  modes  should 
be  included  when  truncating  equation  (A-31 )  and  in  identifying  frequency 
bandwidths  for  higher  resolution  data  acquisition. 
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Figure  A-9.  Measured  acceleration  mobility  of 
the  AH-1G  between  0  Hz  and  200  Hz. 
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SHAKE  TESTING  FOR  GLOBAL  PARAMETERS 


Each  elastic  mode  of  a  structure  is  characteris'd  by  a  natural  frequency, 
n  ,  and  a  damping  coefficient,  gn,  w‘  <ch  are  global  properties  of  the  struc¬ 
ture.  These  are  the  only  constants  that  enter  into  the  mode  frequency 
functions.  They  are  the  same  for  a  given  mode,  regardless  of  the  response 
coordinate.  The  first  stage  of  modal  testing  is  to  determine  the  global 
parameters  of  the  dominant  elastic  modes  which  occur  in  the  frequency 
range  of  interest. 

The  experimental  data  required  for  determining  the  global  parameters  are 
the  continuous  frequency  plots  of  a  number  of  mobilities  which  are  con¬ 
sidered  to  represent  the  global  vibrational  behavior  of  the  structure. 

For  a  selected  set  of  shaking  locations,  e.g.,  tail  vertical  and  tail  rotor 
gearbox  lateral,  the  transfer  functions  between  the  response  coordinate  and 
the  shaking  coordinate  are  measured  over  the  determined  frequency  range. 
Typical  response  coordinates  for  such  measurements  are:  (1)  nose  vertical, 

(2)  wing  (right  and/or  left)  vertical,  (3)  center-of-gravity  vertical, 

(4)  tail  vertical,  and  (5)  horizontal  stabilizer  vertical. 

The  test  setup  for  measuring  frequency-dependent  mobility  functions  is 
shown  in  Figure  A- 10.  The  helicopter  is  suspended  as  a  free  body  by  soft 
rubber  bungee  chords.  The  configuration  shown  has  the  shaker  located  for 
vertical  excitation  at  the  tail  and  the  accelerometer  located  for  measuring 
vertical  response  in  the  vicinity  of  the  horizontal  stabilizer.  Signals 
for  driving  the  electromagnetic  shaker  originate  from  the  signal  generator. 

A  force  gage,  installed  at  the  point  of  force  application,  generates  a 
voltage  signal  which  is  proportional  to  the  applied  force.  This  signal 
is  input  to  the  dual  channel  digital  signal  analyzer.  The  accelerometer 
at  the  response  coordinate  generates  a  voltage  signal,  proportional  to 
the  response  acceleration,  which  is  also  input  to  the  digital  signal 
analyzer. 

The  signal  analyzer  is  capable  of  sampling  the  time-domain  force  and 
response  signals,  digitizing  these  samples,  and  computing  the  real-time 
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Fourier  transforms  of  the  data.  It  also  computes  the  frequency-domain 
transfer  function  between  the  input  and  output  spectra,  which  is  the 
mobility  between  the  response  and  forcing  coordinates.  All  the  frequency 
functions  computed  by  the  analyzer  over  the  specified  frequency  interval 
can  be  stored  on  cassette  tapes  for  future  restoration  and  analysis.  The 
oscilloscope  allows  the  monitoring  of  the  time-domain  signals  emanating 
from  the  force  and  response  transducers.  The  frequency  counter  is  used 
to  precisely  measure  the  frequencies  of  harmonic  signals  when  required. 

The  accuracy  with  which  the  global  parameters  car,  be  estimated  is  criti¬ 
cally  dependent  on  the  quality  of  the  data  acquired  for  this  purpose. 

For  each  pair  of  force  and  response  locations,  a  random  shake  is  done  with 
the  frequency  bandwidth  set  to  span  twice  the  frequency  range  of  interest. 
This  is  done  to  insure  that  the  modes  in  the  frequency  range  of  interest 
are  not  coupled  to  modes  occurring  beyond  the  range,  as  may  be  the  case 
when  a  local  mode  is  present.  For  each  shaking  station,  several  force 
levels  are  tested  until  the  range  of  applied  force  is  reached  where  the 
mobility  plots  are  independent  of  the  force  level.  This  is  one  of  the 
linearity  requirements  for  mobility  measurement.  Having  established  the 
required  force  level  and  the  absence  of  local  coupling  modes  at  higher 
frequency,  another  random  shake  is  done,  this  time  with  the  bandwidth  set 
at  the  frequency  range  of  interest.  The  above  procedure  is  repeated  for 
all  the  accelerometers  which  have  been  selected  for  global  parameter 
testing. 

The  ratio  of  modal  acceleration  coefficient  to  damping  tAjkn/9n)  varies 
not  only  from  mode  to  mode  but  also  from  mobility  to  mobility,  for  a 
given  mode,  and  the  prominence  of  the  various  modes  of  the  structure  will 
be  different  in  each  of  the  mobilities  recorded.  That  is,  a  given  mode 
n  occurring  at  may  appear  very  prominently  on  mobility  Y^D,  while 
the  same  mode  may  not  be  so  significant  in  the  mobility  Yj&,  where  z 
designates  a  response  coordinate  different  from  k.  This  will  especially 
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be  the  case  If  the  mode  shape  associated  with  mode  n  has  a  much  larger 
mode  element  at  coordinate  k  than  coordinate  a.  The  prominence  of  a  mode 
may  also  be  due  to  light  damping.  Thus,  by  examining  the  set  of  broad¬ 
band  mobilities  recorded,  It  is  possible  to  associate  each  mode  n  with 
the  mobility  where  the  mode  appears  most  prominently. 

Although  it  is  possible  to  obtain  rough  estimates  of  the  natural  frequen¬ 
cies  and  damping  of  the  structural  elastic  modes  from  these  broad-band 
mobility  plots,  especially  when  damping  is  very  light  (e.g.,  peaks  of  the 
imaginary  mobility  plots,  and  frequency  separation  of  the  peaks  in  the 
real  mobility  plots),  there  are  a  number  of  specific  considerations  why 
broad-band  mobility  data  is  not  suitable  for  global  parameter  extraction. 
Among  these  considerations  are: 

Measurement  Accuracy  -  The  low-frequency  resolution  associated  with 
broad- band  mobility  measurements  tends  to  introduce  errors  into  the  measured 
mobility  values  due  to  the  phenomenon  of  leakage  —  a  spreading  of  ihe 
energy  contained  at  each  discrete  frequency  over  a  relatively  narrow 
band  nearby.  Although  considerable  effort  Is  exerted  by  equipment 
manufacturers  in  reducing  leakage  effects  by  using  window  techniques, 
the  phenomenon  still  has  to  be  reckoned  with  when  the  frequency  resolution 
gets  below  certain  limits.  For  acceptable  measurements,  bandwidths  of 
about  25%  of  the  center  frequency  have  been  recommended. 

Parameter  Extraction  Accuracy  -  Also  associated  with  low-frequency  resolu¬ 
tion  are  inaccuracies  in  the  parameter  extraction  methods  due  to  the  fre¬ 
quency  spacing  between  successive  data  points.  The  polar  plot  of  mobili¬ 
ties  (see  Figure  A-8),  in  the  vicinity  of  a  mode,  describes  a  circular  arc. 
Most  methods  for  extracting  natural  frequencies,  damping  and  modal  accelera¬ 
tion  coefficients  are  based  on  fitting  a  continuous  circle  through  measured 
data  and  in  some  cases  computing  the  rate  of  change  of  the  arc  length  with 
frequency.  Since  the  frequency  data  is  discrete,  arcs  of  the  circle  are 
necessarily  approximated  by  segments.  The  error  incurred  by  approximating 


a  circular  arc  by  a  straight  line  segment  Increases  as  the  frequency  spac¬ 
ing  between  successive  data  points  Increases.  Narrow  band  data  with  band¬ 
width  less  than  25%  of  the  natural  frequency  of  a  given  mode  have  been 
found  to  yield  sufficiently  accurate  results.  Initial  estimates  of  the 
natural  frequencies  can  be  obtained  from  the  broad-band  data. 

For  sufficient  frequency  resolution  and  to  minimize  leakage,  the  follow¬ 
ing  bandwldths  are  recommended  for  use  in  narrow  band  testing  using  the 
HP5420A  signal  analyzer. 


TABLE  A-l .  BANDWIDTH  RECOMMENDATIONS. 


Natural  Frequency  n 

n 

Equal  to  or  greater  than  But  less  than 

Bandwidth 

2  Hertz 

3  Hertz 

.5  Hertz 

3 

4 

.781 

4 

6 

1.000 

6 

8 

1.5625 

8 

12 

2.000 

12 

16 

3.125 

16 

25 

4.000 

25 

32 

6.250 

32 

50 

8.000 

50 

64 

12.500 

64 

100 

16.000 

In  cases  of  gn  greater  than  ,25  broader  bandwidth  should  be  used  and  in  all 
cases  the  natural  frequency  should  be  the  center  frequency. 


SWEPT  SINE  TESTING 


For  all  the  narrow  band  mobility  measurements,  the  excitation  was  achieved 
by  applying  pure  sine  wave  signals  to  the  electromagnetic  shaker  and  vary¬ 
ing  the  frequency  of  the  sine  waves  over  the  range  spanned  by  the  band¬ 
width.  This  so-called  swept  sine  technique  was  preferred  to  other  excita¬ 
tion  techniques  over  a  narrow  frequency  band.  Other  reasons  for  choosing 
the  swept  sine  technique  are  that: 

1.  The  energy  input  Into  each  measurement  frequency  is  maximum. 

2.  By  choosing  the  right  sweep  speed  (see  sketch  below),  the  steady- 
state  sinusoidal  response  of  the  structure  is  achieved  at  each 
measurement  frequency.  This  is  one  of  the  assumptions  made  in 
the  derivation  of  the  generalized  linear  model. 

3.  Measurements  are  more  accurate  and  reproducible. 

4.  The  sampling  frequencies  and  the  adequate  number  of  averages 
are  more  easily  determined. 

5.  Good  linearity  and  reciprocity  checks  are  obtained. 

6.  High  resolution  of  close  modes  can  be  achieved  by  selecting 
the  right  sweep  speed. 

Consider  an  undamped  single  degree  of  freedom  linear  system,  described 
by  the  following  equation  of  forced  vibrations: 

mx  +  kx  =  Feiu)t  (A-41 ) 

If  the  forcing  frequency  coincides  with  the  undamped  natural  frequency, 
i.e.,  w  =/j|  ,  the  response  of  the  system  grows  linearly  with  time. 
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Schematically: 

F(t) 


System 


Input 


Response 


The  undamped  response  Is  governed  by 


••  X  n2  F  jnt 
x  +  «  x  =  -  e 
m 


or 


*  <*>  *  -1 


(A-43) 


However,  because  of  the  various  dissipative  mechanisms  which  constitute 
damping,  the  oscillations  reach  a  limiting  amplitude  after  some  time  t. 


For  hysteretic  damping,  g,  the  steady-state  response  is  governed  by 


x  +  (1  +  1g)n2x  =  jjj-  e1nt 


(A-44) 
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or  x(t) 


(A-45) 


The  steady -state  amplitude  Is  given  by 

Xs  -  (A-46) 

mgn 


and  the  characteristic  time  for  reaching  steady-state  response  can  be 
estimated  by  equating 


Xt(t)  «  Xs 

(A-47) 

F  F 

which  gives  *-pr  r  -  — *■ 

2mfl  mgn2 

(A-48) 

2  1 

thus  T  *  —  =  - — 

(A-49) 

ga  wgf 


where  f  =  frequency  In  Hz. 

Suppose  there  are  two  neighboring  structural  modes  with  the  natural  fre¬ 
quencies  separated  by  Af,  H2.  To  resolve  these  two  close  modes,  the  speed 
at  which  the  excitation  frequency  is  changing  must  be  of  the  order  of 

Af 

v  =  —  =  Trgf  Aj.  Hz/sec  (A-50) 

where  g  *  damping  coefficient  (lower  bound) 

A^  -  mode  resolution  in  Hz. 

or 

v«  f  (A-51) 

for  a  given  frequency  separation  between  modes  to  be  determined  and  for  a 
given  damping  coefficient. 
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The  relationship  between  the  linear  scale  and  the  logarithmic  scale  on 
the  signal  generator  is 


(A-52) 


where  ■  frequency  In  hertz  (fQ  <_  fHz  <.  lOfo) 
fdec  "  frequency  In  decades  (0  <  fdec  <  1.0) 
f  *  base  frequency  on  the  scale 


From  equation 

(A-52)  f„2 

fo  x 

10fdec 

(A-53) 

Sweep  speed  v 

.  dfHz 
dfdec 

dfH2 

df.er 

x  Hz/sec. 

(A-54) 

From  equation 

(A-53)  » 

dTdec 

(f0  ln 

10)  x  10fdec  =  fHz  An  10 

(A-55) 

and  r—  * 

tHz 

dfdec 

dt 

JlnlO 

(A-56) 

Thus^by  selecting  a  constant  logarithmic  sweep  speed  (dfdec/dt  r  const.  *  a), 
equation  (A-51)  is  automatically  satisfied. 

The  constanta  is  determined  by  substituting  the  desired  value  of  -J- 

THz 

into  equation  (A-56).  For  example,  if  at  2  Hz  we  desire  a  sweep  rate  of 
0.01  Hz/sec,  then 

a  s  dec/mi n  =  .13  dec/min  (A-57) 
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ESTIMATION  OF  GLOBAL  PARAMETERS 


Various  techniques  have  been  developed  for  estimating  the  natural  frequen¬ 
cies  and  damping  coefficients  of  the  elastic  modes  of  a  structure  from 
mobility  data.  In  all  cases,  certain  assumptions  have  to  be  made  about 
these  modes.  The  simplest  case  Is  when  the  mode  Is  well  separated  and 
lightly  damped.  For  such  modes,  the  natural  frequency  can  be  approximated 
by  the  peak  of  the  Imaginary  displacement  mobility.  The  damping  coeffi¬ 
cient  can  be  estimated  as 


* 


1 

7 


2  2 
u2  - 


T 


Wg  ■ 

Q 


(58) 


where  ug  and  are  the  turning  point  frequencies  In  the  real  displace¬ 
ment  mobility.  The  above  simple  case  is  almost  exclusively  reserved  for 
simple  structures  with  uniform  distribution  of  mass,  stiffness,  and  damp 
ing.  Very  few  of  the  modes  of  the  helicopter  can  be  treated  this  way. 


2 

The  ds/df  method  of  Kennedy  and  Pancu  -  The  following  is  a  more  general 
approach  which  has  been  found  to  work  well  for  both  classical  and  complex, 
close  or  separated  modes.  By  analogy  with  equations  (A-32)  through  (A-37), 
the  jkth  displacement  mobility  can  be  expressed  as 


jk 


Ejk  ? 

a  "2  +  i 

l 


j'kn  rn  Mjkn  hn 


+  i 


Mjkn  hn  Mjkn  hn 


(A-59) 
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where  Ajkn  *  Ajkn  +  iAjkn  is  the  jkth  modal  acceleration  coefficient 
of  nth  mode  and  Ejk  is  the  contribution  from  the  rigid  body  modes. 

Recall  that 

F>)  -  f[Jm  +  i  f’(u)  =  F  (u)  (A-60) 

0) 
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which  gives 


In  the  immediate  vicinity  of  the  nth  natural  frequency,  the  displacement 
mobility  can  be  approximated  by 

Yjk  W  .  k\n  F„R  U  -  Ajkn  f'  U  -C* 

“n 

+  1  k\n  Fn<“>  +  Ajkn  Fn«  +  Ci  4  +  dil  (A-63) 

nn 

In  equation  (A-63),  the  sum  of  the  contributions  from  all  other  modes  has 
been  represented  by  a  complex  straight  line 

(c* + 1  ci]  £  +  < + 1  < 

Dropping  the  subscripts  j,  k,  and  n,  and  writing  the  real  and  imaginary  parts 
of  the  displacement  mobility  separately,  gives 


YR  (oj*n)  *  AR  FR  (w)  -  A1  F1  (w)  +  CR  ~  +  dR  (A-64) 

n2 

Y*  (w-n)  -  AR  F1  (w)  +  a1  FR  (to)  +  c1  ^  +  d1  (A-65) 

n 
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If  the  nth  mode  Is  classical  and  well  separated  the  Imaginary  part  of  the 
modal  acceleration  coefficient  vanishes  and  the  contributions  from  other 
modes  are  nearly  Independent  of  frequency.  In  other  words,  A*,  CR,  and 
C1  vanish.  Thus, 

YR  Ml)  «  AR  FR(w)  +  dR 


Y^  (u-ti)  «  AR  F^(w)  +  d1 


The  peak  of  the  Imaginary  mobility  occurs  when 

,l  m  r  _  .R  d 

d (to  ) 


dY 

dT?) 


(A-66) 

(A-67) 


(A-68) 


or 


aR  2/Q2  g  1  -  w2/n2 

Tf  7 

■  I  *  _  V  I 

+  g 


1  -  (//Q* 

.  . 


2  2 

which  is  whon  ur/Sr  «  1,  as  stated  previously. 

The  peaks  of  the  real  displacement  mobility  occur  when 


(A-69) 


i 


dY  8  n  -  flR  d 
d(uT )  d(wz) 


FR(w) 


or 


1  -  (o2/J22 


[l  -  0)W 

which  gives  peaks  at 


=  0 


(A-70) 


(A-71) 


(A-74) 


or 


When  the  mode  is  complex, equations  (A-64)  and  (A-65)  indicate  that  both  the 
real  and  imaginary  parts  of  the  displacement  mobility  contain  linear  com- 
binations  of  F  (<d),  F  (co),  oo  and  constants.  The  peaks  in  the  mobilities 
in  the  general  case  may  not  be  simply  related  to  the  natural  frequency  and 
damping  coefficients.  Naturally,  different  degrees  of  approximations  are 
feasible,  depending  on  how  complicated  the  situation  really  is. 

A  general  technique  that  has  been  found  to  be  applicable  to  the  majority 
of  modes  encountered  on  the  AH-1G  helicopter  is  based  on  the  rate  of  change 

of  the  arc  length  of  the  modal  curve  (plotted  on  the  complex  plane,  i.e., 

I  R 

the  plot  of  the  Y  against  Y  with  frequency  as  a  parameter). 


ds 

d(io2) 


where  s  is  the  arc  length. 


(A-75) 


The  rate  of  change  of  the  arc  length  with  respect  to  the  square  of 
frequency  is  stationary  when 


(A-76) 


(A-77) 
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By  substituting  equations  (A-64)  and  (A-65)  into  equation  (A-77)  and  sim- 
plifying,  the  following  condition  for  the  peak  of  the  —  ■  y,  plot  is 
obtained: 


d 

Od 

lx. 

-0 

1 _ 

2 

+ 

r  dF>  ■ 

2 

+  c» 

A R  d2FR 

A1  *‘1 

d(u>2) 

d(o>2) 

CM 

3 

TJ  , 

'  +  T  ■ 

A  2 

d(u2> 

-  A  0 

d(J) 

where 


R  2  I  2 

(AR)  +  (A1) 


,R  dV 


d(w2)‘ 


+  A 


I  d2FR 


d(u,2)‘ 


=  0 


(A-78) 
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For  a  well  separated  mode,  the  constants  C  and  C  will  be  nearly  zero  and 
equation  (A-78)  gives 


d 

d(u2) 


dF 


d(»2) 


( A-79 ) 


Equation  (A-79)  does  not  involve  the  modal  acceleration  coefficients  A  and 
A1.  Thus,  the  condition  for  the  peaking  of  the  rate  of  change  of  the  arc 
length  with  respect  to  frequency  squared  holds  true  regardless  of  how 
complex  the  mode  is,  as  long  as  it  is  well  separated. 


Equation  (A-79)  can  be  expanded,  making  use  of  equations  (A-61 )  and  (A-62), 
and  the  result  is 


d 

d(w2) 


f  2' 

1  -  — 
2 

2 

-  g2 

2 

+  4g2 

2 

1  - 

2 

_ 

l  n  . 

n 

=  0 


(A-80) 
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Equation  (A-80)  can  be  evaluated  to  yield 


d(»2) 


1  w  ,  _< 
i  — p  +  g 

n 


.  4 

*7* 


1  -  %  +9' 


(A-81 ) 


Thus,  for  a  well  separated  mode,  the  peak  of  the  — plot  will  occur 

d(«Z) 

at  the  natural  frequency,  regardless  of  whether  the  mode  is  complex  or 
classical.  Any  suitable  finite  difference  scheme  can  be  used  to  compute 

--dSo-  from  measured  data  using  equation  (A-75). 
d(u>  ) 

Even  for  modes  that  are  not  well  separated,  the  peaks  of  the  -ds9 -  plot 

d(w  ) 

still  give  good  approximations  to  the  natural  frequencies.  To  establish 
why  this  is  so,  consider  equation  (A-78)  term  by  term;  the  first  term 
vanishes  at  the  natural  frequency,  as  we  have  already  seen,  and  the 
remaining  terms  can  be  rearranged  as 


_L  d2FR  (cRaR  +  CIAI]  +  1_  d2F* 

fl2  d(.2)2  1  J  “2  d(A2 


cIaR  _  Q Ra 1 


(A-82) 


At  the  natural  frequency 


d(ui2)' 


(A-83) 


the  remaining  term  becomes 


Lightly  damped  modes  generally  tend  to  be  classical  and  well  separated. 

Tt,is  is  understandable,  since  in  the  limit  of  zero  damping,  a  classical 

8  3 

undamped  mode  results.  Thus,  the  low  damping  that  drives  2/n  g  up, 
also  drives  [cJAR  -  c'V)  down,  effectively  neutralizing  the  expres¬ 
sion.  This  consequently  reduces  the  error  incurred  by  approximating 
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p 

the  natural  frequency  by  the  peak  of  the  ds/d(w  )  plot.  Experience  has 

2 

shown  that  modes  which  are  too  close  to  be  resolved  by  the  ds/d(w  )  routine, 
may  probably  not  be  resolvable  by  any  other  method  presently  known. 

The  diameter  of  the  modal  circle  that  fits  the  curvature  of  the  displace¬ 
ment  mobility  plot  {on  the  polar  plane)  in  the  vicinity  of  the  natural 
frequency  is 

D  =  (A-84) 

gsr 

At  the  natural  frequency 

dY? 

d(u>  )  a)  =  ft 
and 

dY1 

d  (o)^ )  to  =  ft 

Substituting  equations  (A-85)  and  (A-86)  into  equation  (A-75)  gives 

<*—  *—  =  +  Error  of  approximation  (A-87) 

d(w^)  a)  =  ft  g 

2 

The  plot  of  arc  length  s  against  w  has  a  characteristic  S  shape,  as 
shown  by  the  sketch  below 


S 


g^i?" 


CR 

7 


g2n4 


Q 


(A-85) 


(A-86) 


206 


By  fitting  the  best  straight  line  to  the  Inflection  region  of  the  S  curve, 
one  obtains  a,  the  gradient  of  this  line. 

a  “  |A|/g2ft4  (A-88) 

From  equations  (A-84)  and  (A-88),  the  damping  coefficient  Is  evaluated 
as 


g  »  D/aft2  (A-89) 

where  D  is  the  diameter  of  the  circle  fit  to  mobility  data,  a  is  the  grad¬ 
ient  of  the  line  fit  to  the  S  plot,  and  ft  is  the  natural  frequency. 


TESTING  FOR  ORTHONORMAL  MODES  AND  MODESHAPES 

The  mode  shapes  of  any  structure  are  related  to  the  modal  acceleration 
coefficients  as  shown  in  equation  (A-28);  i.e., 

Ajkn  =  m^  ^jn  ^kn  (A 

where  <}>jn  and  <j>kn  are  mode  shape  elements  at  the  jth  and  kth  coordinates 
of  the  nth  mode,  mn  is  the  generalized  mass  of  the  nth  mode. 

The  orthonormal  mode  elements  are  defined  as 


-L 

S/S  J" 

(A-91) 

and 

J-  ♦. 

^  kn 

(A-92) 

Thus, 

Ajkn  =  ^jn  ^kn 

(A-93) 

and 

Ajjn  =  (^jn)Z 

(A-94) 
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It  follows  from  equations  (A-93)  and  ( A-94 )  that 
*kn 

p  JJ" 

It  can  also  be  deduced  that 


(A-96) 


The  choice  of  using  either  equation  (A-95)  or  equation  (A-96)  to  deter¬ 
mine  ip.  depends  on  the  accessibility  of  the  modal  acceleration  coef¬ 
ficients  involved.  Note  that  two  shaking  stations  are  involved  in 

equation  (A-96),  whereas  only  one  shaking  station  is  involved  in 

equation  (A-95).  It  may  turn  out  that  the  driving-point  data  that  yields 

A..„  are  such  that  accurate  estimations  of  the  A.,  for  a  number  of  the 
jjn  jjn 

modes  are  not  easy.  This  may  be  due  in  part  to  a  strong  local  mode 
coupling,  or  a  residual  effect.  In  cases  where  this  is  so,  it  may  be 
better  to  shake  at  a  number  of  coordinates  and  then  use  schemes  similar 
to  that  in  equation  (A-96). 

Consistency  of  the  phase  angle  in  equation  (A-96)  is  achieved  in  the 
following  manner.  For  an  orthonormal  mode  element,  in  the  nth  mode,  of 
large  magnitude,  say  t^,  let 


where  <j>  is  phase  angle.  For  any  other  orthonormal  mode  element,  say  p. 


Mode  shapes  of  the  structure,  normalized  with  respect  to  the  highest 
mode  element,  can  be  obtained  directly  from  the  modal  acceleration 


coefficients  as 


1 


j ,max,n 


A 


\jln 

j2n 


AjNn 


(A-99) 


where  Aj>max  n  is  the  modal  acceleration  coefficient  with  the  maximum 
amplitude  in  the  column  corresponding  to  the  nth  mode,  when  shaking  at 
the  jth  coordinate.  The  generalized  mass,  corresponding  to  the  mode 
shape  thus  normalized  is  computed  from  equation  (90)  as 


mn  ^jn  ^kn^jkn  ~  Ajjn^Aj ,max,n^ 


(A-100) 


In  order  to  obtain  the  elements  of  the  orthonormal  modes  and  mode  shapes, 
the  modal  acceleration  coefficients  of  all  the  modes  for  the  mobilities 
relating  the  response  coordinates  to  the  shaking  coordinates  have  to  be 
determined.  The  computational  scheme  for  determining  the  modal  accelera¬ 
tion  coefficients  requires  mobility  data  at  discrete  frequencies.  The 
technique,  herein  referred  to  as  the  matrix  difference  method,  well  suited 
to  processing  large  numbers  of  transducers  for  modal  analysis  using  multi¬ 
plexing  data  acquisition  systems  common  in  the  helicopter  industry.  The 
natural  frequencies  and  modal  damping  must  be  determined  beforehand. 


For  two  frequencies  and  idT  in  the  region  of  the  natural  frequency  of 

the  ith  mode,  equation  (A-34)  could  be  written  thus: 


Vjk  '  V"?>  '  V“l>  ■  j,  Ajk„Vn 


(A-101 ) 


where 


A. F  =  F  (toT)  -  F  (u>7) 
in  nv  l '  nv  i ' 


Equation  (m-101)  can  be  written  for  all  the  remaining  modes,  having 
selected  the  corresponding  pairs  of  frequencies.  The  resulting  system  of 
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equations  is  the  matrix  difference  equation? 


’  ••  ' 

"*ipi 

¥2  •  • 

A1FN 

Ajkl 

a2  Yjk 

.  = 

Vi 

•  • 

a2f 2  •  • 

a2FN 

« 

Ajk2 

• 

• 

• 

• 

• 

• 

• 

an  \jk 

hf 

•  • 

Anf2 

** 

anfn 

« 

Aj  kN 

,  . 

or 

Hk}  - 

[AF]  • 

From  which  j 

*4  ■ 

[AF]'1 

(A-1G2) 


(A-103) 

(A-104) 


An  immediate  observation  about  the  matrix  difference  scheme  is  that  all 
contributions  to  the  mobilities  near  a  given  mode  which  are  weakly  vary¬ 
ing  with  frequency,  such  as  the  effects  of  distant  modes  or  rigid  body 
modes,  are  subtracted  out.  Proper  selection  of  wj  and  w7,A..Fn  can 
be  made  such  that  A1-Fi  is  large  and  A^Fj  is  small  for  all  j^i.  Experience 
shows  that 


and 


<4  =  V1  +  g1/2^ 

o)T  =  n.(l  -  g./2) 


(A-105) 

(A-106) 


are  the  most  effective  choices  for  the  upper  and  lower  discrete  frequencies. 
For  these  discrete  frequencies,  the  matrix  [AF]  is  well  conditioned  for  in¬ 
version  since  the  off-diagonal  terms  are  small  compared  to  the  diagonal 
terms. 
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Test  procedure  -  Figure  A-11  shows  the  schematic  of  the  Instrumentation 
setup  for  the  shake  test  for  orthonormal  modes  and  mode  shapes.  Signals 
from  all  the  accelerometers  and  from  the  fo m  gage  a*"e  transmitted  via 
telemetry  to  a  computer  where  the  transfer  functions  between  the  response 
coordinates  and  the  force  coordinate  are  computed  and  printed  out.  The 
excitation  signals  are  sinusoidal  at  the  discrete  frequencies  and  wj’ 
for  1  *  1,  2,  ...  N.  The  same  force  levels  used  for  the  swept  sine 
global  parameter  shake  test  are  also  used  for  the  modal  shake  test  at  the 
corresponding  discrete  frequencies. 

DERIVATION  OF  MOBILITIES 

Underlying  any  technique  of  modal  analysis  Is  the  principle  of  linear  de¬ 
composition  of  structural  response  mobility  Into  contributions  from  the 
natural  modes  occurring  between  a  chosen  frequency  interval.  The  preced¬ 
ing  methods  estimate  not  only  the  natural  frequencies  and  damping  coeffi¬ 
cients  of  each  mode,  but  also  the  modal  acceleration  coefficients  of  each 
modal  contribution  to  the  mobility  between  response  and  forcing  coordinates. 

Subsequent  to  the  determination  of  the  modal  parameters  and  modal  constants, 
the  next  logical  step  Is  to  reconstruct  mobilities  both  between  a  pair  of 
forcing  and  response  coordinates  over  a  continuous  frequency  interval,  and 
at  a  chosen  frequency  between  several  pairs  of  forcing  and  response  coor¬ 
dinates.  By  comparing  the  mot  lity  derived  over  a  continuous  frequency 
range  with  the  measured  mobility  over  the  same  frequency  range,  some  assess¬ 
ment  of  the  accuracy  of  the  global  parameter  estimations  can  be  made.  The 
comparison  of  discrete  frequency  mobilities  for  a  large  number  of  coordin¬ 
ate  pairs  allows  the  assessment  of  the  acceptability  of  the  orthonormal 
mode  and  mode  shape  calculations.  The  results  of  these  comparisons  build 
the  confidence  in  the  mobilities  which  are  derived  but  not  actually  measured. 
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Figure  A-ll.  Schematic  of  test  setup  for  matrix  difference 
method  of  modal  testing. 
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Comparison  of  measured  and  simulated  mobilities  over  frequency  band  - 
Global  parameters  and  gn  of  system  modes  occurring  within  a  specified 
frequency  range  can  be  satisfactorily  estimated  using  methods  based  on  the 
properties  of  the  mode  functions,  F  (u>).  The  matrix  difference  method  can 

D  T 

then  be  used  to  calculate  the  modal  acceleration  coefficients  (A^,  Ajkn) 
of  the  relevant  elastic  modes.  Table  A-2  summarizes  the  parameters  esti¬ 
mated  between  0  and  SO  Hz  from  the  tall  vertical  shake/nose  vertical 
acceleration  data.  Figure  A-12  shows  plots  of  the  mobility  measured 
between  2  and  50  Hz.  Using  the  parameters  of  Table  A-2  and  equation  (A-31 ) 
without  Including  the  rigid  body  coefficients,  the  plots  of  Figures  A-13, 
A-14,  and  A- 15  were  generated. 


The  computed  and  measured  mobilities  are  superimposed  In  Figures  A-14  and 
A-15.  It  is  seen  that  the  two  plots  agree  to  within  a  frequency  indepen¬ 
dent  complex  constant,  which  Is  an  estimate  of  the  contribution  of  the 
rigid  body  modes. 


TABLE  A-2.  ESTIMATED  MODAL  PARAMETERS  BETWEEN 
0  -  50  Hz  (TAIL  VERTICAL  SHAKE,  NOSE 
VERTICAL  ACCELERATION) 


Mode  no. 
n 

Natural 

frequency 

Damping 

coefficient 

Nose/tall  modal  acceleration 
coefficient 

Q 

n 

(Hz) 

% 

Real 

ar 

mZN,  ZT,  n 
(g/ib) 

Imaginary 

AZN,  ZT,  n 
(g/lb) 

1 

7.33 

0.062 

7.26  x 

10‘4 

1.40  x  10"4 

2 

8,09 

0.12 

4.48  x 

10‘4 

-3.55  x  10~4 

3 

13.3 

0.13 

-1.57  x 

10'5 

-4.97  x  10'6 

4 

15.97 

0.085 

2.25  x 

10"5 

1.39  x  10”4 

5 

16.35 

0.05 

5.18  x 

10-5 

-1.04  x  10'6 

6 

17.63 

0.08 

1.37  x 

10"4 

-5.54  x  IQ'5 

7 

22.1 

0.15 

-7.28  x 

10'5 

-4.06  x  10“4 

8 

28.4 

0.11 

-1.06  x 

10‘4 

1.13  x  10'4 

9 

40.7 

0.12 

7.84  x 

10-6 

-2.16  x  10"5 

10 

45.3 

0.026 

4.02  x 

10"5 

1.42  x  10'4 

Figure  A-12.  Measured  acceleration  mobility 
between  2  Hz  and  200  Hz. 
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Figure  A-13.  Numerical  simulation  of  the  elastic  component  of  the 
acceleration  mobility  between  0  Hz  and  50  Hz. 


Figure  A-14.  Simulated  and  measured  real  part  of 
acceleration  mobility  2  -  50  Hz. 


Figure  A-15.  Simulated  and  measured  imaginary  part  of 
acceleration  mobility  2  -  50  Hz. 


Modal  Series  Method  -  Equation  (A-31 )  can  be  rewritten  in  the  following 
form:  N 

»jk  *  EJk  *  RJk(“>  +  AJkn  En<“>  +  Rjk<“)  'A-'07> 

where  , 

Rjl^(w)  is  the  low-frequency  mobility  residual;  i.e.,  contribu¬ 
tions  to  the  mobility  by  elastic  modes  which  occur  at 
frequency  below  the  lower  test  frequency  limit. 

U 

Rjk(w)  is  the  high-frequency  mobility  residual. 

The  rigid  body  acceleration  coefficient,  E^,  is  determined  from  geometry 
and  weights  data. 

If  the  lower  test  frequency  limit  is  near  zero,  it  follows  that  RL  =  0. 

The  higher  test  frequency  limit  is  usually  selected  high  enough  so  that  RH 
can  be  safely  assumed  to  vanish  for  all  but  certain  driving-point  mobili¬ 
ties  wh\ch  may  suffer  either  from  local  mode  effects  or  from  high-frequency 
mode  residuals. 


When  the  natural  frequencies,  the  damping  coefficients,  and  the  modal  accel¬ 
eration  coefficients  have  been  determined,  the  acceleration  mobilities  be¬ 


tween  pairs  of  motion  coordinates  which  uo  not  include  the  shaking  coor¬ 
dinate  can  be  computed  from 


A 


£kAkn 

Akkn 


Fn(w) 


(A-108) 


where  k  is  the  coordinate  of  the  shaking  station  for  the  data  which  gener¬ 
ated  A£)<  and  A^.  It  is  necessary  to  select  the  shaking  station  k  such 
that  there  is  no  local  mode  or  high-frequency  residual  effect  cn  the  esti¬ 
mated  value  of  A^kn. 


If  only  of  the  modes  are  well  defined  by  shaking  at  k,  while  the  remain¬ 
ing  N  modes  are  better  defined  by  the  shake  at  p,  then 


Y  =  E  + 
TSjti 
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n=l 


^ftkn^mkn 

Akkn 
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p  AtpnAmi 
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(A-109) 
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